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Abstract 

^ Shadowing is a quantum phenomenon leading to a non-additivity of electroweak cross sections 

on nucleons bound in a nucleus. It occurs due to destructive interference of amplitudes on different 
nucleons. Although the current experimental evidence for shadowing is dominated by charged- 
lepton nucleus scattering, studies of neutrino nucleus scattering have recently begun and revealed 
unexpected results. 
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Part I 

Theory of Nuclear Shadowing 

1 Introduction 

The term shadowing was naturally incorporated into quantum mechanics as a tag for one of the basic 
quantum phenomena related to the destructive interference between elastic amplitudes on different 
scattering centers. The popular observable sensitive to shadowing in a process, is the normalized ratio 
of the cross sections of this process on nuclear and nucleon targets, 

Shadowing manifests itself in a suppressed value of this ratio, Ra/n < 1, i-e. the nuclear cross section is 
less than the sum of the cross sections on all bound nucleons^J Intuitively, the source of the suppression 
is the survival probability for the projectile particle or its fluctuations to pass through the nuclear 
medium and reach a bound nucleon deep inside the nucleus. In other words, the interaction with a 
given bound nucleon is shaded by the probability of having preceding interactions with other nucleons. 
This is easy to see on the example of the total hadron-nucleus interaction cross section. It is clear that 
a strongly interacting particle has no chance to pass through a heavy nucleus without interaction. This 
means that the partial elastic amplitude at impact parameter less than the nuclear radius, b < Ra, 
nearly saturates unitarity bound, lmf^ A (b) ^ 1, i.e. the of f ~ 2nR A , or Ra/n oc A^ 1 ^ 3 , because the 
nuclear radius Ra oc A 1 / 3 . In the following sections we calculate these quantities more accurately within 
the Glauber-Gribov approach. 

hereafter we neglect the isospin corrections, i.e. do not make difference between target protons and neutrons, unless 
specified. We also neglect the effect of binding, unless it is important. 
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Notice that in exclusive reactions one should discriminate between shadowing and final state ab- 
sorption. Both effects lead to a suppression of the ratio (IT]), but have quite different origins. 



2 Shadowing in soft interactions 
2.1 Shadowing in the Glauber model 

The Glauber model [1] was the first theoretical approach that correctly calculated the effects of shad- 
owing in hadron-nucleus interactions. The probability for a hadron to interact with the nucleus is one 
minus the probability of no interaction with any of the bound nucleons. So the HA elastic amplitude 
at impact parameter b has the eikonal form, 



T hA {b- {s v z 3 }) = 1 - J] [l " r ftJV (& - s k ) 

k=l 



(2) 



where {sj, Zj} denote the coordinates of the target nucleon Nj. iT hN is the elastic scattering amplitude 
on a nucleon normalized as, 
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In the approximation of single particle nuclear density one can calculate a matrix element between 
the nuclear ground states. 



T hA (b;{s J}Zj }) 



I __ i d 2 sT hN (s) I dzp A (b-s } z) 
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J i=2 



(4) 



(5) 



is the nuclear single particle density. 

The magnitude of shadowing depends on the process. Here we calculate the shadowing effects within 
the Glauber model for several basic processes. 

Total cross section. 

The result Eq. Q is related via unitarity to the total hA cross section, 
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where p pp is the ratio of the real to imaginary parts of the forward pp elastic amplitude; 
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(6) 
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and 

/oo 
dzp A (b,z) , 
-oo 

is the nuclear thickness function. We use exponential form of T hN (s) throughout the paper, 

ReT hN (s) 
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(9) 



where is the slope of the differential hN elastic cross section. Note that the accuracy of the 
optical approximation (the second line in (pj) is quite high for gold, ~ 10~ 3 . We use the optical form 
throughout the paper for the sake of simplicity, although for numerical evaluations always rely on the 
accurate expression (the first line in (|6|). The effective nuclear thickness Eq. ([7]), implicitly contains 
energy dependence, which is extremely weak. 

In what follows we also neglect the real part of the elastic amplitude, unless specified, since it gives 
a vanishing correction ~ p^ p /A 2 ^. 

In Table jlj the result of the calculation of the shadowing factor R^/n f° r t ne total proton-lead cross 
section at y/s = 5.5 TeV is shown. We use a™ = 93.1 mb and elastic slope B^ N = 20GeV -2 , which 



Table 1: The shadowing suppression factor Eq. Q for the total, elastic, inelastic, quasielas- 
tic and production cross sections. 



Model 


TDtot 

^A/N 


TDel 

^A/N 


pin 

A/N 


A/N 


Glauber 


0.22 


0.4 


0.17 


0.03 



agree well with the recent measurements of the TOTEM experiment [2] at y's = 7 TeV. We see how 
strong the shadowing effect for the total cross section is. 



Elastic cross section. 

As far as the partial elastic amplitude is known, the elastic cross section reads, 

(10) 



JiA 



d 2 b 



From the numerical result depicted in Table [T] we see that the elastic cross section is less suppressed by 
shadowing than the total one. This shows that even at the energy of LHC the proton-lead interaction 
does not reach yet the unitarity bound, where the shadowing ratios Eq. ([T]) are expected to be alike for 
both processes. 

Total inelastic cross section. 

is given by the difference between the total and elastic cross sections, 

* = o% - a h / = Jd%{\- exp [-o™ T h A (b)] } . (11) 

This cross section covers all inelastic channels, where either the hadron or the nucleus (or both) break 
up. According to Table [T] inelastic interactions are shadowed stronger than the total cross section. 
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Quasielastic cross section. 

As a result of the collision the nucleus can be excited to a state \F). Summing over final states of the 
nucleus and applying the condition of completeness, one gets the quasielastic cross section, 
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(12) 



Here we extracted the cross section of elastic scattering when the nucleus remains intact. 
Then in the fist term of this expression we make use of the relation, 
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and arrive at, 
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(13) 



(14) 



The large "volume" terms proportional to A cancel in this expression, so one should expect a very strong 
nuclear suppression, as is confirmed by the result of evaluation presented in Table [T] Strictly speaking, 
this suppression is due to a combination of shadowing in the initial state and the survival probability 
of the scattered hadron in the nuclear matter. Final state interactions are particularly important for 
exclusive channels of hadron production. In the case of electro-weak interactions final state interactions 
can imitate shadowing even when it is absent, e.g. when the coherence time is short. 



Production cross section. 



Quasielastic scattering is a part of the inelastic cross section Eq. (11), and the only process with no 



production of new particles. Therefore it should be subtracted to get the total production cross section. 
The result is rather simple, 
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2.2 Inelastic shadowing corrections 
2.2.1 Intermediate state diffractive excitations 

The Glauber model is a single-channel approximation, therefore it misses the possibility of diffractive 
excitation of the projectile in the intermediate state as is illustrated in Fig. [TJ These corrections called 
inelastic shadowing were introduced by Gribov back in 1969 [I]. The formula for the inelastic corrections 
to the total hadron-nucleus cross section was suggested in [SJ, 
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(16) 



where is the cross section of single diffractive dissociation hN — > XN with longitudinal momentum 
transfer 



M 2 - m\ 
2E h 



(17) 



This correction makes shadowing stronger and nuclei more transparent [6] , because it is added with 
positive sign to the probability amplitude of having no interaction, exp(— It takes care of 
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Figure 1: Diagonal and off-diagonal diffractive 
multiple excitations of the projectile hadron in 
intermediate state 



Figure 2: Data and calculations [3] for the total 
neutron-lead cross section. The dashed and solid 
curves correspond to the Glauber model and cor- 
rected for Gribov shadowing respectively. 



the onset of inelastic shadowing via phase shifts controlled by and does a good job describing data 
at low energies [3J E] , as one can also see in Fig. [2j Notice that the higher order off-diagonal transitions 
neglected in ( Jl6| ), including the diagonal transitions (or absorption of the excited state), are important, 
but unknown. Indeed, the intermediate state X has definite mass M, but no definite cross section. It 



was fixed in (|16|) at with no justification. We will come back to this problem later. 



There is, however, one case which is free of these problems, shadowing in hadron-deuteron interac- 
tions. In this case no interaction in the intermediate state is possible, and knowledge of diffractive cross 
section hN — > XN is sufficient for calculations of the inelastic correction with no further assumptions. 



In this case Eq. (16) takes the simple form [U [8], 

MN 



2.2.2 Inelastic shadowing in the eigenstate representation 

If a hadron were an eigenstate of interaction, it could experience only elastic scattering (as a shadow of 
inelastic channels) and no diffractive excitation were possible. In this case the Glauber formula would 
be exact and no inelastic shadowing corrections were needed. This simple observation suggests to switch 
from the basis of physical hadronic states, which are the eigenstates of the mass operator, to the basis of 
a complete set of mutually orthogonal states which are eigenstates of the scattering amplitude operator. 
This was the driving idea of description of diffraction in terms of elastic amplitudes [HI [TO], and becomes 
a powerful tool for calculation of inelastic shadowing corrections in all orders of multiple interactions 
[llj . Hadronic states (including leptons and photons) can be decomposed into a complete set of such 
eigenstates \k), 

l*> = X>Sl*>, ( 19 ) 

k 

where are hadronic wave functions in the form of Fock state decomposition. They obey the orthog- 
onality conditions, 

£ ^ = Shhl ; 

k 
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(20) 



We denote by = io\^ t /2 the eigenvalues of the elastic amplitude operator / neglecting its real 
part. We assume that the amplitude is integrated over impact parameter, i.e. that the forward elastic 
amplitude is normalized as l/^l 2 = da^ /dt\ t=Q . We can express the elastic f e i and off diagonal 
single-diffractive f sc t amplitudes as, 
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hN 



kN 
tot 



(21) 
(22) 



Notice that if all the eigen amplitudes were equal, the diffractive amplitude (22) would vanish due 
to the orthogonality relation, (20). The physical reason is obvious. If all the f^ N are equal, the 



interaction does not affect the coherence between the different eigen components \k) of the projectile 
hadron \h). Therefore, the off-diagonal transitions are possible only due to the differences between the 
eigen amplitudes. 



Summing all final states and making use of the completeness condition (20), then, excluding the 
elastic channels one arrives at [TTJ [T21 H3] , 
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(23) 



If the lifetimes of the different eigenstate components of the hadron are sufficiently long, so that 
they do not mix with each other during propagation through the nucleus, the cross sections for different 
processes on nuclei can be written as, 
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(24) 

(25) 
(26) 



Notice that the last expression for cr^ 4 is already free from the diffraction contribution. Although 
only elastic and quasielastic cross sections were subtracted from in the Glauber model, Eq. (15), 
after averaging over eigenstates it turns out that diffraction is subtracted as well [Hj. 

The difference between the cross section Eq. (24) and the Glauber approximation Eq. ([6]), is in the 
way of averaging. In the former case the whole exponential is averaged, while in the Glauber approx- 
imation only the exponent is averaged. The difference should correspond to the Gribov corrections 
summed in all orders. Indeed, the first order terms in the expansion of (24) and ^ cancel and in the 
second order using the relation (23) we get, 



hA 
tot 



dot 



d 2 b - [(a^) 2 - ((^) 2 >J T» 2 = -4vr J d%T h A {bf J dM 2 ^L. ^ . ( 27 ) 

This result is identical to Eq. (16), if we neglect there the phase shift vanishing at high energies, and 
also expand the exponential. 



Note that since the inelastic nuclear cross section in the form Eq. (15) is correct for eigenstates, one 



may think that averaging this expression would give the correct answer. However, such a procedure 
includes the possibility of excitation of the projectile and disintegration of the nucleus to nucleons, but 
misses the possibility of diffractive excitation of bound nucleons which is not a small correction. We 
introduce a corresponding correction in the next section. 
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2.2.3 Dipole description of shadowing 

The light-cone dipole representation was proposed in [13J as an effective tool for calculation of hadronic 
cross sections and nuclear shadowing, relying on the observation that color dipoles are the eigenstates 
of hadronic interactions at high energies, and the eigenstate method [H] allows to sum up the Gribov 
inelastic corrections to all orders. 

The key ingredient of this approach is the cross section of the dipole-nucleon interaction, ^(r^, s), 
which is an universal and flavor independent function, depending on transverse separation tt and 
energy. One cannot calculate it reliably, since that would involve nonperturbative effects, but should 
fit it to data. Once this cross section is known from data on a proton target, the nuclear effects can 
be predicted. Still, the results of calculations remain model dependent, because the hadronic wave 



functions participating in averaging Eq. (23) are poorly known. 

Applications of the dipole formalism to nuclei are especially simple, if the energy is sufficiently high 
to freeze the fluctuations of the dipole size during its propagation through the nucleus. Otherwise one 
should rely on the path-integral technique [15], [T6J [17] , which takes care of these fluctuations (see below, 



Sect. |3.2.2[ ). 

Due to color screening a colorless point-like dipole cannot interact with an external color field. Since 
the underlying theory is non-abelian, the interaction cross section for such dipoles vanishes at — > as 
^qqi^r) oc 7"t [IS], a phenomenon called color transparency]^] At high energies nuclei are transparent for 
small-size fluctuations of the incoming hadron, therefore the strong exponential attenuation suggested 
by the eikonal Glauber formula cannot be correct and the nuclear medium should be more transparent, 



as was already mentioned above. For a dipole of a fixed size the eikonal form of Eqs. (24)-(|26j) is 
exact. In this case the role of the eigenvalues of the cross section a\^ t is played by the dipole one (Jq q (r T ), 
and the total hadron-nucleus cross section has the form H31. 



a% = 2 I d 2 b I d 2 r T \V h (r T 



■±<Tq q (r T )T A (b) 



(2* 



It is interesting to point out that although nuclear transparency, which is the probability of no- 
interaction, exponentially falls with nuclear thickness, after averaging over dipole sizes for a long path 
in the nuclear medium the resulting transparency drastically changes p2] , 

e -* ?q (r T )T A ^ f ^ |^( rT )|2 e -<x qq (r T )T A ^ J_ _ ^ 

J T A 

Here we assumed a Gaussian shape for the /^-dependence of the hadronic wave function \E , / l (rr), and 
small-rY regime for the dipole cross section, <Jq q {rT) oc rf,, which is justified for a large nuclear thickness 
T A . 



The result (29) should be compared with the exponential attenuation in the Glauber model, given 
by exp(— T4). Apparently, the difference cannot originate from the lowest order inelastic correction 



Eq. (16), which has the same exponential dependence on T4. This color transparency effect includes 
all the higher order inelastic corrections. 

The Gribov inelastic shadowing corrections were calculated in [T4"l ITSf [H?| |2T)] for various soft inter- 
action processes. Although the nuclear transparency significantly increases for heavy nuclei, the related 
variation of shadowing is pretty mild. This happens because for heavy nuclei the transparency term 
in the cross section (the second term in Eq. ([6])) is a small correction to the big first (volume) term. 
Even a considerable variation of this small correction does not affect much the total cross section, i.e. 
shadowing. Indeed, the Gribov correction evaluated in Fig. [2] at low energy is several percent. With 
rising energy the hadron-nucleus amplitude approaches the black-disc limit (unitarity bound), where 

2 Actually, the cross section behaves as r\ ln(ry) |13| . but with a good accuracy one can fix the logarithm at an effective 
value of tt typical for the process under consideration. 
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Gribov corrections vanish, because all eigen amplitudes become equal. We found that at the energy of 
LHC the Gribov corrections to the numbers in Table [T] are only minus one percent. So we skip this 
comparison here, but one can find the details of the calculations in [TU [18J fT9] . 

2.3 Shadowing in photo-nuclear reactions 

At first glance the weakly interacting particles, like photons or leptons, should interact with nuclei with 
no shadowing. Indeed, applying the Glauber model expression ^ to the photoabsorption cross section 
cr^of one gets a vanishingly small shadowing effect. However, data presented in Fig. [3] clearly show that 
the photoabsorption cross section is strongly shadowed as function of energy. The nuclear ratio 
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Figure 3: Nuclear ratio R\/ N for lead as function of photon energy. The curves correspond 
to a fit and model calculations [2"Tj . 
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(30) 



is significantly suppressed by nuclear effects, which should be interpreted as shadowing. 

The mechanisms of shadowing for photon interactions are well explained in the comprehensive review 
|22j . The observed shadowing suppression is a manifestation of the hadronic properties of the photon. 
Namely, the photon interacts with hadrons via its hadronic fluctuations. The observed smallness of the 
photoabsorption cross section is related to the smallness of the fluctuation probability (it is proportional 
to a em ), while the hadronic cross section is large. Although the fluctuation life time t ~ 1/rrih is quite 
short, it is subject to Lorentz time dilation, and at sufficiently high energies it maybe much longer 
than the nuclear size, as is illustrated in Fig. |4j In this regime the photon-nucleus interaction is indeed 
shadowed as much as in a hadron-nucleus collision. The transition region between the short and long 
regimes can be described adding up the two contributions from the single- and double-step interactions, 
as is illustrated in the left part of Fig. |4} In the case of vector meson dominance, i.e. dominance of one 
pole in the dispersion relation in Q 2 for the amplitude, the sum of the two contributions to the photon 
nucleus amplitude, pictorially shown in the left part of Fig. Ill reads [22J, 



p7 

A/N 



8tt 



Aa 



■yN 
tot 



p T =0 



d 2 b / dz 1 p A (b,z 1 ) / dz 2 p A (b,z 2 )e^ {z2 - z ' ) e-^ TA{b ' Z2 ' Zl \ (31) 
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Figure 4: The sum of the two possibilities of either elastic *yN scattering on a bound nucleon, 
or with preceding diffractive production of an intermediate state X (left), at high energies 
is equivalent to photon interaction via a hadronic fluctuation X (right). 



where q\ = my/2E 1 ; T^b, Z2,Zi) = dz Pa(P, z). Here we neglected the real parts of the elastic and 
diffractive amplitudes. 

This mechanism of shadowing is a particular case of Gribov corrections, and its evaluation is suffering 
of similar problems. First of all, even in the lowest order correction, the second term in Eq. (31), the 
value of aY f is unknown and is fitted to data. The higher order multiple scattering corrections contain 
off-diagonal diffractive transitions 7 — > X — > X' — > X"..., which are also unknown. These higher 
order corrections become more important with increasing photon virtuality. Indeed, the size of the qq 
fluctuation of the photon decreases as (r^) ~ l/Q 2 and the magnitude of shadowing should diminish 
accordingly. In hadronic representation it does not look that simple and can be achieved only by a 
specific tuning of diagonal and off-diagonal diffractive amplitude, which have to cancel each other. 
This is usually modeled within generalized vector dominance models, but the predictive power of such 
approaches is very low. 

The alternative color dipole approach naturally incorporated the color transparency features, how- 
ever the qq distribution amplitudes cannot be calculated perturbatively at low Q 2 and can only be 
modeled. These amplitudes were calculated within the instanton vacuum model in [23] and the cross 
section of Compton scattering on protons and nuclei was evaluated within the dipole formalism in 



2.4 PCAC and shadowing of soft neutrino interactions 

Neutrinos are known as an significant source for the axial current, which is subject to nuclear shadowing, 
similar to that for the vector current, discussed in the previous section. There are, however, specific 
features of shadowing of the axial current, related to partial conservation of axial current (PCAC), 
which are absent in the case of vector current. Indeed, the vector current is "trivially" conserved, 
oc (m n — m p ) rs 0, since the proton and neutron masses cancel. At the same time conservation 
of the axial current seems to be heavily broken, oc m n + m p . In order to get the axial current 

conserved one has to introduce into additional terms besides 757^, and this is how the massless 
Goldstone pseudo-scalar pion appears. In reality the pion has a small mass, this is why the current is 
conserved partially, 

d,tf = mlf n <p n , (32) 

where and f n ~ 0.93771^ are the pion mass and decay coupling, and is the pion field. 

A beautiful manifestation of PCAC is the Goldberger-Treiman relation [23], which bridges weak and 
strong interactions. It surprisingsly connects the pion decay constant with the pion-nucleon coupling, 
which seem to have very little in common. Indeed, the former depends on the pion wave function, while 
the latter is controlled by the wave function of the nucleon. Nevertheless, data on /3-decay and muon 
capture confirm this relation between very different physical quantities. 



11 



In the case of high energy neutrino interactions PCAC leads to the Adler relation (AR) between the 
cross sections of processes initiated by neutrinos and pions [2"6] . 



d 2 a{vp -+IX) 



dQ 2 dv 



G 2 ,2 E 



2tt 
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Ev 



o(ixp ->■ X). 



(33) 



where E is the neutrino energy; G = 1.166 x 10~ 5 GeV~ 2 is the electro-weak Fermi coupling; Q 2 = —q 2 , 
where = — k'^ and v = E — E' are the 4-momentum and energy transfer in the v — > I transition 
(the same notation as for neutrinos should not cause confusion). 

The Adler relation expresses the axial current interaction with that of the pion, similar to the 
vector dominance model, which relates the interactions of the vector current and p-meson. This make 
it tempting to interpret the Adler relation as pion dominance. It turns out, however, that a neutrino 
cannot fluctuate to a pion, v — > ttI, because the pion pole in the dispersion relation in Q 2 for the axial 
current does not contribute to the interaction of the neutrino at high energies [271 1211 [29] . Indeed, the 
axial current J^(Q 2 ) can be presented as, 



Q 2 + m\ 



T(ixp ->■ X) + 



fa 



Q 2 + m\ 



M^{ap ->X) + 



(34) 



Here the second and following terms represent the contributions of the a\ meson and other heavier 
axial-vector states like p — it and other multi-particle states. 



The first term in (34), corresponding to the pion pole, contains the factor q^, which then terminates 
its contribution to the cross section, Eq. (33). Indeed, the amplitude of the reaction is 



A{vp ->■ I X) oc L^J 



(35) 



where L M = l{k')^^{l + j 5 )u(k) is the lepton current, which is transverse, i.e. q^L^ = (up to the 



lepton mass, assumed to be zero in what follows, unless specified). Therefore, the pion term in (34) 
does not contribute to the amplitude Eq. (35), and this is true at any Q 2 . 



Thus, PCAC does not mean pion pole dominance, but connects the contribution of heavy axial states 



(the second term in Eq. (34)) with the nonexistent pion contribution [23 [2HI EH]- One can see that 
in the Q 2 -dependence of the neutrino cross section at small Q 2 . Fitting the measured Q 2 -dependence 
with the parametrization (Q 2 + M 2 ^)~ 2 , one gets M e ff pa 1.1 GeV [29], which is a clear evidence of the 
dominance of heavy states in the axial current. 



2.4.1 Shadowing in the total neutrino cross section 

Nuclear shadowing of neutrinos was predicted a long time ago [30] based on the simple observation that 
shadowing in the pion-nucleus cross section immediately means that neutrino cross section is shadowed 
as well, because both are connected by the Adler relation. However, this connection may be affected 
by coherence phenomena, which are discussed below. 

As far as the effective mass of a typical hadronic fluctuation of a neutrino, v — > I + a, is large, quite 
a high energy, v J> 10 GeV, is needed to make the fluctuation lifetime, 

2E 

tfluct - (36) 
1V1 eff 

comparable with the radii of heavy nuclei. Then one might jump to the conclusion that there should 
be no shadowing at low energies. 

However, this conclusion is not correct. It is based on the usual wisdom that the fluctuation lifetime 
and the coherence time are equivalent quantities, which is usually correct, but not in this case. The 
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Figure 5: Left: Direct inelastic interaction of the axial current with a bound nucleon. 
Diffractive pion production preceding the inelastic collision. 



Right: 



amplitude of inelastic neutrino-nucleus collision is shown schematically in Fig. [5] The left picture 
corresponds to a direct inelastic interaction of the projectile axial current with a bound nucleon, while 
in the right picture the diffractively produced pion experiences an inelastic collision. This looks similar 



to the calculation of shadowing for the photoabsorption cross section done in Sect. in the case of 
vector dominance. In fact, it seems to be even better justified, since the pion pole is much closer to 
the physical region, than the p pole. However, as we have just learned, it is not the pion pole which 
is behind the Adler relation, but heavier axial-vector singularities, which conspire together and imitate 
the pion pole. Nevertheless, the diffractively produced pion is on mass shell, but the details of the 
dynamics are hidden in the amplitude of diffractive neutrino-production of pion. 

Since both contributions depicted in Fig. [5] lead to the same final state, they can interfere. Squaring 
the inelastic amplitude we get the following expression for the ratio Q for the total neutrino-nucleus 
cross section 



S7T 



m diff 



Arr uN 
n(T tot 



dp\ 
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d 2 b / dz lPA (b, Zl ) / dz 2 p A (b,z 2 )e- iq ^ Z2 - Zl) e~i a ?°t TA{b > Z2 > Zl \ (37) 



was defined in Eq. (31). We neglected here the higher order terms in multiple 



where T A (b,z 2 ,z 1 ) 
interactions. 

At very low energy where (f L is large, the second term in (37) is suppressed, and the first term, which 
corresponds to a nuclear cross section proportional to A, dominates. At high energies g£ 1 can be 
neglected and the integrations over zi t % can be performed analytically |3T]. Due to the Adler relation 
the first term in (37) and the volume part of the second term cancel, and the rest is the "surface" term 
oc A 2 '\ 

It turns out that the high-energy regime actually starts at rather low energies if Q 2 <m 2 . Indeed, 
for elastic neutrino-production of pions, up — > I ttN, the longitudinal momentum transfer, (f L = (m 2 + 
Q 2 )/2u, i.e. the coherence time, 

1 1v 

tl 



mi 



+ Q 2 



(3* 



is very long even at low energy of few hundred MeV. This is actually what matters for the onset of 



shadowing. As for the fluctuation lifetime Eq. (36), it is indeed much shorter. 



This is a result of the nontrivial origin of PCAC. The impossibility for the axial current to fluctuate 
into a pion leads to the dominance of off-diagonal processes, like v — > fiai and a±N — > irN. The 
same happens for the vector current, if one considers, for example, p photoproduction via intermediate 
excitation p'\ 7 — > p' and p'N — > pN. Such an off-diagonal contribution is negligibly small for the 
vector current, but is a dominant one for neutrinos. Only for diagonal transitions the fluctuation 
lifetime and the coherence time are equal, tfi uct = t c . For neutrino interactions the former controls the 
Q 2 dependence of the cross section, while the latter governs shadowing. 
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The results of numerical evaluation of the shadowing effect Eq. (37) for neon are plotted in Fig. [6j 
The calculations [3T] done within the Glauber approximation, Eq. (37), and also including the Gribov's 




Glauber 
Gribov 



0.7 
10 



v (GeV) 



10 



Figure 6: The neon to nucleon ratio of to- 
tal neutrino cross sections at different Q 2 [3H 
32] . Dashed and solid curves correspond to the 
Glauber and Gribov corrected calculations. 




v (GeV) 

Figure 7: The neon to proton ratio of the total 
neutrino cross sections, calculated in [311 132] f° r 
x < 0.2 and Q 2 < 0.2 GeV 2 . The data present 
the BEBC results 



inelastic corrections (important at high energies v > m 2 ai RA) are plotted in Fig. [6] by solid curves as 
function of energy, for different Q 2 . 

The calculated shadowing effects are compared with BEBC data [33] in Fig. [7j As was anticipated, 
the shadowing exposes an early onset, and a significant suppression occurs at small Q 2 in the low energy 
range of hundreds MeV. This is an outstanding feature of the axial current. This seems to be supported 
by data, although with rather poor statistics. 



2.4.2 Diffractive neutrino interactions: Breakdown of PCAC by shadowing 

The process u + p — > I + tt +p offers probably the most stringent test of PCAC in neutrino interactions. 
Indeed, the analysis performed by Piketty and Stodolsky [28J revealed a potential problem related 



to the dispersion representation, Eq. (34). Assuming the dominance of the a\ pole in the dispersion 
relation they arrived at an equation connecting the elastic and diffractive pion-nucleon cross section, 
a(np — > dip) = a(7rp — > Tip). This relation strongly contradicts data: diffractive production of a± meson 
is more than an order of magnitude suppressed compared with the elastic cross section. 

This puzzle was relaxed in [311 EE] by pointing out its shaky point, namely, the a\ pole cannot 
dominate in the axial current, since it is quite a weak singularity compared to the p pole in the 



vector current case. In fact, the main contribution to the expansion Eq. (34) comes from the p-ix 
and 37r cuts, related to diffractive pion excitations. The invariant mass distribution for diffractive pion 
excitations peaks at M w 1.1 GeV [32] and is well explained by the so called Deck mechanism [SB] of 
diffractive excitation tt — > p-n . The interpretation of this peak has been a long standing controversy, 
until a phase-shift amplitude analysis (see references in [37]) eventually revealed the presence of the 
weak a\ resonance having a similar mass. Summing up all diffractive excitations (excluding large 
invariant masses corresponding to the triple-Pomeron term), one concludes that the magnitudes of 
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single-diffractive and elastic pion-proton cross section are indeed similar. This helps to resolve the 
Piketty-Stodolsky puzzle. 

Nevertheless, it was demonstrated that absorptive corrections break the relation between the heavy 



states and pion in the dispersion relation (34) imposed by PC AC. The deviation from the Adler relation 
was estimated in [38j ES EQ] at about 30%. 

Even more dramatic breakdown of PCAC caused by nuclear shadowing was found in [38] for diffrac- 
tive neutrino-production of pions on nuclei. These processes are usually classified as coherent or incoher- 
ent, which according to the conventional terminology correspond to processes which leave the nucleus 
intact or the nucleus breaks up into fragments, respectively. 

In what follows we assume the validity of the Adler relation for a nucleon target, in order to identify 
the net nuclear effects. In the coherent pion production process the amplitudes on different nucleons 
interfere, and the interference is enhanced by the condition that the nucleus remains intact. The 
coherence effects can lead to substantial deviations from the AR and from simplified expectations, as 
is demonstrated below. 



In addition to the pion coherence time, Eq. (38), another time scale related to heavy axial- vector 
state is important, 

^ = ~ = (39) 

The two scales control the amplitude of coherent neutrino-production [38], whose imaginary part has 
the form, 

M vA ^ A (y, Q\ b) = M x {u, Q\ b) - M 2 (z/, Q\ 6), (40) 

where 

oo 

Mi(i/,Q 2 ,6) = M vN ^ N {v,Q 2 ) J dze^ z PA (z,b)e~^T A (b,z). (41) 

— oo 

oo 

M 2 (z/,Q 2 ,6) = M uN ^ laN {v,Q 2 )M aN ^ N {v) [ dz e ^~^ z p A (z,b) e~^ T ^ 



x / dz 1 e l ^ PA {z 1 ,b)e-^o N AT A {b, Zl )-T A ^)]_ (42) 



The structure of the amplitude is similar to the one in Eq. (37). While the first term, M 1; corresponds 
to pion production by the neutrino without any preceding interaction, the second term M 2 corresponds 
to diffractive production of the heavy state a preceding the pion production. This is the first order 
Gribov inelastic shadowing correction [4] to the coherent pion production amplitude. 

The results for coherent production of pions on lead are depicted in Fig. M as function of transferred 
energy v, for several values of Q 2 . The suppression predicted at low energies is related to the shortness 
of t* <C R A and lack of coherence. At higher energies the nuclear ratio at Q 2 = forms a plateau 
in the energy interval 0.5 < v < 5GeV, which corresponds to the condition t£ 3> R A , but t a c R A . 
The height of the plateau corresponds to the Adler relation, which holds in this regime of short l a c . 
Remarkably, at higher energies R^ N (u, Q 2 = 0) steeply falls from the A-dependence A 2 ^ 3 down to A 1 / 3 , 
exposing a dramatic violation of the Adler relation. The source of the breakdown is shadowing, which 
begins at long 

As function of Q 2 the plateau considerably shrinks leaving almost no room for the Adler relation 
(extrapolated to nonzero Q 2 ) to hold, as one can see from the few examples depicted in Fig. [8] 

A similar behavior is expected for incoherent production of pions v + A — > l + ir + A*, except at the 
low energy range where no suppression is expected, because pions are produced on different nucleons 



15 




Figure 8: Nuclear ratio R^ N (u,Q 2 ) of p T - 
integrated cross sections of coherent neutrino- 
production of pions on lead at different Q 2 . 



Figure 9: The same as in Fig. [8] for incoher- 
ent production on lead at Q 2 = 0, 0.2, 0.5 and 
1 GeV 2 . 



incoherently. Therefore the plateau corresponding to the Adler relation starts at very low energies, 
as is demonstrated in Fig. [9] for several values of Q 2 . At higher energies, when t a c < Ra initial state 
shadowing causes an additional suppression, and similarly to the coherent case the A dependence drops 
from A 2 / 3 down to A 1 ^ 3 as is demonstrated in Fig. [9] (the details of calculations can be found in [38]). 
Again, shadowing led to a dramatic breakdown of PCAC. 



3 Shadowing of deeply virtual photons 

The sizable shadowing effects observed above were caused by the large cross sections of soft interactions 
of light hadrons, or soft hadronic fluctuations of photons and leptons at a low virtuality. It has been 
observed, however, that considerable shadowing also suppresses the cross section of hard reactions on 
nuclei, like deep-inelastic lepton scattering (DIS), Drell-Yan reaction, etc. In what follows we demon- 
strate that the dominant contribution to the shadowing effects comes from the soft component of hard 
processes, while the hard component is not shadowed. 

The cross section of deep inelastic lepton scattering I + iV — )■ I' + X has the form, 

*° 47raL / ( n y _^Mi) t^l +y^Q*)), (43) 



dxdQ 2 Q 4 (A Q 2 J x 

where the invariant structure functions can be expressed in terms of the total cross sections <jt{x,Q 2 ) 
and ai{x, Q 2 )oi transversely and longitudinally polarized virtual photons, 

Fi(x,Q 2 ) = Q *}\~ x) (a T + a L ); 



'em 

F^Q 2 ) = f ( / " X) a T . (44) 
8ir 2 a em x 

We use the standard notations, y/s and W are the cm. energies; x = Q 2 /(2qP), y = Q 2 /(2kP), where 
k, q and P are the 4-momenta of the initial lepton, virtual photon and the target nucleon, Q 2 = —q 2 . 
Notice that the nuclear ratio of the structure functions, 

+ L (45) 



TpN ' -y*N . y*N ' 
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which is independent of the lepton energy, is not the same as the measured ratio of the lepton-target 
cross sections, 



M F A{l + ^ L/T ){l + Rf /T ) 



a DIS F 2 (! + - R L/t)( 1 + £> R L/t) 

where 

4(1 -y)-^ 



(46) 



A(l-y) + 2y 2 + 

is the photon polarization parameter, and 



Z = ' V^" (47) 



R A(N) _ V Q 2 J L ^ Q 2 u t ^ a± 

L/T j*A(N) ~ 7*A(AT)' ^ C 

O rp (Jrp 



According to (46) the lepton nucleon cross section ratio is only equal to the structure function ratio, 
if there are no nuclear effects in R^ T or £ = 1. This is assumed in all NMC data and supported by 
experimental observations. In the kinematical region relevant for NMC, data for R^ T show practically 
no A dependence [UJ. Nevertheless this problem needs further study and more precise data. 

Already the first DIS measurements at SLAC jl2] showed that the structure function F2(x,Q 2 ) is 
nearly constant as function of Q 2 at fixed x. An explanation of this phenomenon was given by Bjorken 
[13] and by Feynman [44J. If the transverse momenta of the partons are neglected, then the cross 
sections for transverse and longitudinal photons scattering off spin-1/2 partons, i. e. quarks, are given 
by 

A-n 2 a em Z 2 f ( x 
°t = 7^ v 5 1 



Q 2 (l-x) \ x q 
o L = 0, (49) 

where x q is the momentum fraction of the proton carried by the struck quark and Zf is the flavor 
charge in units of the elementary charge. The 5-function arises from momentum conservation and gives 
a physical meaning to the Bjorken variable. In the Breit frame [15], x is the momentum fraction of the 
proton carried by the struck quark. For massless quarks, the longitudinal cross section is zero due to 
helicity conservation. Introducing the density qf(x) of quarks of flavor /, inside the proton, one obtains 
a simple partonic interpretation of the structure functions, 

F 2 (x) = xJ2z 2 f (qf(x) + q f (x)), (50) 
/ 

F L = 0. (51) 

In the naive parton model, the structure functions depend only on x and not on Q 2 , the longitudinal 
structure function vanishes and one obtains the Callan-Gross relation |4*6"1. 



F 2 - 2xF l = 0. (52) 

These equations are the basic results of the parton model and they are approximately confirmed by 
experiment. 

One of the great achievements of QCD is the successful description of the deviations from the naive 
parton model seen in experiment. In particular at low x deviations from Bjorken scaling become quite 
pronounced. In the QCD improved parton model, perturbation theory is applied to calculation of 
the corrections to the parton model predictions. Increasing the photon virtuality one enhances the 
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resolution of the parton probe, and one sees more partons. New partons appearing at higher resolution 
are generated by the splitting processes G — > 2G, G —> qq, q — >■ qG. The evolution of the structure 
function with Q 2 is described by the DGLAP equations [HJ HH SSI EQ], which for the singlet parton 
densities read, 



d 




2^ 



( 



P 



ff 



P fG 



P 



GG 



J 



<lf(xuQ 2 ) 
G(x 1 ,Q 2 ) 



(53) 



where qf is a quark of a given flavor, and the splitting functions P(x) are calculated perturbatively [5l] . 

All the soft physics is contained in the parton distributions. These are essentially contaminated by 
nonperturbative interactions and have to be parametrized at some input scale Qq. Different parametriza- 



tions have been provided by several collaborations EHl El] performing global analyses of data in 
leading and next to leading orders. With the parton distributions as input, one can then calculate F 2 
at a higher value of Q 2 . The other essential property is that the parton distributions are universal, i.e. 
they do not depend on the process under consideration, but only on the hadron state. This is based on 
QCD factorization proven for some processes and up to higher twist effects [55J. 



3.1 Phenomenology of shadowing in DIS 



The space-time picture of an interaction within the parton model varies significantly depending on 
the reference frame. Only observables are Lorentz invariant, but not our theoretical ideas about the 
dynamics of the interaction. In particular, what looks like absorption of the virtual photon by a parton 
in the nucleon in the Breit frame, looks very different in the rest frame of the proton. Here the same 
process looks like fluctuation of a high-energy virtual photon into a colorless qq dipole, with a subsequent 
interaction of the dipole with the target via gluonic exchanges. 

Correspondingly, shadowing also looks quite different. In the infinite momentum frame of the nucleus 
both the nucleus and each bound nucleon are Lorentz contracted by the 7-factor. A Lorentz-boosted 
nucleus looks like a pancake, as well as the bound nucleons. So if the nucleons do not overlap in the 
nuclear rest frame, they are still separated in the Breit frame also. However, the Lorentz 7-factor of 
the small-x partons is reduced by the factor x, so they contract much less and start overlapping, like 
is illustrated in Fig. [lOj This leads to a fusion of partonic chains originated from different bound 




Figure 10: Overlap of parton clouds originated 
from different nucleons in a fast nucleus due to 
the reduction of the Lorentz factor at small x 



nucleons [HB], and to a reduction of parton densities, and this produces shadowing. 

In the nuclear rest frame the physics of shadowing is more intuitive and corresponds better to its 
optical analogy. In this reference frame the gg-dipole fluctuation of the photon propagates through 
the nucleus and experiences multiple interactions, which lead to a reduction of the photon flux, and 
eventually to a suppression of the cross section. This is the same phenomenon as parton fusion, but 
seen from different reference frames. 
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The results of DIS measurements on nuclear targets are usually presented in the form of a ratio, 
Eq. for the structure functions 

The shadowing suppression factor R^/ N (x, Q 2 ) < 1 can be estimated using the same two-step correction 



as for real photons, as is illustrated in Fig. |4j It can be calculated with the same formula (16) replacing 
a sd ^ a ld, N > although the value of the absorption cross section remains problematic. 

Shadowing in DIS on nuclei begins at sufficiently high energies W 2 3> Q 2 , i.e. at small x <C 1. In 
this case the coherence time, also called Ioffe time, in the target rest frame reads, 



1 2u 

Ql Q 2 + M\ 2m N x' 



tc = ~= 7^TT72 « TiZT^- (55) 



One should expect the onset of shadowing at t c > 2fm, i.e. at x < 0.05, and saturation at t c 3> R^. 
Actually the saturation is not exact due to the contribution of large masses, increasing with 1/x. 



Expression (24) written in the "frozen" approximation, i.e. with no mixing of different eigen components 
during propagation through the nucleus, includes all the inelastic shadowing corrections and can be 
applied to DIS. Expanding the exponential one gets in the lowest order in multiple interactions 



Ra/n(x, Q 2 ) = l- \<Je ff (T A ) + ... , (56) 



where (T4) = (1/^4) / d 2 bT\(b) is the mean value of the nuclear thickness function, and 



(i a tot) 2 ) 1 r-, 

°zff = j-m ■ (57) 

\ a tot I 



It is interesting to notice that in this ratio the numerator and denominator are controlled by different 
scales. Let us classify the hadronic fluctuations of a highly virtual photon as either hard or soft, as is 
presented in Table [2] Naturally, hard eigenstate fluctuations interact weekly with the cross a\ ot ~ 1/Q 
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Table 2: Contributions of soft and hard fluctuations of a virtual photon to the DIS cross 
section and to nuclear shadowing. 



Fluctuation 


weight: Wi 




Wi a\ ot 


Wi (al ot ) 2 


Hard 


~ 1 




-1/Q 2 


~i/Q 4 


Soft 


~ v 2 IQ 2 


~i//i 2 


-1/Q 2 


~ i/^ 2 Q 2 



while soft ones have a large cross section a\ ot ~ l//i 2 , where /x is a soft hadronic scale. At the same time, 
hard fluctuations dominate in a highly virtual photon, while soft ones appear rarely with a probability 
suppressed at least as 1/Q 2 (1/Q 4 for longitudinally polarized photons). So smallness of the hard 
cross section is compensated by the large weight, and vice versa. Thus, both contributions, hard and 
soft, behave like 1/Q 2 , and as Table [2] shows, the DIS cross section (a\ ot ) gets a finite contribution 
from the soft interactions at any high Q 2 . However the shadowing term in (56), as well as diffraction, 
turn out to be dominated by soft interactions, whose contribution scales as 1/Q 2 similar to the DIS 
inclusive cross section. Therefore the x-dependence of a e ff(x,Q 2 ) should be pretty mild, (1/x) 5 , where 
5 = 2Ajp(/i 2 ) — Ajp(Q 2 ). The Pomeron intercept a_p(0) = 1 + Ajp(Q 2 ) rises with Q 2 , so the two terms 
in 5 essentially compensate. 
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Figure 11: Logarithmic Q 2 - derivative of the ratio 
of structure functions for tin to carbon. The data 
are from j5T]. 
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Figure 12: Nuclear shadowing versus scaling vari- 
able n(x, Q 2 ,A) (see text). The data for Li, C 
and Ca are from [SU [BO] ■ 



Another interesting consequence of the presence of two scales in the nuclear shadowing term in (56)- 
(57) is the absence of Q 2 dependence in the nucleus-related denominator, (cr l tot ). The Q 2 dependence 



of Ra/n(x,Q 2 ) comes from the nucleon structure function F^(x,Q 2 ) in the denominator. Comparison 



with NMC data presented in Fig. 11 demonstrated good agreement. This fact makes questionable the 
possibility of extracting the nuclear gluon density from Q 2 dependence of Ra/n(%,Q 2 ) via a DGLAP 
based analysis. Indeed, the DGLAP equations relate the logarithmic Q 2 derivative of the structure 
function at high Q 2 to the gluon density [61]. Correspondingly, one can extract gluon shadowing from 
the Q 2 variation of Ra/n(%, Q 2 ) as was done in [62J. However, data presented in Fig. 
with the lowest qq Fock component of the virtual photon, containing no gluons. 
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can be explained 



Furthermore, one can calculate the effective cross section Eq. (57) since the numerator ((cr^J ) 



according to Eqs. (21) and (27) is related to diffraction 



16tt 



eff 



a 



7 *iV 
tot I 



dM\ dp 2 T 



(5* 



pt=o 



So far we neglected in (56) only the higher order multiple interactions, hoping that cr e jj(T y i) is 



small. If, however, x is not sufficiently small for the "frozen" approximation, the phase shifts diminish 



the shadowing effect. These corrections should be done similar to Eq. 16, and another simplifying 
approximation is the one made in (55), where we fixed M\ = Q 2 . In this case Eq. (56) is generalized 



for the regime of shadowing onset as [57, 58J, 

^ A (x,Q 2 ) ^ 



R 



A/N 



a tot 

7* N 
Vtot 1 



x,Q 2 



-a eff (T A ) F 2 (q L ) 



n(x,Q 2 ,A), 



(59) 



where the longitudinal nuclear form factor 



A(T A ) 



d 2 b 



dz pA(b, z) e 



(60) 
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takes into account the effects of the finite coherence time t c w Eq. (55). At large qp > ^-/Ra, the 
nuclear form factor (60) vanishes and suppresses the shadowing term (59). This is easily interpreted: 



for large qp the fluctuation lifetime and its path in nuclear medium are short, and shadowing is reduced. 
In expression (59) all ingredients are calculable for different nuclei, and as functions of x and Q 2 . 



This invites the introduction of a scaling variable n(x, Q 2 , A) defined in (59) [53 EH], which should make 
shadowing universal for all nuclei, and values of x and Q 2 
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NMC data points [59 



60] for 

Data 



In Fig 

Ra/n(x, Q 2 ), each having specific values of x and Q 2 are plotted against the variable n(x, Q 2 , A) 
confirm the predicted scaling with good accuracy. 

In terms of parton distribution functions (PDF) one can disentangle shadowing for quark and gluons, 
based on the triple- Regge phenomenology. The cross section of the single-diffractive process, 7* + iV — » 
X + N, can be expressed in terms of the triple-Regge graphs. Indeed, summing up all final state 
excitations X, one can apply the unitarity relation to the 7*-iP amplitude, or to the total cross section 



o"7 ot , as is shown in Fig. 13 The latter, according to (pit) , is proportional to the structure function of 



E 

X 





N N N N 



Figure 13: The cross section of single diffraction, 
7*iV — > XN summed over all excitation chan- 
nels at fixed effective mass Mx- The 7*-iP total 
cross section at high collision energies is a sum 
of Reggeon and Pomeron exchanges. 





Figure 14: The Gribov correction with qq (left) 
and qq + g (right) intermediate states with large 
invariant mass, which is related to small frac- 
tional momenta either of a quark, or the gluon 
respectively. 



the target, i.e. the Pomeron, so one can say that this way one can measure the PDF of the Pomeron 
[63]. 

Provided that the effective mass of the excitation is large (but not too much), so <C M\ <C s, one 
can describe the Pomeron-hadron elastic amplitude via exchange of the Pomeron or secondary Reggeons 
in the t-channel. Then one arrives at the triple-Regge graphs, which lead to the cross section [64J, 



da 



sd 

dxF dt 



Gjpjp^I-xfY^- 2 ^ 



So 



Qi(0)-1 



(61) 



Here xp is the Feynman variable for the recoil nucleon in the cm. of collision, xp ~ 1 — M x /s; Gjpp]p(t) 
and G]p]PM(t) are phenomeno logical triple-Regge vertices. 

In terms of quark-gluon intermediate states the Reggeon and Pomeron exchanges correspond to the 
photon fluctuating either to heavy qq or qqg fluctuations respectively, as is depicted in Fig. [14| In the 



latter case what is important is that the large invariant mass M x is due to a large difference between 
the light-cone momenta of the qq and the g. In terms of PDFs, the Reggeon and Pomeron parts of aj f 
correspond to measurements of the quark, or gluon distribution functions. 

In the case of Gribov shadowing corrections the two Pomerons couple to different bound nucleons. 
Notice that the Gribov correction in the elastic ^*-A amplitude, Fig. |4j comes from different unitarity 
cuts [65J corresponding to single and double inelastic interactions in the nucleus, or single diffraction, 



as shown in Fig. 13 
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3.1.1 Models for Gribov corrections 



The Gribov correction to the nuclear PDF for the parton j was calculated within the usual approach 
[5], but formula (16) was presented in the form [661 E7J 168] . 

oo oo 

5f j/A (x } Q 2 ) = -87rRe J d 2 b J dz x p A (b, z x ) J dz 2 p A (b } z 2 )e~ 1 ^f T f T ^ b ' zl ' z ^ 

— oo z\ 

x f X "'° dxjp e^M^-^fDW (/?) g2> ^ Pt)pt= ^ (g2) 

J X 

where xjp = (M x + Q 2 ) / s; Mx is the invariant mass of the diffractive excitation; /3 = xjxp\ fj/$ is the 
diffraction PDF differential in four variables. The small real part of the amplitude and terms ~ 1/A 
are neglected. 



This approach suffers from the same problems as Eq. (16), namely, the sum over the higher order 



Gribov corrections, depicted in Fig. [TJ is replaced by the eikonal attenuation exponential with the 



effective cross section [66J written in the same way as Eq. (|58|) 

CXpfi 

Jj/N' /V "' 



= 7 ~T~7m / dX]p ff/N (Z 3 ' <? 2 ' t *nin) ( 63 ) 



There are, however, several concerns related to this assumption: 

• Within the range of Bjorken x where data on DIS on nuclei are available, this cross section 
is significantly varying during propagation through the nucleus. Only at very small x < 1CT 3 
hadronic fluctuations of a photon are "frozen" by Lorentz time dilation. Still, this never is precise 
enough for the fluctuations containing gluons, which are a source of gluon shadowing. 

• Even if a fluctuation is "frozen" , it is not correct to use the averaged cross section in the exponent, 



the whole exponential should be averaged. The example given above in Eq. (29) shows that the 



nucleus becomes much more transparent compared to the exponential attenuation assumed in ( 62 ) 



The difference is exactly the Gribov corrections. Notice that a better result was achieved in [68] , 
based on the Schwimmer model [69]. The nuclear transparency was found to depend on nuclear 



thickness as in (29). However the way how the phase shift was introduced in the calculations 
had no justification, and the validity of the Schwimmer model at rather large x, where data are 
available, is doubtful. 

Gluon shadowing is self-quenching. Once it becomes strong, it affects and suppresses the value of 



cr e ff, which is defined in (63) on a free nucleon. Correspondingly, gluon shadowing gets reduced. 
Such a self-quenching process leads to a nonlinear equation derived in [TDJ [7T] in the frozen dipole 
approach. 



The calculations with Eq. (62) and data for the diffractive structure function led to a magnitude of 
gluon shadowing, which considerably exceeds the results of other approaches [721 EH ESI ESI E7J and 
that of the dipole approach described below. 

A much more simplified modeling of the Gribov corrections was performed in [78J. Although this 
model is outdated compared to the contemporary state of art in the field, we highlight here several 
features of this description of shadowing, because the model is frequently used by experimental groups 
analyzing their data. The model exhibits all the problems listed above, and creates even more troubles. 

• The lowest order Gribov correction is calculated in |78j in analogy to the vector dominance model 
[22]. Namely, the invariant mass distribution of the diffractively produced intermediate states is 
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replaced by a delta function, 5(M — Q ), reducing the Gribov corrections to a single-channel 



problem. On the contrary, the above discussed analyses perform mass integration, Eqs. (16), (62) 
in accordance with the measured mass-dependence of the diffractive cross section. 



The cross sections of diffractive processes hp — > Xp in (16), or 7*p — >• Xp in (62), are well 
measured and fitted to data. Instead, in [7SJ this cross section is fitted to data on shadowing, 
which have incomparably lower statistics than in diffraction. 

The falling Q 2 dependence of the effective absorption cross section is motivated by the dominant 
contribution of heavy intermediate states of mass M 2 = Q 2 , which have smaller cross section. 
That is not correct, since the heavy hadrons have larger radius and larger cross section. This is 
the basis of the so called Bjorken puzzle [79], which is solved by a nontrivial cancellation between 
diagonal and off-diagonal diffractive amplitudes [80J. 



Shadowing at large Q 2 is known to be dominated by the aligned-jet hadronic components of the 
virtual photon [81] (see below Sect. 3.2.3), which lead to a logarithmic, rather than power 
Q 2 -dependence of shadowing 



We conclude that this model is oversimplified and has a very low predictive power compared with 
other approaches. 



3.2 Dipole representation 



The problems of Eq. (62) listed above are typical for a description based on the hadronic representation, 



since hadrons are the eigenstates of the mass matrix, but they do not have a certain size, controlling the 



nuclear transparency [13]. Nevertheless, this description has the advantage of having in Eq. (62) definite 



phase shifts, controlled by xjp related to Mx- On the other hand, switching to the dipole representation, 
one gets certainty with the dipole sizes, but the mass and phase shifts become uncertain. Later on, in 



Sect. |3.2.2| we will explain how to deal with the phase shifts within the dipole approach, but here we 
start assuming that the energy is high enough to neglect the phase shifts. 

The dipole description of [13] was applied to DIS in [SI] . The inclusive DIS cross section of a proton 
gets the form, 



t T,l( x iQ ) 



da / d 2 rT 



^ L (a,r T ,Q 2 ) 



cr, 



(r T ,x), 



(64) 



where the \l/^ L (a, r^, Q 2 ) are the light-cone (LC) wave functions for the transition 7* — > qq. The LC 
wave functions can be calculated in perturbation theory and read in first order in the fine structure 
constant a em [85 | l86 | 184], 

I / '.- ; ■ .' i (' / ! 

qq 



C/(f T ,«,Q 2 



2tt 



where x an d X are the spinors of the quark and antiquark respectively, 
momentum of the photon carried by the quark; 



a(l 



a)Q 2 + m 2 



(65) 

Here a is the light-cone 

(66) 



Ko(erT) is the modified Bessel function. The operators O t,l for transversely and longitudinally polarized 
photons have the form, 



3 T = m q a ■ e + i(l - 2a) (a ■ n) (e • V T ) + (<? x e) ■ V T 
6 L = 2Qa(l-a)a-n , 



(67) 
(68) 
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where the two-dimensional operator Vt acts on the transverse coordinate f^; n — q/q is a unit vector 
parallel to the photon momentum; e is the polarization vector of the photon. 

For the universal total cross section of dipole-proton interaction we rely on the parametrization [87] 
fitted to data for the proton structure function F 2 (x, Q 2 ) measured at high energies (small x) at HERA, 



cr, 



(r T , x) 



0-q 



exp 



"tQo 



A(x/xq) x ) 



(69) 



where Qq = 1 GeV and the three fitted parameters are o"o = 23.03 mb, xq = 0.0003, and A = 0.288. This 
dipole cross section vanishes oc r\ at small distances, as implied by color transparency [13], and levels 
off exponentially at large separations, which reminds eikonalization. Although this parameterization 
might be unrealistic at large separations (see discussion in [EE]), and does not comply with the DGLAP 
evolution (see improvements in [EH]) we will use it, because of its simplicity, and because DIS and 
Drell-Yan (see Sect. [4]) data are all available in the kinematical range where (69) works rather well. So 
in what follows we employ this parametrization for numerical calculations (unless specified). 



3.2.1 Shadowing and the lifetime of qq fluctuations of the photon 

Nuclear shadowing is controlled by the interplay between two fundamental quantities. 

• The lifetime of photon fluctuations, or coherence time. Namely, shadowing is possible only if the 
coherence time exceeds the mean inter-nucleon spacing in nuclei, and shadowing saturates (for a 
given Fock component) if the coherence time substantially exceeds the nuclear radius. 

• Equally important for shadowing is the transverse separation of the qq. In order to be shadowed 
the (/(/-fluctuation of the photon has to interact with a large cross section. As a result of color 
transparency (T3J |90j [91], small size dipoles interact only weakly and are therefore less shadowed. 
The dominant contribution to shadowing comes from the large aligned jet configurations 
of the pair. 



The lifetime of a hadronic fluctuation given by Eq. (55) can be presented as 

P 



tr 



xttln 



Pt 



where P = (1 + M 2 /Q 2 )~ l , and t™ ax = l/m^x. The usual approximation is to assume that M 2 rj 
since Q 2 is the only large dimensional scale available. In this case P 
For a noninteracting qq the coefficient P has a simple form, 



1/2. 



(70) 

Q 2 



P(k T ,a) 



Q 2 a(l 



a: 



k 2 +e 2 



(71) 



where m„ and kx are the mass and transverse momentum of the quark respectively. To find the mean 



value of the fluctuation lifetime in vacuum one should average ( 71 ) over kx and a weighted with the 
wave function squared of the fluctuation, 



\ qq 


P(k T ,a) 


qq 


\ qq 


qq J 



(72) 



The normalization integral in the denominator of ( 72 ) diverges at tt — > for transversely polarized 
photons, therefore we arrive at the unexpected result {P T ) vac = 0. This means that in vacuum a 



24 



transverse photon fluctuates mainly to heavy qq pairs with large kr, which have a vanishingly short 
lifetime. However, they also have a vanishing small transverse size tt ~ an d interaction cross 

section. Therefore, such fluctuation cannot be resolved by the interaction and do not contribute to the 
DIS cross section. To get a sensible result one should properly define the averaging procedure. We are 
interested in the fluctuations which contribute to nuclear shadowing, i.e. they have to interact at least 
twice. Correspondingly, the averaging procedure has to be redefined as, 



shad 



(f s d(j* ->■ qq) 


P{k T ,a) 


fsdd* ->■ qq)) 


(fsd(j* -> qq) 


fsd{l* ->■ qq)} 



(73) 



where f s d(l* — > qq) is the amplitude of diffractive dissociation of the virtual photon on a nucleon 
7*iV ->■ qqN. 

Thus, one should include in the weight factor for the averaging the interaction cross section squared 
CTqq(rT, s), where s = 1m N v — Q 2 + m 2 N . Then, the mean value of the function P(a, kr) reads, 



(p t - l ) 



Jq 1 da J d 2 r\d 2 r 2 


^ T q q L {r 2 ,a) 


°qq ( r 2, s) P (f 2 - fi, a) ^ L (fx, a) erg (n, s) 


Jo da I d2r 


<i L (r>X-(r, S ) 


2 



(74) 



with 



P(r 2 - fi,a) 



(2vr 



2exp ^-zA>r- (r 2 -ri)JP(a, fc T ) = K (e \r 2 -ri|) 



(75) 



which is a Fourier transform of Eq. (71). 



The results of numerical calculations for (P) , for transversely and longitudinally polarized photons, 
are depicted in Fig. 15 as function of Q 2 , and in Fig. 16 versus Bjorken x. One can see that (P L ) > (P T ), 



which means that a longitudinal photon produces lighter qq fluctuations than a transverse one. This 
is related to the suppression of the very asymmetric qq pairs with a — > or a — > 1 in the distribution 
amplitude of a longitudinal photon. Such asymmetric pairs have the largest invariant mass. 



3.2.2 The path integral technique 

As far as the lifetime of partonic fluctuations of a photon significantly exceeds the nuclear size, the dipole 
approach is very suitable and easy tool in order to calculate the shadowing effects in DIS. Indeed, in 



this case one can rely on the "frozen" approximation Eq. (24), which for interaction of a virtual photon 
has the form, 



a^{x,Q 2 ) =2 I d 2 bj da / d 2 r T ^{a,r T ,Q 2 

o 



1 _ e -l^ T ,x)T A {b) 



(76) 



where the ^/^ L (a,r^) are given by Eq. (65), and in the quadratic form read 



K(a,rr)| a 



2N c a e 



lN r a 



J2 Zf {[1 - 2«(1 - a)} e 2 K\{er T ) + m 2 f K 2 (er T )} , 



= ^^V(l-«) 2 Kg(er r ) 



(77) 
(78) 



25 




0.1 1 10 100 





1.0 r 




0.9 - 




0.8 - 




0.7 - 




0.6 - 








0.5 - 




0.4 - 




0.3 - 




0.2 - 




0.1 - 




0.0 L 





<P L > 












<P T > 



















0.0001 



0.001 



0.01 



0.1 



Figure 15: Q 2 dependence of the factor (P) de- 



fined in (70), at x = 0.01, for the qq fluctuations 
of transverse and longitudinal photons, and for 
qqG fluctuation, from top to bottom. Dotted and 
dashed and solid curves correspond to different 
forms of the dipole cross section and distribution 
functions (see details in [T7]). 



Figure 16: Bjorken x dependence of the factors 
(P T ' L ) and (P G ) defined in (70) corresponding 



to the coherence length for shadowing of trans- 
verse and longitudinal photons and gluon shad- 
owing, respectively. Solid and dashed curves cor- 
respond to Q 2 = 4 and 40 GeV 2 . 



The advantage of the dipole description is clear, Eq. ( 76 ) includes Gribov inelastic shadowing cor- 
rections to all orders multiple interactions [13 



what is hardly possible in the hadronic representation. 
On the other hand, the dipoles having a definite size, do not have any definite mass, therefore the 



phase shifts between amplitudes on different nucleons cannot be calculated as simple as in Eq. (16). A 
solution for this problem was proposed in 



If t c <Ra the "frozen approximation is not appropriate 
and one should correct for the dipole size fluctuations during propagation through the nucleus. This 
can be done with the path integral technique [94] . which sums up different propagation paths of the 
partons. 



For a qq component of the photon Eq. (76) should be replaced by 



x 



1 oo oo 

<7tot A ^j = A (crtot*^) ~~ J ^ J ^ a J J ^ Z2 J < ^' ri J < ^' r2 

— oo zi 

^fq 1 (e, r 2 ) Pa (b, z 2 ) <r£J (s, f 2 ) G (f 2 , z 2 \ F u z x ) p A (b, z x ) (s, n) (e, fx) 



(79) 



The Green's function G (f^, z 2 \ fi, z\) describes propagation of a qq pair in an absorptive medium, 
having initial separation f\ at the initial position zi, up to the point z 2 , where it gets separation f 2 , as 
is illustrated in Fig. [17| It satisfies the evolution equation, 



'. d A 

dz 2 



(r 2 ) 



2ua(l 
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U(r 2 ,z 2 ) 



G (f 2 , z 2 1 fi, zi) = iS (z 2 - zi) 5 {2) (f 2 - fi) 



JO) 
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Figure 17: Propagation of a gg-pair through a 
nucleus between points with longitudinal coordi- 
nates Z\ and z 2 . The evolution of the qq separa- 
tion from the initial, r*i, up to the final, f 2 , due to 
the transverse motion of the quarks, is described 
by the Green's function G (p 2 , z 2 ', Pi, Z\). 



Figure 18: Propagation through a nucleus of the 
qq — g fluctuation of a longitudinally polarized 
photon. Neglecting the small, ~ 1/Q 2 size of 
the color-octet qq pair, the effective octet-octet 
dipole propagation is described by the Green's 
function G gg (r 2g , z 2 \ r lg , Zx). 



The light-cone potential term in the left-hand side (l.h.s.) of this equation describes nonperturbative 
interactions within the dipole, and its absorption in the medium. The real part the potential responsible 



for nonperturbative quark interactions was modeled and fitted to data of F$ in 



Here we fix 



Re U(r 2 , z 2 ) = 0, and treat quarks as free particles for the sake of simplicity. The imaginary part of the 
potential describes the attenuation of the dipole in the medium, 



Im U(r, z) 



— - cr, 



'r)p A {b,z) 



At small x when the coherence length substantially exceeds, Zj ,L ^> {z 2 — z\) (the nuclear radius) the 
solution of Eq. (80) much simplifies, G(f 2 , z 2 | r*i, z\) oc 5^ (f 2 — r*i), so Lorentz time dilation "freezes" 
the variation of transverse qq separation. Correspondingly, the total cross section gets the simple form 
of Eq. (f76|) 

Eq. (80) can be solved analytically if the medium density is constant, p^ 
cross section has the simple form, cr q q(r) = Cr 2 
with a complex frequency 



p , and then the dipole 
The solution is the harmonic oscillator Green function 



G (f 2 , z 2 ; fx, zi. 
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where 



2tc sinh (coAz 



Az 
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[rl + rf) coth (wAz) 



2f 2 ■ n 
sinh (uAz) 



z 2 — Z\, 

iCp 
va (1 — a) ' 
— % Cpova (1 — a) 



(82) 

(83) 
(84) 
(85) 



This formal solution and Eq. (80) properly account for all multiple scatterings and for the finite lifetime 



of the hadronic fluctuations of the photon, as well as varying transverse separation of the qq pair during 
propagation through the medium. It effectively sums up all the Gribov corrections including nonzero 
phase shifts between different amplitudes. 

For practical applications these results can be extended to a realistic nuclear density varying with 
coordinates as is described in [15]. Also the real part of the light-cone potential U(r,z) in (80) was 
modeled in [881 Ej to incorporate nonperturbative effects. 

The best data available today for nuclear shadowing in DIS are for the structure function ratio tin 



over carbon from NMC [95]. These data are shown in Figs. 19 and 20 as function of Bjorken x and 
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Figure 19: The Bjorken-x dependence of nuclear 
shadowing in DIS for the structure function ra- 
tio of tin relative to carbon. The data are from 
the NMC experiment [95]. The solid and dashed 
curves are calculated including or excluding the 
nonperturbative effects respectively. 
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Figure 20: The Q 2 dependence of nuclear shad- 
owing in DIS for the ratio of the structure func- 
tions of tin to carbon. The data 195] are the same 



as shown in fig. [19j and the curves have the same 
meaning. 



Q 2 respectively. The numerical results of the calculations, which are performed either disregarding or 
including the nonperturbative effects, are plotted in Fig. 19 by dashed and solid curves respectively. 
One can see that inclusion of the nonperturbative effects does not lead to a significant change of the 
magnitude of shadowing. Comparison with NMC data E5] shows pretty good agreement. Note 
however that inclusion of the antishadowing effect might shift the curves upwards. Indeed, the low-x 
tail of the EMC suppression observed at at large x should be an enhancement, as is suggested by the 
momentum conservation sum rule. Whether this will happen at all values of x or only around x ~ 0.1, 
depends on the model. 

As an example of the A dependence of shadowing, the numerical results for the ratios of various 
nuclei to carbon are depicted in Fig. 21 in comparison with NMC data [951 196]. Calculations including 
no interaction within the qq pair are shown by dashed curves. The calculations are parameter free and 
agree with data rather well. Introduction of a nonperturbative interaction between q and q described 
by the real part of the light-cone potential Re<7(r, z), parametrized in the oscillator form (see details 
in [SSI EH]), does not produce any significant change. These results, plotted in Fig. 21 by solid curves, 
also agree with the data. 

We also compare our calculations with data from the E665 experiment [H7J EH] , which covered much 
lower values of x, Fig. 22 The agreement with E665 data is not as good as with the NMC data. Indeed, 
the two datasets seem to be somewhat inconsistent [95] HI], so it looks challenging to reproduce both. 

The disagreement between NMC and E665 vanishes however, in the F 2 ratio relative to carbon. 
This observation might give a clue to the origin of the "disagreement" between the results of the two 
experiment. The modification of the properties of bound nucleons compared with free ones is a popular 
interpretation of the EMC effect, which is the observed suppression of nuclear F 2 (x, Q 2 ) at large x. This 
effect should propagate down to small x as an enhancement (momentum conservation sum rule), which 
is indeed observed at x ~ 0.1 — 0.2. Such a few percent enhancement should have similar magnitudes 
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Figure 21: Comparison between calculations for 
shadowing in DIS and experimental data from 
NMC [951 EB] for the structure functions of dif- 
ferent nuclei relative to carbon. Q 2 ranges within 
3 < Q 2 < 17GeV 2 . The curves have the same 
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Figure 22: E665 data for shadowing in DIS for 
various nuclei [971 EE]- Full circles show data 
taken with electromagnetic calorimeter cuts, 
while the open circles result from the FERRAD 
radiative correction code. Hadronic cuts were ap- 
plied for the triangles in the lower left plot. Q 2 
ranges within 0.15 < Q 2 < 22.5 GeV 2 , except 
for the xenon data, 0.03 < Q 2 < 17.9 GeV 2 . The 



meaning as in Fig. 19 



curves have the same meaning as in Fig. 19 



for carbon and heavier nuclei, because the nuclear density does not vary much. However for deuterons 
the EMC effect and related small-x enhancement are much weaker. Therefore the ratios A/D should 
be shifted by few percent upwards compared with A/C ratios. This is what happens and can be seen 



in Figs. [21] and [22} Our calculations of shadowing have not been corrected for the EMC effect of 
antishadowing, and this is probably why they agree with the NMC data for A/C but underestimate 
the E665 data for A/D. 

Notice that the measurement of nuclear shadowing in inclusive DIS, where only the final lepton is 
detected, is a very difficult task because of large radiative corrections [99J. These corrections occur, 
e.g. when the incident lepton reduces its momentum due to electromagnetic bremsstrahlung and does 
not produce a deep inelastic event jlOOj . In the older publications [591 HOI] , the NMC collaboration 
calculated the radiative corrections employing the computer program FERRAD, which relies on the 
theoretical analysis |102j . For a reevaluation of the shadowing data [HE] three different codes were 
tested, FERRAD, an improved version of FERRAD and TERAD. The last code relies on calculations 
|103[ I104[ 1105] and was used for the published NMC result. A more detailed discussion on radiative 
corrections for NMC and HERMES can be found in |106j . One may say that the whole shadowing 
effect in data is calculated. Without radiative corrections, the cross section ratio would be around unity 
[106J. However, the correctness of the calculation was checked by comparing to the "hadron tagged" 
data, making sure that really a deep inelastic event was measured. 

A different way to identify the deep inelastic events was chosen by the E665 collaboration [971 EE] , 
where cuts in the electromagnetic calorimeter were applied. In [97] the xenon data obtained in this way 



were compared to an evaluation with hadron tagging, see Fig. 22. The two methods give consistent 
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results. In |98j, the E665 collaboration also applied the FERRAD code for a comparison with NMC 
The results are depicted by open circles in Fig 
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One recognizes a systematic discrepancy between 
the two evaluation methods. The radiative corrected ratios are lower by ~5% than the points from the 
calorimeter analysis. 



3.2.3 Shadowing of longitudinal vs transverse photons 



Eqs. (64) and (80) allow to calculate the absolute values of ai and cry, as well as the magnitude of 
nuclear shadowing for each of them. The results for R p L i T = cr p L /a T calculated at x — 0.01 are plotted 
as function of Q 2 in Fig. 23 The solid and dashed curves, like previously, correspond to calculations 
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Figure 23: Ratio R p L j T on a proton as function of 




Q 2 calculated with Eqs. (76)-(78). For the curve 



notations see Fig. 19 The NMC data |107] are 
from the left to the right for x = 0.0080, 0.0045, 
0.0125, 0.0175, 0.025, 0.035, 0.05, 0.07, and 0.09. 



Figure 24: Ratio of longitudinal to transverse 
cross sections Rl/t = o"l/o"t calculated for nitro- 
gen 14 N divided by the same ratio for a proton, 



~R L f T l 'R P L / T - For the curve notations see Fig. 
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with or without inclusion of the nonperturbative effects in the qq distribution amplitude of the photon. 
Theory agrees well with NMC data |107j measured at close values of Bjorken x. 

The ratio turns out to be quite small, R L / T < 0.2. This is understandable intuitively as a result 
of a smaller transverse size of qq fluctuations in longitudinally polarized photons. Indeed, the mean 
transverse separation squared, 



(4) 



Q 2 a(l — a) + m 



2 - 



36) 



is minimal for the symmetric fluctuations, a ~ 1/2, but reach the maximal size of the order of the 
confinement radius for very asymmetric pairs at a — > 0, 1. The latter configurations are suppressed by 



the distribution function of longitudinal photons, Eq. (78), so the qq fluctuations are more symmetric, 
i.e. have smaller size, while fluctuations of transverse photons get a considerable distribution from 
asymmetric, large size fluctuations. 



This observation has important consequences for the relation Eq. (46) connecting the shadowing 
effects in the nuclear structure function and nuclear DIS cross section (J^dis- The observed nuclear 
ratio ([T| for the DIS cross section equals to that for the structure functions F 2 only if (i) either Rl/t ^ 1 
and can be neglected in (46); (ii) or £ — 1 ■ 1 D "" " ,,/ ' ,! ..-!,„,,,,> d \ 

<p 



1; (iii) or Rl/t has no nuclear dependence, R^j T 



-dN 
1X1 LjT- 



Regarding the condition (i), indeed R P L / T is rather small, but not that small to neglect it in (46). 
Regarding (ii), the value of £, according to (47), at small x is mainly controlled by the transferred 
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fractional energy y = vjE\. Apparently, at the same of x and Q 2 the value of v is also the same in 
high-energy (NMC, E665) and low-energy (JLAB, HERMES) experiments. Therefore the values of y 
and £ are very different. The low-energy experiments at small x ~ 0.1 are mainly done with large 
y — > 1, i.e. small £ <C 1, while at high energies ?/ C 1 and £ ~ 1. In the latter case, according to (46) 
one can safely treat the ratio of the DIS cross sections as that for F<i- 

However, at smallest x even in large energy experiments the transferred fractional energy reaches 
large values, so £ is small, and one should have a look at the A-dependence of R^/ T . This ratio was 

calculated in [TTJ for nitrogen 14 iV and hydrogen at x = 0.01, and the ratio of the results J^l/t/^l/t 



is plotted in Fig. 24 as function of Q 2 . We see that the difference between Rl/t on 14 N and p ranges 
within about 10/CThis correction is diluted by the small value of Rl/t ~ 0.2 down to about 2% even 
in the case of £ — > where the nuclear corrections are maximal. 

Notice that the weak nuclear dependence of Rl/t predicted in [T7], was at the time of publication 
in dramatic disagreement with HERMES data |108j . which showed a nuclear enhancement of Rl/t by 
a factor of several. However, few years later the HERMES Collaboration discovered a mistake in their 
measurements [109J. 



3.2.4 Nuclear shadowing for valence quarks 

Nuclear shadowing for valence quarks is usually believed to be small |110j . if it occurs at all. It was 
demonstrated in [1 llj . however, that shadowing for valence quarks is quite sizable, even might be 
stronger than the shadowing of sea quarks. Notice in this regard that the nuclear structure function 
F 2 (x) is different from the quark distribution function in an essential way; namely, the former contains 
shadowing effects and therefore the baryon number sum rule is not applicable to it |112j . a difference 
that might explain the discrepancy of present results compared to Ref. [110] . 



Note that we relate the nuclear cross section Eq. (76) to shadowing for sea quarks because the 



dipole cross section a^(rx, x) was fitted to HERA data at very low x < 0.01, so it includes only 
the part corresponding to gluonic exchanges in the cross- channel. Therefore, this is the part of the 
sea generated via gluons (there are also other sources of the sea, for instance the meson cloud of the 
nucleon, but they steeply vanish with x). The fact that the color-dipole cross section includes only the 
part generated by gluons is the reason why it should not be used at larger x. This part of the dipole 
cross section can be called the Pomeron in terms of Regge phenomenology. In the same framework, one 
can relate the valence quark distribution in the proton to the Reggeon part of the dipole cross section, 
which has been neglected so far. So, to include valence quarks in the dipole formulation of DIS, one 
should replace 

Vqq(r T , x) => (Tq q (r T , x) + a^(r T , x) , (87) 

where the first (Pomeron) term corresponds to the gluonic part of the cross section, responsible for 
the sea quarks in the nucleon structure function. The second (Reggeon) term must reproduce the 
distribution of valence quarks in the nucleon; this condition constraints its behavior at small x. One 
can guess that it has the following form, 



(7^(r T ,x) = Nr^Vx, (88) 

where \fx should reproduce the known x dependence of valence quark distribution (as, in fact, motivated 
by Regge phenomenology), and the factor r\ is needed to respect the Bjorken scaling. The factor N 
will cancel in what follows. 

We are now in a position to calculate shadowing for valence quarks by inserting the cross section 



Eq. (|87|) into the eikonal expression Eq. (|76|). If one expands the numerator in powers of cr^ and picks 
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out the linear terrrj^J then one arrives at the following expression for nuclear shadowing of the valence 
quarks, 



J d 2 bT A (b) da J d 2 r T V qq (r T , a, Q 2 ) a*(r T , x) e -l*T,*)T A (&) 



A da J d 2 r T ^ qq {r T ,a,Q 2 ) 



S9) 



The results of numerical calculations with this expression are compared with the results of the global 
fit |llUj in Fig. 25 We see that shadowing for valence quarks is stronger than for sea quarks, and is 
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Figure 25: Shadowing for sea and 
valence w-quarks in DIS off gold at 
Q 2 = 20 GeV 2 . Solid lines are cal- 
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much stronger than in the parameterization of [110] . Note that the global fit of [110] assumes that 
the nuclear valence quark distribution has to satisfy the baryon number sum rule, which takes it for 
granted that the structure function is a measure of the number of quarks in the target. However, the 
very meaning of shadowing is a reduced access of the probe (the virtual photon) to some of the bound 
nucleons. So the reduction of the effective number of valence quarks in the target, does not mean that 
baryon number is not conserved, it only shows that the probe is not sufficiently hard. Increasing Q 2 
one does not get rid of the soft component of DIS [58] . 

Unfortunately, it will be impossible to extract the low-x valence quark distribution of a nucleus from 
DY experiments, because the nuclear structure function is dominated by sea quarks. Maybe neutrino- 
nucleus scattering experiments and accurate measurement of shadowing for Fg N + F£ N (see Sect. [5J 
can provide information on shadowing of valence quarks. 



3.2.5 "Antishadowing" 

Shadowing is a quantum-mechanical phenomenon, which is impossible for hard reactions in classical 
physics, since multiple interactions have too small cross sections to shadow each other. Only if a hadronic 
quantum fluctuation has a sufficiently long lifetime to undergo multiple interactions (at least twice), can 
two amplitudes on two different nucleons interfere, either destructively (shadowing), or constructively 
(antishadowing). This is why the term antishadowing is labelled by quotation marks in the title of this 
section, because this small enhancement of the nuclear ratio F^(x, Q 2 )/F^(x, Q 2 ), which was observed 
[113J at x ~ 0.1 — 0.2, is well outside of the coherence region. Therefore no coherence phenomena, 
either shadowing, or antishadowing are possible. Indeed, the fluctuation lifetime given by the Ioffe 



time, Eq. (55), for this x-interval is very short 0.5 — lfm compared to the mean nucleon spacing in 



3 The small size of a^(rT,x) at small x motivates such an expansion; however, one should note that it would not be 
proper to include the higher powers of the Reggeon cross section. Indeed, the Reggeons correspond to planar graphs. 
These cannot be eikonalized since they lead to the so-called AFS (Amati-Fubini-Stangelini) planar graphs, which vanish 
at high energies @]. 
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nuclei, ~ 2fm. Although the observed nuclear enhancement cannot be related to antishadowing, this 
terminology is widely adopted so we will use it. 

No consensus has been reached so far regarding the mechanisms responsible for this effect, as well 
as for the EMC effect, which is a nuclear suppression at larger x ~ 0.2 — 0.7 [113j . As far as both of 
these collective effects cannot be caused by coherence, they may be only related to a modification of 
the bound nucleons compared with free ones. This subject goes far beyond the scope of the present 
paper focused on coherence phenomena, but a detailed discussion of the medium effects can be found 
in the comprehensive review [114] . However, those effects might propagate down to smaller x and also 
contribute to the observed nuclear suppression, therefore we briefly discuss here possible implications 
at small x. 

Suppression of the quark PDF of a medium-modified bound nucleon at large x, should cause an 
enhancement at smaller x due to the baryon number (number of valence quarks) conservation sum 
rule. This enhancement would propagate down to small x, however shadowing onsets at x < 0.1 and 
overtakes the enhancement. In this qualitative picture the nuclear enhancement observed in the interval 
0-1 ~ x 0-2 is a result of the interplay of two phenomena, which have quite different origin. One, the 
EMC effect, has no relation to coherence, and results from medium modification of the properties of 
bound nucleons. Another, the shadowing suppression, is caused by destructive interference of the DIS 
amplitudes on different bound nucleons, regardless whether these nucleons are, or are not modified by 
the medium. A quantitative analysis based on the model of " swelling nucleons" was performed in [115] 
in good agreement with data. 



3.2.6 Coherence time for gluon radiation and gluon shadowing 

Shadowing in the nuclear gluon distributing function at small x, which looks like gluon fusion gg — >■ g 
in the infinite momentum frame of the nucleus, should be treated in the rest frame of the nucleus as 
shadowing for the Fock components of the photon containing gluons. Indeed, the first shadowing term 
contains double scattering of the projectile gluon via exchange of two i-channel gluons, which is the 
same Feynman graph as gluon fusion. Besides, both of them correspond to the triple-Pomeron term in 



diffraction which controls shadowing (see Eq. (61) and Fig. 13). 

The lowest Fock component containing a gluon is the \qqg). The coherence time controlling shad- 
owing depends on the effective mass of the \qqg), which should be expected to be heavier than that for 
a \qq). Correspondingly, the coherence time (if) should be shorter and the onset of gluon shadowing is 
expected to start at smaller x. 



For the coherence time one can rely on the same Eq. (55), but with the invariant mass of the 
fluctuation, 

^^ + T^> <») 

where a g is the fraction of the photon momentum carried by the gluon, and Mq q is the effective mass of 
the qq pair. This formula is, however, valid only in the perturbative limit. It is apparently affected by 
the nonperturbative interaction of gluons, which was found in [HE] to be much stronger than that for a 
qq. Since this interaction may substantially modify the effective mass M qqg , we switch to the formalism 



of Green's function described above, which recovers Eq. (90) in the limit of high Q . 

We treat gluons as massless and transverse. For the factor P defined in (70), in the case of gluon 
shadowing one can write, 

N 9 

where 



n = i>i • < 9i > 



N ° = mNX /^^^**-<*«.<''»K)»k,«..* 



qqi Qqi C^g) 
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dZ2 Gggg (f 2 Q, ?2qqi z 2 , ^Igi ^lqq, Z\) ^ qqg { 



(92) 



D 9 



d 2 r lg d 2 r lqq d 2 r 2g d 2 r 2qq da q d\n(a g ) V\ qg (f 2g , f 2qq , a q , a g ) 



$ {2) (r 2g - f lg ) 5 {2) (f 2qq - f lqq ) ^ qqg (fi g , f lqq , a q , a g ) 



(93) 



Here we have introduced the Jacobi variables, r, 



Rn 



R q and f g 



R g - (a q R q + a q R q )/(a q + a q 



R g ,q,q are the position vectors of the gluon, the quark and the antiquark in the transverse plane and 
otg,q,q are the longitudinal momentum fractions. 

Differently from the case of a \qq) Fock state, where we found that at high Q 2 perturbative QCD 
can be safely used for shadowing calculations, the nonperturbative effects remain important for the 
\qqg) component even for highly virtual photons. High Q 2 squeezes the qq pair down to a size ~ l/Q, 
while the mean quark-gluon separation at a g <C 1 depends on the strength of gluon interaction which 
is characterized in this limit by the parameter 6 ~ 0.65 GeV |88j . The presence of such a semi-hard 
scale, which considerably exceeds Aqcd, is confirmed by various experimental observations [116], in 
particular by the observed strong suppression of the diffractive gluon radiation |88j . In nonperturvative 
QCD models this scale is related to the instanton size [117[ 1118] , tq ~ 1/bo ~ 0.3 fm. 

For Q 2 ^> 6q the qq is small, r? <C r 2 and one can treat the qqg system as a color octet-octet dipole, 



as is illustrated in Fig. 18 Then the three-body Green's function G qqg factorizes 
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(94) 



The color octet-octet Green's function G gg , describing the propagation of a glue-glue dipole through 
the medium, satisfies the simplified evolution equation [HE], 



d_ 

dz 2 



2u 



+ 



A± (r 2g 



ft 4 r 2 

u ' 2g 



2ua g (1 - a g ) 



1va g (1 



ol. 



Ggg {f 2 g, Z 2 \ fg, Z X ) = 5 (z 2 ~ Z X ) 5 {2) (f 2g - 



Correspondingly, the modified qqg wave function simplifies too, 
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where the nonperturbative quark-gluon wave function has a form 

Vqg{r g ) = lim Jf gg {a g ,r g ) ■ 
and the color-octet dipole cross section reads, 

N ( \ 9 N / \ 
°gg \ X ^ T a) = l a qq\ X i r 9) 



°gg (X, Tg) 



bl 2 
exp ( - — r g 



(95) 
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(98) 



Within these approximations we can evaluate the factor (P ) = (l 9 c )/l™ ax given by Eq. (91). The 



results are depicted in Figs. 15 and [16] With the approximations made above, the calculations cannot 
cover the low Q 2 region and are perform at Q 2 > 1 GeV 2 . The gluon coherence length turns out to be 
much shorter than both /J and for \qq) fluctuations. This observation corresponds to delayed onset 
of gluon shadowing, shifted to smaller x compared with quark shadowing, as was predicted in [88J. 

Although gluon shadowing is related to the higher Fock component of the photon, 7* — > qqg, 
shadowing might be related also to a large qq separation, or g — qq. The former is related to quark 
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shadowing, which is a higher twist effect, but even increasing Q 2 one cannot get rid of this contribution. 
As was demonstrated above, the main contribution to shadowing of transversely polarized photons 
comes from very asymmetric sharing of light-cone momentum, a ~ m 2 /Q 2 . To get a net effect of 
gluon shadowing one should suppress this term by taking either heavy flavors |119[ I120j . or DIS with 
longitudinally polarized photons [88]. Here we rely on the latter. 

Longitudinal photons can serve to measure the gluon density because they effectively couple to 
color-octet-octet dipoles. This can be understood in the following way: the light-cone wave function 
for the transition 72 — > qq does not allow for large, aligned jet configurations [86]. Thus, unlike the 
transverse case, all qq dipoles from longitudinal photons have size 1/Q 2 and the double-scattering term 
vanishes like oc 1/Q 4 . The leading-twist contribution for the shadowing of longitudinal photons arises 
from the \qqg) Fock state of the photon. Here again, the distance between the q and the q is of order 
1/Q 2 , but the gluon can propagate relatively far from the gg-pair. In addition, after the emission of the 
gluon, the pair is in an octet state. Therefore, the entire |ggg)-system appears as a gg-dipole, and the 
shadowing correction to the longitudinal cross section is just the gluon shadowing we want to calculate. 

One can also see that from the expression for the cross section of a small size dipole [93j 1121] , 



a 



-A,N, 

qq 



7T 
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where <7at(x, Q 2 ) is the gluon density at Q 2 ~ 1/rf,. Thus, we expect nearly the same nuclear shadowing 
at large Q 2 for the longitudinal photoabsorption cross section and for the gluon distribution, 



v L a(x,Q 2 ) _ 9a(x,Q 2 ) 
aj^(x,Q 2 ) g N (x,Q 2 ) 



The shadowing correction to o\[x, Q 2 ) has the form (compare with (80)) 
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Assuming Q 2 3> 6q we can neglect r g . The net diffractive amplitude F 7 *^q qg (fi g , r lqq , a q , a g ) takes 

the form of Eq. (96), and we can rely on the factorized relation (94) for the 3-body Green's function, 



with equation (95) for the evolution of the gluonic dipole. The latter has the solution, 
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and other notations are the same as in Eq. (82), except C gg = \C according to (98). 



The results of numerical calculation of ( 



R 9 (X,Q 2 ) 



01 ) for the ratio 

g A (x,Q 2 ) 



Ag N (x,Q 2 Y 

are depicted in Fig. 26 as function of Bjorken x for Q 2 = 4 and 40 GeV 2 . 



(102) 



(103) 



(104) 



35 



0.7 



' C 



10 



10 



X 



0.6 



x=10"° 

6 8 10 12 14 
L(fm) 



Figure 26: Ratio of the gluon distribution func- 
tions in nuclei (carbon, copper and lead) and nu- 
cleons at small Bjorken x and Q 2 = 4 GeV 2 (solid 
curves) and 40 GeV 2 (dashed curves). 



Figure 27: Gluon shadowing vs. the length of the 
nuclear medium L = 2^/B? A — b 2 . All curves are 
for Q 2 = 20 GeV 2 but for different values of x. 



Our results for gluon shadowing as a function of the length of the nuclear medium at impact 



parameter b are shown in Fig. 27 The calculations are performed for lead with a uniform nuclear 
density of pa = 0.16 fm~ 3 . The small size of the gg dipole leads to a rather weak gluon shadowing. 
For most values of x, gluon shadowing increases as a function of L as one would expect. At the largest 
value of x = 0.01, however, gluon shadowing becomes smaller as L increases, and R g approaches 1. 
Although this behavior seems to be counterintuitive, it can be easily understood by noting that at 
x = 0.01 the coherence length of the \qqg)-Fock state becomes very small and the form factor of the 
nucleus suppresses shadowing [T7] . 



4 Drell-Yan process 

The Drell-Yan (DY) process, 



h 1 + h 2 ^ll + X, (105) 



in the kinematical region where the invariant dilepton mass M is small compared to the center of mass 
energy -y/s is of similar theoretical interest as DIS at low Bjorken x. Moreover, the cross sections of these 
processes are related by the factorization theorem In contrast to DIS, where only the total cross 
section can be measured, there is a variety of observables which can be measured in the DY process, 
such as the transverse momentum distribution or the angular distribution of the lepton pair. 
The fractional light cone momenta of the dilepton relative the colliding hadrons, 



x x = %] x 2 = (106) 
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satisfy the relations 
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X\ — x 2 = xp, (107) 
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where M is the dilepton invariant mass, and xp is the Feynman variable. 

Within the collinear approximation of the lowest order parton model pt = 0, and the fractional 



light-cone momenta Eq. (107) coincide with the standard Bjorken X\p of the q and q annihilating into 



the heavy dilepton, qq — > 11. In this lowest order approximation the cross section of the process (105) 
is expressed in terms of quark and antiquark densities, 

da 47ro£ 



em 



dr 3N r M 2 



[ dx i y2z 2 f {q f (x 1 )q f (r/x 1 ) + (1^2)}. (108) 
Jo f 

Here r = M 2 /s, and the factor N c (number of colors) appears in the denominator, because quark and 
antiquark must have the same color in order to annihilate. 

Some features of dilepton production cannot be understood in the lowest order approximation. 



Cross sections calculated straightforwardly from (108) are too small by a factor of 2 — 3 compared 



to the measured value. This discrepancy is usually treated by the introducing of a so called 
-fT-factor. The .ff -factor is approximately independent of M 2 , but it is process dependent. 

Dileptons with large transverse momenta, of the order of few GeV, are observed in experiment. 
There are however no transverse momenta in the naive parton model, which treats the process 



in the collinear approximation, Eq. (108). Phenomenologically, one can introduce a primordial 
momentum distribution of the quarks. However, the width necessary to describe data is much 
larger than what one would expect from Fermi motion of the quarks. 

These problems can be settled by taking into account the first order QCD corrections. Due to the 
radiation of the gluon, the quark acquires a transverse momentum. In this way, the pQCD corrections 
provides the missing mechanism for the production of lepton pairs with large transverse momentum pp 
[1221 CL23J [TM]. 

The next-to-leading (NLO) correction solves most of the problems of the naive parton model. It 
explains how large pr-dileptons are produced and account for almost all of the i^-factor |125j . Yet, 
not all the problems are settled by this correction. Since it is numerically large, one has to investigate 
how much higher order corrections affect the result. Furthermore, the transverse momentum spectrum 
is not well described even qualitatively. The theoretical result agrees with data only at p\ ~ M 2 , and 
even diverges at pt — > 0: da jdp\ oc p^ 4 , while the observed cross section is of course finite. The reason 
for this behavior is that large logarithms ln(M 2 /p^) occur in higher order corrections and one has to 
re-sum all these terms. This is possible within pQCD [126] I127[ I128j . by a re-summing soft gluons 
radiated from the quark or the antiquark. 

4.1 Dipole description of heavy dilepton production 

Although cross sections and all observables are Lorentz invariant, the partonic interpretation of the 
process depends on the reference frame. We have seen this in the case of DIS, and similar considerations 
can be made for DY. It was pointed out in |129j that in the target rest frame, DY dilepton production 
looks like bremsstrahlung, rather than parton annihilation. A quark (or an antiquark) from the projectile 



hadron radiates a virtual photon on impact on the target, as is illustrated in Fig. 28 The radiation can 
occur before (not shown) or after the quark scatters off the target. The sum of the two contributions 
gets the form in impact parameter representation similar to the dipole interaction in DIS. Such a dipole 
description of the DY process was proposed in |129] and developed in many publications including 

pnu una nnni nn em]. 

The cross section for radiation of a virtual photon from a quark after scattering on a proton, can be 
written in factorized light-cone form, 

da{qp^YX) f , 2 I ,t,t,l / 



d\n a 



J d 2 p\^ L q (a,p)\ 2 a q q(x 2 ,ap), (109) 
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Figure 28: Feynman diagrams squared for the 
cross section of di-lepton production in a q-N 
collision. In the target rest frame the quark (or 
antiquark) scatters off the target color field and 
radiates a massive photon, which decays into the 
lepton pair. The photon can be radiated before 
or after the quark hits the target. G(x) = xg(x) 
is the gluon distribution function in the proton. 
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Figure 29: Data from E772 [130J and solid curves 
are for the DY cross section in pD collisions at 
y/s = 38.8 GeV, dashed curves are for y/s = 
500 GeV. For each energy, the upper and lower 
curves correspond to massless quarks, or mf = 
200 MeV. 



similar to the case of DIS. Here, a q q is the dipole-proton cross section, which depends on the qq separation 
ap, where p is the photon-quark transverse separation and a is the fraction of the light-cone momentum 



of the initial quark taken away by the photon. In ( 109 ) T and L stand for radiation of transverse and 
longitudinal photons respectively. 



An interesting feature of (109) is the appearance of the dipole cross section, although there is no 



physical gg-dipole in Fig. 28 The physical interpretation of (109) is similar to DIS. The projectile quark 



is expanded in the interaction eigenstates. We keep only the first two, 



\q) = VZ: 



h\qbare) + ^f q \qi*) + • 



(110) 



where Zi is the wave function renormalization constant for fermions. In order to produce a new state the 
interaction must resolve between the two Fock states, i.e. they have to interact differently. Since only 
the bare quarks interact in both Fock components the difference arises from the relative displacement in 
transverse plane of the quark radiating the photon. If p is the transverse separation between the quark 
and the photon, the j*q fluctuation has a center of gravity in the transverse plane which coincides with 
the impact parameter of the parent quark. The transverse separation between the photon and the center 
of gravity is (1 — a)p and the distance between the quark and the center of gravity is correspondingly 
ap. A displacement in coordinate space corresponds to a relative phase factor in momentum space, 
— exp(iap • kj_), which produces the color screening factor exactly like in the dipole cross section. 
The LC distribution functions for the |?7*) component can be written in a form similar to DIS, 



where x are two component spinors, and 
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The operators 1 * |131] look similar to, but are different from the analogous operators Eqs. (67) and 



(68) in the case of DIS, 
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The two dimensional gradient V p acts on the transverse coordinate p\ n is the unit vector parallel to the 
momentum of the projectile quark; and a are the Pauli spin-matrices. The LC distribution functions 
for the transition q — > 7*g read explicitly, for a given flavor of unit charge 
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(114) 
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Comparing (115) and (116) with their DIS counterparts (77) and (78) shows that the factor N c is no 



longer present in the q — > 7*g LC wavefunctions. This is due to the N c in the denominator of (108). 
Furthermore, we see that \^ tq (a, p)\ 2 has an extra factor of 2, because in the DY process, one has to 
sum over the transverse polarizations of the photon, rather than average, as in DIS. 

Eventually, we are in a position to calculate the DY cross section in hadron-proton collisions, relying 



on the elementary cross section Eq. (109). According to the definition (106) the photon carries away 



the LC momentum fraction x\ from the projectile hadron, therefore the hadronic cross section reads, 
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(117) 



Employing in (109) the universal dipole cross section in the parametrization [87J, which we have 



already used in DIS, one can perform a numerical evaluation in a parameter free way. As an example, 



we compare in Fig. 29 the results of numerical evaluation of the DY cross section, plotted by solid 
curves, with data from the E772 experiment |130] at y/s = 38 GeV. To show the sensitivity to the 

(upper curve) and rrif 
500 GeV. 



200 MeV (lower 



quark mass in (112) the calculation was done with mj 
curve). The dashed curves show predictions for the energy of RHIC \fs 

Remarkably, the parameter free calculations explain well the observed cross section without any 
adjustment, and without any i^-factor. This is the advantage of the dipole approach: the phenomeno- 
logical dipole cross section is fitted to DIS data, I + p — y I' + X, at small Bjorken x which include all 
higher order effects, all possible gluon radiation. Therefore, one should not try to add to the above 
results higher order corrections, which would be double counting. This approach is also an alternative 
to the twist decomposition, since all higher twist effects are included by default. A numerical test of 
the dipole approach for Drell-Yan reaction vs NLO parton model calculations was performed in p.35j, 
and no visible deviation was found. 

As we already noticed, the DY process provides more observables than DIS, in particular the di- 
leptons are produced with a certain Feynman xf and transverse momentum q?. The xp distribution 
was calculated above in (117). The dipole approach also allows to calculate the differential DY cross 
section. The qx dependence comes from the interference of the DY amplitudes with different impact 
parameters and has the form |131[ 1136] , 

da(qN 7 *X) 



dlnad 2 qx 



1 

(2tt) 2 
1 



d pid p 2 e 



a"Jap 1 ,x 2 ) + <JqJap 2 ,x 2 ) - o"^(a|pi - p 2 \,x 2 ) 



-^#;. 9 (a,pi)tf 7 . 9 (a,p2) 



N i 



'118) 



Replacing the (^-integrated quark-nucleon cross section da(qp — > 7*A)/c?lna in (117) by the dif- 



ferential cross section (118), one gets the differential g^-dependent DY cross section in hN collisions. 
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4.2 Drell-Yan process on nuclear targets 



4.2.1 Shadowing in di-lepton production 

The total cross section of the DY process on a nuclear target can be calculated with the same equation 



(117) with a replacement da(qp — > 7*X)/dlna =>- da(qA — > 7*X)/c?lna. 
quark- nucleus DY cross section is also rather straightforward. One can employ Eq. (76) 



The calculation of the 
which was 

used for DIS on nuclei, replacing the photon distribution functions \l/^ L (a,r T ) by the 7*g distribution 

Examples of numerical results are presented in Fig 



functions ^l L q {a,r T 
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given by Eq. (Ill) 

^-dependence of the g^-integrated DY cross section is predicted [11 lj for the energies of RHIC and 
for proton and deuteron collisions with different nuclei. Calculations are done in the "frozen-dipole" 
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Figure 30: Shadowing for the total DY cross sec- 
tion as function of A 1 ^ 3 at M = 4.5 GeV and 
xf = 0.5. The lower pair of curves is for deuteron 
- gold scattering and the two upper curves are for 
proton - gold collisions. 
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Figure 31: Shadowing for the total DY cross sec- 
tion in proton - gold (upper curves) and deu- 
terium - gold (lower curves) collisions at the en- 
ergies of RHIC and LHC at xp = 0.5 as function 
of di-lepton mass M 2 . 



approximation which is sufficiently accurate at these energies, and xf = 0.5. The scale dependence 
of these calculations is presented in [3TJ As in DIS, the magnitude of shadowing slowly decreases with 
rising scale. Notice that gluon shadowing is important at these high energies and was included in the 
calculations (see details in |lllj ). 

The results demonstrate a considerable difference between the nuclear effects for proton and deuteron 
beams, which is related to isotopic effects. The radiation of a heavy photon is sensitive to the quark 
charge. 



4.2.2 Nuclear modification of the transverse momentum distribution 



The quark-nucleon differential DY cross section, Eq. (118) is easily converted to the quark-nucleus one 



if the conditions for the "frozen" approximation are fulfilled. Then, the cross section reads |131l 1 1 1 l~j . 



da(qA ->• 7*X) 
d\nad 2 q T 



x 



p J d 2 b J ^ft^e*^-^*;,^,^)^^,^ 

1 _ e -^(a Pl ,x 2 )T A (b) _ e -^K(a P2 ,x 2 )T A (b) + e -^(a\^-p 2 \,x 2 )T A (b) 



(119) 
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Convoluting this cross section with the quark distribution function of the beam hadron, one arrives at 
the differential hadron-nucleus DY cross section, 



dM 2 dxFd 2 qT 3ttM 2 x\ + X2 



da 



a 



F 2 h 



da{qA -> j*X) 
din a d?qT 



(120) 



j-i 



The numerical results of calculations of the nuclear ratio (JI]) with equations (117) and (120), at the 
energies of RHIC and LHC are depicted by dashed curves in Fig. [32} The ratio exhibits the so called 
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Figure 32: The influence of gluon shadow- 
ing on the nuclear ratio Eq. Q of the DY 
differential cross sections. Dashed curves 
are calculated without gluon shadowing, 
while solid curves include gluon shadow- 
ing. The influence on the longitudinal DY 
cross section is shown separately in the two 
left-hand plots (L). The two plots on the 
right show the DY ratio for the sum of 
the transverse and longitudinal cross sec- 
tions (T + L). Calculations are done at 
xp = 0.625, M = 6.5 GeV, and the ener- 
gies of RHIC and LHC. 
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Cronin effect (e.g. see Cronin effect for hadrons |137j ). a nuclear enhancement at medium values of qr- 
This can be interpreted as a result of increased mean transverse momentum of the di-lepton produced 
on a nucleus compared to a proton target. Intuitively is can be understood as a result of multiple 
interactions of the quark which "shakes off" a heavy photon with increase qp after getting a stronger 
kick from multiple rescatterings. 



Gluon shadowing was also added in these calculations, and the results are plotted in Fig. [32] by 
solid curves. One can rely on the small-p approximation, (Jq q (p) oc p 2 , which is quite accurate for DY 
processes at qx > 1 GeV [133]. In this approximation the dipole cross section is proportional to the 
gluon density |138j . 

7T 2 / A \ o „ / A 



<C0>> x ) . = - a s[ j 2 j P 2 G N [x, - ) , ( 121 ) 



where A ~ 10 is a phenomenological parameter. On a nuclear target one should replace Gjy (x, 
G a /A, and correspondingly in Eq. (119), 



<r»{p, x) <4(p, x)R G (x,^,b) , (122) 



The gluon contribution is appreciable, especially at the energies of LHC, and it practically terminates 
the Cronin enhancement for radiation of transversely polarized photons. 

4.2.3 A word of caution: incoherent effects imitating coherence 

In the case of DIS we saw above that the smaller x is, the longer is the coherence time t c ~ (2xmjy) _1 . 
Correspondingly, shadowing increases with 1/x. This is not true anymore for DY processes. According 
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to (107) x 2 reaches the minimal value when x\ — » 1. This is the kinematical limit and one needs to be 



cautious. The coherence length (same as time) for a fluctuation of a quark g — >■ ZZ g is given by, 

2E n 



where 



Mf lq 



M 2 



1 
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m 



'if 
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a 



a 1 



a) 



(123) 



(124) 



Apparently a > xi, therefore when x\ — >■ 1, then a — > 1 and t c — > 0, contrary to the naive expectation 
based on the experience with DIS. 



The mean coherence length can be evaluated similar to how it was done in Sect. 3.2.1 for DIS, i.e. by 



averaging l c weighted with the 7*g distribution functions and dipole cross sections. The results obtained 



in |139l 1 140] for pp collisions at 800 GeV are plotted by solid curves in Fig. 33 Although at small x\ 
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Figure 33: The mean coherence 
length (123) vs X\ at different 
dimuon masses. Dashed lines show 
the predictions of QCD factorization 
dependent only on x 2 . 



Figure 34: Shadowing in DY reaction on carbon, iron and 
tungsten as function of x 2 for M = 4.5 GeV (left) and 
7.5 GeV (right). Nuclear shadowing disappears at large 



and small x 2 because the coherence length, Eqs. (123) and 



(124), vanishes in these limits. 



{lc ) rises, because x\ is decreasing, upon approaching the kinematic limit X\ — V 1 is steeply falls down. 
Correspondingly, no shadowing is possible in this regime. Dashed lines linearly rising with 1/x 2 oc X\ 



correspond to the usual behavior (70) for DIS, and would be valid also for DY, if QCD factorization 
were correct. 

As a result of these kinematical constraints the interval of x 2 available for shadowing in DY data of 
the E772 /E866 experiments is rather narrow and the magnitude of shadowing [140] depicted in Fig. 34 
is considerably smaller that what one would expect based on QCD factorization. 

As long as shadowing in DY processes is weak at large x\, one may wonder what causes the con- 
siderable nuclear suppression observes at x% — Y 1 in DY data of the E772 and E866 experiments? 
An example of data is depicted in Fig. 35 The data are selected with large di-lepton invariant mass 



M > 6 GeV. According to Fig. [33] the coherence length for such a large scale is shorter than the mean 
nucleon spacing 2fm at all x±, therefore shadowing is impossible, while the data demonstrate a strong 
nuclear suppression steadily growing towards x 1 — 1. 

Qualitatively, one should expect any hadronic system propagating through a medium to dissipate 
energy |1391 11401 1142] . Apparently, this makes it more difficult to produce any particle, in particular a 
di-lepton, carrying a large fractional LC momentum x% — > 1, i.e. one should expect a nuclear suppression 
of leading particles. In terms of the Fock state decomposition this can be seen as nuclear enhancement 
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Figure 35: Ratio of DY cross sections on 
tungsten to deuterium as a function of x\. 
Data |141j are selected to have large di- 
lepton masses to eliminate nuclear shad- 
owing. The curves are calculated in [142] 
including energy loss effects. 



of higher components containing more gluons |142j . Correspondingly, should expect an xi-scaling for 
the nuclear suppression. The results on numerical evaluation [142] is plotted in Fig. 
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Notice that these results of the E772 and E866 DY experiments are included in most of the global 
DGLAP based analyses aiming at the extraction of the nuclear PDFs from data. The incorrect treatment 
of these DY data as shadowing certainly causes a bias in the results of such analyses. 

The effect of increase nuclear suppression at the edge of the kinematic range is a general phenomenon 
and has been also observed in hadron-nucleus collisions (see discussion and data in |143j ). 

The differential cross section of the DY reaction was found above to be modified by the coherence 
effects on nuclei. However, even if the coherence length is short, a similar modification is possible. Fig. 36 
presents data from the E866 experiment |144|. 1145] . which exhibit a pronounced Cronin effect, although 
the coherence length is too short to explain this by coherence. Even if l c is short, the projectile quark 
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Figure 36: Comparison of theoretical pre- 
diction of Rw/Beipr) vs. Pt (in GeV/c) 
to experiment for X2 = 0.05. Data 
are from the FermiLab E772/E866 experi- 
ments |144[ 1145] . The curve shows the re- 
sult of the calculation of broadening in the 
color dipole approach [146J. 



experiences multiple interactions in the nucleus and p^-broadening prior radiation of the di-lepton. The 
increase pr of the quark is transferred to the radiated dilepton as qx = apx- The quark broadening can 
be also calculated in the dipole approach [1361 1147] 



Ap 2 T 



2T A (b) 



(125) 



r T =0 



The result of numerical evaluation |146j of the Cronin effect is compared with data in Fig. 36 



5 DIS of neutrinos 



5.1 Parton model 



The charge current (CC) neutrino interaction v{y) + N — > 1(1) + X probes three different types of 
structure functions [148J, 



da u{D) -aa 2 m N E v 



dxdy 2sm\d w ){M 2 v + Q 2 ) 2 



i/ 2 o T?v{V)f , ft X7Tl N \ v (p), 
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± 



y 



2 



(x,Q 2 



(126) 



where the sign ± corresponds to either neutrino or antineutrino beams, respectively. 

The new interesting features comes with the structure function F$, which in the leading order 
carries information about the number of valence quarks N v = 3 in the nucleon. Namely, measuring the 
difference between the neutrino and antineutrino cross sections one arrives at a relation given by the 
Gross-Llewellyn Smith sum rule [149, 148J, 



dx 

x 



{xFf 



xFl N ) = N v 



TV 



c_ 
Q 2 



o(g- 4 ) 



(127) 



The current value of the left-hand side of Eq. (127) measured in the CCFR experiment |150[ I151j 
is 2.5 ± 0.18 ± 0.078, which is 1.8 standard deviations below the expected value |148j . Notice, however, 
that the CCFR data were taken on an iron target, therefore the structure function F% is subject to 
nuclear shadowing, reducing its value. 



5.2 Shadowing: dipole description 



As usual, the dipole representation is the most effective tool for prediction of nuclear shadowing. The 
PDFs describing light and heavy flavors at small x in the nucleon correspond to fluctuations of the W- 
boson into quark- ant iquark dipoles, which interact with the target. The Cabibbo enhanced transitions 
are to light quark dipoles, or strange-charm pairs. The perturbative light cone distribution functions, 



compared to Eqs. (77)-(78), also include the axial parts, and read |152j . 
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Using these distribution functions and the eikonal form (Eq. (28)) of the nuclear amplitude, the 
shadowing effects for neutrino-nucleus interactions were calculated in |1521 1153] . The results of [153J 
for 56 Fe are depicted in Figs. 37 and 38 for the nuclear ratios of F2 and xF^ respectively, plotted vs Q 2 
at fixed Bjorken x. Nuclear shadowing for F3 demonstrates a remarkably weak scale dependence. 

Unfortunately, none of these calculations is sufficiently realistic to be compared with data: (i) At 
small x < 10~ 2 gluon shadowing provides substantial corrections. It reduces the cross section of neutrino 
interaction on bound nucleons, and at the same time make the nuclear medium more transparent, (ii) 
At larger x > 10 -2 gluon shadowing can be neglected (see Sect. 3.2.6), however the "frozen" dipole 



size approximation employed in th e above calculations is not valid, and should be replaced by the path 
integral technique (see Sect. 3.2.2). (iii) At low Q 2 < 1 GeV 2 the perturbative distribution functions 
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(128 )-( 129 ) are not to be trusted. Nonperturbative distribution functions were derived for vector 
and axial [3H1 HO] currents in the instanton vacuum model. 

Notice that the expected magnitude of nuclear suppression, xF^/AxF^ « 0.85 explains well the 
observed deviation from the Gross-Llewellyn Smith sum rule, Eq. (127). Power and quark mass correc- 
tions to this sum rule were discussed and evaluated in [154J. The dipole description is an alternative 
way of calculation, which is supposed to contain all such terms. 



Part II 

Experimental Study of Shadowing in 
Lepton-Nucleus Scattering: Recent 
Developments 

6 Introduction 

The experimental study of shadowing and anti-shadowing in lepton nucleus scattering has been domi- 
nated by charged-lepton nucleus scattering. /i-A scattering at CERN and Fermilab and e-A scattering 
at SLAC provided early evidence that the value of F 2 (x, Q 2 ) per nucleon in nuclei was measurably dif- 
ferent than F 2 (x, Q 2 ) measured in deuterium in the low-x shadowing region. Most of these experiments 
were performed in the '80s and '90s, and there have been only a few such experiments, notably the 
HERMES experiment |155j . since the shutdown of muon DIS experiments at CERN and Fermilab and 
e-A experiments at SLAC. Several experiments at Jefferson Laboratory studied mid-to-high x nuclear 
effects in more detail, but did not contribute significantly to the study of shadowing. With the up- 
coming availability of the Jefferson Laboratory 12 GeV program offering a wider kinematic range, it is 
expected that charged-lepton nucleus studies of shadowing and anti-shadowing will begin again. 

On the other hand, studies of nuclear effects with neutrinos have only recently become possible 
with high-statistics DIS experiments that also provide detailed covariant error matrices, allowing the 
inclusion of correlated errors into the analysis. Neutrino nucleus scattering introduces a new facet 
in the study of shadowing (and all x-dependent nuclear effects) with the presence of the axial-vector 
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current. As shown in the theoretical treatment of shadowing, differences between I A and vA scattering 
in the shadowing region are expected. More detailed and varied programatic studies of nuclear effects 
with neutrinos are now being initiated with the availability of intense neutrino beams designed and 
commissioned for neutrino oscillation measurements. The need for this study is further reinforced by 
the necessity of understanding these nuclear effects in detail, to better estimate the systematic errors 
on measured neutrino oscillation parameters. These new studies of neutrino nucleus scattering will 
dominate this part of the review. 



7 Charged-lepton Nucleus scattering 

The present understanding of nuclear effects in lepton-nucleus scattering is mainly based on charged- 
lepton-nucleus (£A) DIS data. In the early 80s, the European Muon Collaboration (EMC) |156] found 
that the per-nucleon structure functions F 2 for iron and deuterium differ as a function of x B j orken = x. 
This intriguing result initiated a series of follow-up experiments [157, 1158^ I159[ I160[ 11611 I162[ \V6'i>\ 
ITM H651 ITUT1 HM IToTl H681 IMfl H7D] to investigate the nuclear modifications of this ratio, R[F$ A ] = 
F 2 A I (A F 2 N ), over a wide range of nuclear targets with atomic number A. These experiments established 
that when scattering off nucleons within a nucleus in the deep-inelastic region with Q 2 > 1 GeV 2 , the 
ratio of cross section per nucleon in nuclei to that in deuterium varies considerably in the kinematic 
range from relatively small Bjorken x ~ 10" 2 to large x ~ 0.8. The behavior of the ratio can 
be divided into four regions: 

• the shadowing region - R < 1 for x < 0.1, 

• the antishadowing region - R > 1 for 0.1 < x < 0.25 , 

• the EMC effect - R < 1 for 0.25 < x < 0.7, 

• and the Fermi motion region - R > 1 for x > 0.7. 

There is no single inclusive model explaining the nuclear modifications across the whole x region. 
The shadowing suppression at small x is the topic of this review, and the previous sections have 
summarized the current status of the theoretical understanding of shadowing. The anti-shadowing 
region is theoretically less well understood but, as indicated, might be explained by the application of 
momentum, charge, and/or baryon number sum rules. The modifications at medium x (the so-called 
"EMC effect" ) are still lacking a convincing, community-accepted explanation, but are usually described 
as nuclear binding and medium effects |171j . It has also been shown [172] that the EMC effect persists 
at lower Q and W in the resonance/transition region. The rapid rise of the ratio at large x is attributed 
to the Fermi motion of the nucleons within the nucleus. 

Notice that the currently available experimental data from charged lepton-nucleus DIS scattering 
have been obtained with fixed target experiments, where the values of x and Q 2 are strongly correlated. 
With the additional requirement that Q 2 be large enough in order to apply perturbative QCD (e.g., Q 2 > 
1 GeV 2 ), the majority of available experimental data covers a relatively limited region in x, where the 
effect of nuclear shadowing is far from saturation. The review by Geesaman, Saito and Thomas [171] 
summarized the experimental situation in 1995, and since then there has been limited amounts of 
new data on shadowing from charged-lepton nucleus scattering in the DIS region. Although limited 
recent data from charged-lepton nucleus scattering on shadowing are available, recent analyses of the 
shadowing region to extract nuclear parton distributions have been performed. Current work on this 



topic will be covered in Section 8.3 



With these qualifications, the evidence for nuclear effects in charged-lepton nucleus scattering can be 



summarized in Figure 39, which displays the F 2 Fe / ' F 2 D structure function ratio, as measured by both the 
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SLAC eA and the BCDMS fiA collaborations. The SLAC/NMC curve is the result of an A-independent 
parametrization fit to calcium and iron charged-lepton DIS data |157j 11611 1162^ \170\ \1TS\ I174j . This 
has often been used as a standard nuclear correction factor to convert data from a nuclear target to 
a free-nucleon target, for both charged-lepton and neutrino interactions. Although not displayed in 
this figure, there was also available experimental evidence for shadowing that is shown in the following 
figures. Notice also in those figures that, particularly in the shadowing region, a single A-independent 
curve does not reflect the experimental situation. 
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Figure 39: Nuclear correction factor, F^ e jF® , 



as a function of x. The parameterized curve is compared 



to SLAC and BCDMS data |157[ 1161] I162[ 1170} 1173] . The details of the shadowing regime are shown in 



the following Figure 40 and Figure 41 for other nuclei since Fe data in the shadowing region did not 
survive other applied cuts. 



A more detailed look at shadowing experimental results, for specific nuclei such as C, N, Ca, Ag, 

These data have been used in the comparison with the 



Sn and Xe, is shown in Figures 40 and 41 



parameter-free calculations shown in Figure [22] and also appear in Figure [12] Even for deuterium itself, 
data from /i-nucleus scattering in the experiments E665 at Fermilab and NMC at CERN show clear 



evidence of shadowing in the ratio of deuterium to hydrogen, Figure 42 
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Figure 40: Ratios of the deep inelastic cross section on targets of (a) carbon and nitrogen and (b) 
calcium to those of deuterium |170t 11561 H651 1 
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Figure 41: Ratios of the deep- inelastic cross section on targets of tin, xenon and silver to those of 
deuterium [1701 CESE1 H55] 




Figure 42: Nuclear correction ratio, F^/F^ 1 , as a function of x [§5"1 U69J. 



All experimental results in the shadowing region involving charged-lepton nucleus scattering are 
consistent with the onset of shadowing occurring at a value of the Ioffe time (55) roughly equivalent to 
the separation between nucleons in a given nucleus. Shadowing increases with increasing A and decreases 
with increasing Q 2 . The Q 2 -dependence in the shadowing region, a small but non-zero slope in log Q 2 , 
as seen by the NMC experiment, is consistent with descriptions |175] of shadowing prevalent some 



years ago, of low-Q scattering being ascribed to vector meson dominance, as illustrated in Figure 20 



An understanding of the role of the longitudinal cross section in shadowing has been provided by the 
HERMES experiment |155] that examined the ratio of Ra/ Rd with R = <jl / <Jt the longitudinal to 
transverse DIS cross sections. They found that Ra/ Rd is consistent with 1.0 for A up to 84 Kr and the 
kinematic range down to x = 0.01 and Q 2 = 0.5 GeV 2 . Furthermore, the HERMES data were quite 
consistent with earlier measurements by the NMC collaboration [166J 1176] accumulated at considerably 
higher energies. No significant Q 2 -dependence was observed over the wide Q 2 range covered by the 
combined HERMES and NMC experiments. 



8 Neutrino Nucleus Scattering 

The measurement of neutrino deeply inelastic scattering, besides providing valuable information on 
the weak current, can be a significant aid in determining the parton distribution functions (PDFs). 
Neutrino and antineutrino DIS data directly measure specific quark flavors including the strange quark 
and strange anti-quark distributions. Neutrino scattering provides the only current input to determine s 
and s in the proton via the strange-to-charm quark transition and the subsequent detection of a charmed 
particle decay. Combining neutrino Charged Current (CC) W ± interactions with other probes, such as 



48 



charged lepton (j,Z), allows us to disentangle the separate PDF flavor components. 

While for charged-lepton scattering the use of a nuclear target is a choice, for high statistics neutrino 
scattering experiments it is a necessity, due to the weak nature of the {W^, Z} interactions. This 
complicates the interpretation of neutrino scattering data, since these data are used for the extraction 
of free nucleon PDFs, which depend on nuclear correction factors to account for the nuclear effects of 
the nucleon being bound in a nucleus. 

As in charged-lepton nucleus scattering, these nuclear correction factors essentially convert the 
experimentally-measured structure functions in a bound nucleon within a nucleus to the corresponding 
structure functions within a free nucleon. For charged-lepton nucleus scattering these factors are usually 
given in terms of curves, such as in Figure 39, measured for a given nucleus at the relevant Q. In earlier 
QCD global fits that attempted to include neutrino DIS data, the charged-lepton nuclear correction 
factors were simply applied to neutrino nucleus scattering results. From the theoretical side it was 
thought that these factors would be the same for neutrino nucleus scattering as for charged-lepton 
nucleus scattering, with the possible exception of the shadowing/anti-shadowing region. 

Early attempts to include neutrino-nucleus DIS scattering data, corrected with charged-lepton nu- 
clear correction factors, introduced such tension in the shadowing region at low-x in global QCD fits 
(described below) that the low-x neutrino data was excluded early-on in CTEQ nucleon PDF global 
analyses. In a more recent examination of high-x parton distribution functions, carried out by the CTEQ 
collaboration [174] and |177j . indications began to accumulate that the nuclear correction factors for 
neutrino nucleus scattering could indeed be different than those for charged-lepton nucleus scattering, 
and not only in the shadowing region. 

Although we will be concentrating on the shadowing region, in this section we will describe current 
studies of nuclear effects in neutrino nucleus interactions across the full x range, since the behavior 
of the cross section as it approaches the shadowing region from higher x is important. Using recent 
neutrino- Fe and neutrino-Pb scattering data to directly extract the neutrino nucleus nuclear correction 
factors, we will compare them to those extracted from charged lepton scattering, as well as to specific 
model-dependent corrections provided by Kulagin and Petti [78J. The limitations of the Kulagin- Petti 
model have been emphasized in section |3.1| However, it is the only "model" that has been available 
for v nucleus scattering, and therefore it has been heavily used by experimental and phenomenological 
analyses. Finally, we will concentrate on what neutrino scattering data off nuclei is telling us about 
shadowing and anti-shadowing, in interactions involving the axial- vector current. 



8.1 Early Neutrino Nucleus Scattering Results 

Soon after the discovery that muon interactions off bound high-A nuclei exhibited a different x depen- 
dence than muon interactions off deuterium, this observation was tested with neutrino nucleus interac- 
tions. The earliest tests were dominated by bubble chamber results |178l 11791 11801 11811 1182] . with one 
notable exception involving the placement of a hydrogen target upstream of the CDHS (mostly) iron 
detector [T83] . 

None of the bubble chamber experiments was able to find a statistically significant effect such as that 
found by the EMC collaboration. The combined CERN WA25 and WA59 analysis of v and V interactions 
off deuterium and neon [179J concluded that their results were consistent with the EMC higher Q 2 
data for the x range between 0.3 and 0.6, and consistent with the SLAC lower Q 2 data |159l 1158] for all 
x. They also concluded that there was no Q 2 -dependence in the low-x shadowing region, in the range 
of 0.25 to 14 GeV 2 . However, within the shadowing region they did see a depletion in the ratio of Ne 
to D, as indicated earlier with Figure |6j Nevertheless, the statistics are limited and the claim of an 
overall normalization uncertainty of 1% between the two independent exposures of WA25 and WA59, 
based on the comparison of "muon" distributions in the alcoves, may be underestimated based on more 
recent understanding of the signal vs background in these alcove measurements. The main challenge 
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of these early v A experiments was the accumulation of sufficient statistics with bubble chambers to 
make significant comparisons between cross sections of different nuclei. An additional difficulty in the 
comparison across experiments was the > 20% uncertainty on the flux of incoming neutrinos for each 
independent experiment. As mentioned above, even attempts at comparison of experiments in the same 
beam were hampered by uncertainties in signal to noise in the measurements of muon alcove intensities 
used for normalization. 

A procedure was developed by the Fermilab E-745 heavy liquid bubble chamber collaboration |18 1] to 
use the bubble-chamber's ability to isolate nucleon recoils ("dark tracks") as indication of interactions off 
tightly-bound nuclei. By taking the ratio of events with and without these dark tracks, the collaboration 
found an x-dependence of this ratio similar to that found by the earlier-referenced EMC collaboration. 
However, a follow-up analysis by a BEBC collaboration [182J showed that one could get a similar 
x-dependent ratio distribution by applying the same criteria to neutrino scattering off both deuterium 
and hydrogen, thus suggesting this criteria was not indicative of scattering off a tightly bound nucleus. 

Although the conclusion of the CDHS study of the ratio of iron to hydrogen cross sections and struc- 
ture functions was that there were no statistically significant differences indicative of nuclear binding 
effects in iron, their results, figure la of their publication |183j . certainly follow the trend displayed in 
the more statistically powerful CERN muon and SLAC electron results. 

8.2 Recent Neutrino Nucleus Scattering Data: The ratio of v-K to v- 
Nucleon Cross sections 

The more contemporary study of v nucleus scattering using high-statistics experimental results with 
careful attention to multiple systematic errors began with the CDHSW, CCFR/NuTeV v Fe, the NO- 
MAD v C and the CHORUS v Pb experiments. Whereas the NuTeV pH] and CHORUS [US] Col- 
laborations have published their full data sets, NOMAD [186] has not yet done so. Although NOMAD 
displayed a small part of their preliminary data in a 2004 publication |187j . using their data would 
require extracting their measured points with errors from figures, an inaccurate procedure at best. The 
more detailed analyses presented in the following will be made with mainly the CDHSW, CCFR and 
NuTeV v Fe results and the CHORUS v Pb results. 

The most direct comparison between scattering neutrinos off nuclei and off nucleons is a comparison 
of the cross sections. This involves none of the assumptions that will be needed to extract structure 
functions from these cross sections nor the procedures needed for global fits. Members of the CTEQ 
collaboration have studied this direct cross section comparison |174l 1177] . The eventual goal of this 
study was to extract the nuclear correction factors from v nucleus scattering, in order to compare them 
with the similar nuclear correction factors from charged-lepton nucleus scattering, and to use them in 
incorporating v nucleus data in global analyses. This, in turn, requires taking a ratio of the measured v A 
cross sections to the corresponding "free-nucleon" quantities. For the charged-lepton nuclear correction 
factors, the measured per-nucleon cross sections or structure functions off nucleus A were compared 
directly to measured charged-lepton scattering results off deuterium. For neutrino scattering, a similar 
ratio using scattering off deuterium is impossible, since there exist only low-statistics bubble chamber 
data (see |188] and references therein). The statistical and systematic accuracy of these early bubble 
chamber neutrino deuterium results is not sufficient for determining a neutrino nuclear correction factor. 

Determining "free" Cross sections and Structure Functions for Comparisons To address 
this, the nCTEQ collaboration decided to approximate the deuterium denominator in the neutrino 
correction factors by creating "free" neutrons and protons from PDFs determined using a new reference 
fit that involved data off mainly hydrogen and deuterium targets. This reference fit starts with the 
CTEQ6.1M PDFs [54J and removes any neutrino nucleus data, keeping data involving proton and 
deuterium targets. Nuclear corrections for deuteron targets were included where appropriate. Details on 
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Figure 43: Comparison between the as measured Chorus and NuTeV neutrino cross section data off 
nuclei, without any nuclear corrections, and the reference free-nucleon fit ("Theory"). 



the treatment of heavy quarks, which primarily affects the low- a; region for neutrino induced processes, 
and therefore are significant when examining the shadowing region, as well as high-x corrections, can be 
found in reference |174j . The resulting data for the reference came from the BCDMS Collaboration [189J 
for F$ and F 2 d , from the NMC Collaboration [Mj for F% and from HI [190] and ZEUS [191] 

for F 2 P , from CDF |192j and D0 |193] for inclusive jet production, from CDF |194j for the W lepton 
asymmetry, from E-866 |195] for the ratio of lepton pair cross sections for pd and pp interactions, and 
from E-605 |196j for dimuon production in pN interactional 



The CTEQ determination of the ratio of v-Fe and z/-Pb to z/-Nucleon Cross sections In 

their analysis of vA scattering data, the CTEQ group used two of the recent sets of experimental 
results mentioned; the CHORUS v Pb and the NuTeV v Fe data sets, examining the behavior of the 



as measured cross sections displayed as a ratio to the "free-nucleon" cross section. Figure [43] shows 
the results with only statistical errors displayed on the measured points. In general, there is agreement 
between the different data sets, which is somewhat surprising. Aside from a few low-statistics CHORUS 
V data points, in the approach to the shadowing region from high-x (0.6 to 0.08) there is quite good 
agreement between v and V and between Fe and Pb. 

The CHORUS V Pb data displays a very "traditional" x-behavior with (from high- to low-x) an 
antishadowing region starting at x = 0.3, with a maximum and turnover toward the shadowing region 
at x ~ 0.15. On the other hand, the CHORUS v Pb data, with significantly smaller errors, have a 
much smaller antishadowing region shifted toward lower-x. Similarly, in approaching the shadowing 
region from high-x, the NuTeV V data follow the CHORUS v Pb data closely. The NuTeV v data, 
again having much smaller errors than the V data, are again consistent with both the CHORUS v Pb 
and the NuTeV V data in the approach to the shadowing region. 



However, in the low-x shadowing region the results are quite different, as shown in Figure 44 (low-x 



details of Figure 43). Only the lower-statistics CHORUS v and v cross section ratios turn over and 



4 Note that the E-605 data were taken on a copper target, but the nuclear corrections in the relevant kinematic region, 
outside the shadowing region, have been measured [198j to be consistent with A 1 , so no nuclear corrections were included 
for these data 
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Figure 44: Low-x details of the ratio of the Chorus and NuTeV measured cross sections, without any 
nuclear corrections, to the "free-nucleon" reference cross section. 



approach the expected R <1.0 axis. From the same figure, the NuTeV data sets show a somewhat 
different behavior. Only at the very lowest x = 0.015 point (with the largest error) does the NuTeV z7 
data turn over toward the shadowing region. The considerably more accurate NuTeV v data does not 
exhibit a shadowing turn-over even at the lowest x = 0.015 point. 

It is reasonable to ask if this observed behavior of v and V off Fe and Pb is consistent with expec- 
tations. One of the few models available that has been specifically designed to predict the v A nuclear 
corrections and is widely used by the experimental neutrino community is the model of S. Kulagin and 
R. Petti (K-P) |199[ I200[ |201j. Their approach includes a QCD description of the nucleon structure 
functions at high x, as well as the treatment of Fermi motion and nuclear binding, off-shell correction 
to bound nucleon structure functions, nuclear pion excess and nuclear shadowing. Their model involves 
several mechanisms in specific regions of x#j, to enable covering the entire x^- range. This model has 



limitations, particularly in the shadowing region, emphasized in 3.1 However, it predicts a shadow- 



ing turnover that is not seen in the NuTeV Fe data. How significant is this deviation of the NuTeV 
measurements from expectations can be determined by carefully considering the errors of the measured 
points. 

A way of understanding how flexible the CHORUS and NuTeV measured points are within their 
quoted errors is to apply a correction, like the K-P correction, to these cross section ratios, allowing the 
so-corrected central value to vary, with an applied penalty term for excursions, within the systematic 
errors. This then also serves as a test of the K-P model's predicted nuclear correction factors for v 
A scattering. CTEQ corrected all four data sets with the Kulagin-Petti factors, while bringing in the 
systematic errors of each data set into the fit. Within the larger systematic errors of the CHORUS 
data, the K-P nuclear correction factors successfully correct the CHORUS data, bringing the corrected 



data to the R = 1.0 axis (within errors), as shown in Figure 45, On the other hand, contrary to 



expectations the uncorrected NuTeV data ratios are already at values of R > 1.0 with small errors in 
the "shadowing" region. The K-P correction factors, being modeled to raise the values of R (toward 
1.0) in the shadowing region, force the NuTeV values even further away from 1.0 when applied to the 
NuTeV data points. This is a very different result than the one obtained in £Fe scattering. 



Evidence for Shadowing from v A Cross sections We then conclude that the CHORUS data, 
within their errors, are consistent with the predictions of a more conventional model in the shadowing/ anti- 
shadowing region. The NuTev data is certainly NOT consistent with conventional expectations in the 
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Figure 45: As an indication of the size of their systematic errors and the effect of the K-P nuclear 
correction factors, the low-x details of the comparison between the Chorus and NuTeV neutrino cross 
section data and the free-nucleon reference cross sections are shown, after the K-P corrections have 
been applied and the measured points have been allowed to vary within the quoted systematic errors 
to reduce the x 2 ( see text). 



shadowing/antishadowing region and the NuTeV v cross section data do not show any shadowing be- 
havior at all, down to the lowest x value of x = 0.015. 

8.3 Recent Neutrino Nucleus Scattering Data: Global QCD Fits and Nu- 
clear Parton Distribution Functions 

To appreciate the more recent analyses of nuclear effects in neutrino nucleus scattering, it is necessary 
to understand the concept and application of global QCD analyses of various experimental data. It has 
been through such global fits, attempts to find a single fit-solution describing multiple data sets, that 
tension between neutrino nucleus scattering and other data sets have suggested that nuclear correction 
factors applied to the neutrino scattering data need to be re-examined. 

The QCD Parton Model provides a comprehensive framework for describing general high energy 
processes. In this framework, the cross-section for a hadron-hadron collision process is written as the 
convolution of a set of universal non-perturbative Parton Distribution Functions (PDFs fi{x,Q"j)) and 
a perturbatively-calculable hard cross section da. 



da i+N->i>+x = J2fi N (Q 2 )®da l+ ^ l ' +i '(fi 2 ,Q 2 ). (131) 

i=q,q 

The parton distributions describe the structure of hadrons, and are the link between the physically 
measured cross-sections da l+N ^ 1 +x and the basic processes of the theory. Obviously, the precise deter- 
mination of these functions is of fundamental importance for the interpretation of experimental results 
within the Standard Model. These parton distribution functions should be applicable to all processes. 
A goal of global QCD fits is to extract a set of parton distribution functions, while testing this concept 
of universality. It is this concept of universality that will be considered in the discussions that follow. 

Although global analyses began in the late 1970's with the work of Gluck and Reya |202] , and con- 
tinued with the analyses of Duke and Owens [203J and Martin, Roberts and Stirling [204J, it wasn't until 
the global analyses of Morffn and Tung |205j . which were initiated in the 1988 Snowmass Workshop [206J, 
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that increased emphasis was placed on careful consideration of experimental systematic errors. Since 
then there have been many subsequent global analyses, accommodating an increasing ensemble of types 
of data with increasingly sophisticated treatment of experimental errors. The leading groups working on 
global analyses of free-nucleon PDFs today are the CTEQ (works described in this review and |207j ). 
the MSTW [208], and the NNPDF [209J collaborations. There has been no reason to suspect that 
independent global analyses could not be applied to experimental data off nuclear targets as well as 
nucleon targets. The challenge has been to perform global fits to a combination of nucleon-target and 
nucleus-target data, which call for nuclear correction factors to compensate for nuclear effects. The 
results of fits to nuclear targets and the determination of nuclear correction factors will be discussed 
below. 

Global analyses adopt functional forms appropriate for parton distributions at all Q, so the evolved 
distributions can then be given in a simple analytic form, which varies smoothly over the entire range of 
x. The goal is to minimize the role of the chosen initial point of evolution. A simple form that satisfies 
the above considerations is, for example: 



where the c-coefficients are called the shape parameters. An advantage of this parametrization for a 
study of shadowing is that it provides a way to emphasize the small-x behavior of the parton distributions 
via the value of c\. 

Once the parameterization of the PDFs has been chosen at the initial scale Qq as indicated above, 
the procedure for a "global fit" is then: 

• Evolve the PDFs from the initial Q$ to the relevant larger scales Q 2 > Qq of the experimental 
data by the DGLAP equations. 

• Compute the physical observables to be tested. 

• Determine how well the computed observables match the experimental ones via a x 2 -function that 
compares the measured experimental values and their associated statistical and systematic errors 
with the computed values. 

• Vary the c-coefficients until the best agreement with the data is reached by observing a minimum 
in the x 2 -function. 

The actual process of finding the minimum x 2 i n practice is a non-trivial task, since the x 2 involves both 
statistical and systematic experimental errors. General-purpose minimization packages like MINUIT [210J 
can be used to initiate the process, but global fitting is also an art that requires additional considerations. 

8.3.1 Global Fits: Nuclear Structure Functions, Nuclear Parton Distributions and Cor- 
rection Factors 

A comparison of results from neutrino scattering as well as charged-lepton scattering off nuclear targets 
can be done with global fits to extract nuclear parton distributions. With these nPDFs, one can then 
study both the magnitude and the shape of the nuclear correction factors Rf(x, Qo), for both neutrino 
and charged-lepton scattering. The experimental input can be, preferably, the measured differential 
cross sections or the extracted nucleon/nuclear structure functions. 

Using the ratio of extracted u-A to z/-Nucleon Structure Functions A possible experimental 
input to these global fits involves the extracted nuclear structure functions from the differential cross 
sections. With sufficient statistics and control of systematic errors the structure functions F 2 (x,Q 2 ), 




(132) 
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xF 3 (x, Q 2 ) and, Rl(x, Q 2 ) = o L jo T can be determined independently for neutrino and antineutrino, by 
fitting the y-distribution of events in a given x-Q 2 kinematic "box" . This procedure leads to a total of 
six structure functions. 

Somewhat less demanding in statistics and control of systematic errors, the "average" structure 
functions F 2 (x,Q 2 ) and xF 3 (x,Q 2 ) can be determined from fits to combinations of the neutrino and 
antineutrino differential cross sections, with several assumptions. This is the procedure used by the 
NuTeV collaboration. The sum of the v and V differential cross sections, yielding F 2 can then be 
expressed as (see also 126): 



cPa v d 2 a v _ G 2 F ME 



dxdy dxdy n 



A Mxy y 2 l + AM 2 x 2 /Q 2 \ T? / y\ 

2 { 1 ~ y -^E + 2 1 + R L ) F2 + V V~ 2) AxF3 \ ■ (133) 



where now F 2 is the average of F% and F%, and the last term is proportional to the difference in xF 3 for 
neutrino and antineutrino probes, AxF 3 = xF 3 — xF 3 . At leading order, assuming symmetric s and c 
seas, this is 4x (s — c). The cross sections are also corrected for the excess of neutrons over protons in 
the iron target (5.67%), so that the presented structure functions are for an isoscalar target. 

A significant step in the determination of F 2 (x, Q 2 ) in this manner, that surely affects the low-x 
(shadowing) values, is the assumed AxF 3 and Rl(x,Q 2 ). NuTeV uses a NLO QCD model as input 
(TRVFS [21 lj ) . and assumes an input value of Rl(x, Q 2 ) that comes from a fit to the world's charged- 
lepton measurements [212J. It is important to emphasize that the NuTeV F 2 (x,Q 2 ) extracted points 
in the shadowing region with x < 0.1 are extremely dependent on these two assumptions. In addition, 
from the NuTeV publications it is difficult to determine if the assumed R and AxF 3 have been corrected 
for nuclear effects up to Fe and, if so, how has this been done. For this reason it is preferable to directly 
use the measured differential cross sections in these fits rather than the extracted average structure 
functions. 



Determining the Nuclear Correction Factors There are two alternative approaches for deter- 
mining these nuclear correction factors. They both lead either directly or indirectly to the concept of 
nuclear parton distribution functions. One approach that will be referred to as method 1, is to use 
the basic global analysis formalism described above, implemented in a global PDF fitting package of 
experimental results off nuclei, with no nuclear corrections applied to the analyzed data. If the nuclear 



modifications of the x-shape in this analysis are modest, the initial parametrization of Equation 132 



should provide enough flexibility to also accommodate these now nuclear parton distribution functions 
(nPDF). Notice that the resulting nPDFs are frequently referred to as the PDFs of a bound nucleon. 
Such an interpretation is reasonable at large x, because the lifetime of hadronic fluctuations is short, 
and the incoming lepton reaches a bound nucleon deep inside the nucleus without any preceding inter- 
actions. However, this is not true at small x in the shadowing region. A long living hadronic fluctuation 
coherently interacts with a whole chain of nucleons, rather than with a separate bound one, therefore 
the DIS process probes the collective PDFs of the nucleus. These nPDFs are different from the ones of 
a free nucleon, and are subject to shadowing, i.e. they rise with the number of nucleons slower than lin- 
early. This can also be interpreted in the infinite momentum frame of the nucleus, where as is explained 
in Sect. 3.1 and illustrated in Fig. 10, the parton clouds originating from different nucleons overlap and 
can fuse at small x, where they form a common parton distribution for all involved nucleons. These 
collective nPDFs cannot and should not be identified with any individual bound nucleon. 

Given then both nPDFs and free-nucleon PDFs, the nuclear correction factors can be determined for 
a given nucleus. Within the parton model, the nuclear correction factor Rf(x, Qo) for an experimentally 
measurable quantity Mj can be defined as follows: 
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rA(x x _ Mf(x,Q ) 



where Mf-(x,Qo) represents the measurable quantity computed with nuclear PDFs, and M^(x,Q ) is 
the same quantity constructed for the average nucleon, (n+p)/2, out of the free-nucleon PDFs. 

An alternative approach, method 2, is to assume a set of free-nucleon PDFs, calculate f£(x,Qo) 
and f™(x,Qo), and then use global fitting techniques to fit the nuclear correction factors, also called 
weighting-factors that modify the free nucleon PDFs to provide agreement with the measured nuclear 
quantity. This method also yields nPDFs as a by-product of the fit: 

M A (x, Q Q ) = R?(x, Q ) j [Z fi(x, Q ) + (A - Z)f?(x, Q )) , (135) 

In either method, Rf(x,Qo) depends on the observable under consideration since different observables 
may be sensitive to different combinations of PDFs. For example the nuclear correction factor Rf(x, Qo) 
for FA and F% will, in general, be different. Additionally, the nuclear correction factor for F A will yield 
different results for the u-Fe process as compared with the charged-lepton £ -Fe process. In both cases 
this is due to the different parton contributions to each process. 

Contemporary attempts to use global fitting techniques to extract these nuclear parton distribu- 
tion functions and/or the nuclear correction factors began with the pioneering work of Eskola and 
Paukkunen |213[ 12151 1214j and later in collaboration with C. Salgado |216j . The study continued with 
the works of Hirai, Kumano and Nagai |217[ I218[ 1219] and deFlorian and Sassot [2201 1225j . all who 
used variations of method 2. The latest group to compute nPDFs, the nCTEQ group, employs method 
1 [2211 12221 12231 1224] . 



8.3.2 Global Fits: The nCTEQ Nuclear Parton Distributions and Correction Factors 

In contrast to the procedure employed by other nPDF global analysis groups, the nCTEQ group0 
employs method 1, that is to use the basic global analysis formalism, implemented in a global PDF 
fitting package of results off nuclei, with no nuclear corrections applied to the analyzed data and no 
assumed set of free-nucleon PDFs. This then is a direct fit of the nPDFs themselves. 

To confirm that this procedure that uses nuclear global fitting techniques and does not assume a set 
of nucleon PDFs to start with, is reasonable, the procedure is first applied to charged-lepton scattering 
results. The goal of the nCTEQ group is to determine an independent evaluation of the charged-lepton 
nuclear correction factor, F^/F^ , that is preferred by the data. It can than be directly compared with 
the SLAC-NMC factor and with the predictions of other models and global charged-lepton fits. After 
cuts on Q (> 2.0 GeV) and W (> 3.5 GeV), the fit used over 700 data points over the complete x range, 
which are listed in Tables I, II and III of publication |222j . The nuclei included in the fit ranged from 
He to Pb in ratios to deuterium, and Be to Pb in ratios to C. In particular the data surviving the cuts 
in the ratios of C, Ca and Cu to deuterium were the main drivers to the turnover into the shadowing 
region at low x. There were no Fe to deuterium ratios for x < 0.7 (the shadowing region) that survived 
the cuts. 



The nCTEQ result is displayed in Figure (46 -a) for Q 2 = 5 GeV 2 , and in Figure 46 -b) for Q 



20 GeV 2 . The figures of the ratio F^JF® are in good agreement with the SLAC/NMC parameterization 



of Figure 39 and with the fits from Hirai-Kumano-Nagai (HKN07) [219J. The data points displayed in 



5 This section is a summary of work performed by a CTEQ/Grenoble/Karlsruhe collaboration consisting of: K. Kovarik 
(Karlsruhe, Inst. Technology), I. Schienbein (Joseph Fourier U.), J. Y. Yu (Southern Methodist U.), C. Keppel (Hampton 
U), J. G. Morfi'n (Fermilab), F. I. Olness (Southern Methodist U.), J. F. Owens (Florida State U.) and T. Stavreva (Joseph 
Fourier U.) 
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Figure 46: Using a fit to all charged-lepton data (i^A) and DY data, off a variety of nuclei, that passed 
their cuts, the nCTEQ computed nuclear correction factor, 



F 2 ^ 



/Ff (solid line labeled fit B) 



SIS db 

function of x is shown in Figure-(a) for Q 2 = 5 GeV 2 and Figure-(b) for Q 2 = 20GeV 2 . Both fits are 
compared with the SLAC/NMC parameterization, as well as the Hirai et al. (HKN07), (Ref. [219]) fit 
to charged-lepton nucleus scattering data and the Kulagin-Petti model predictions. The data points 
displayed in the Figures are the measured Fe to D ratios that survived the cuts. 



Figure [46] are the data available off Fe that satisfied the cuts. The success of this comparison suggests 
that this method can also be successfully applied to u-A scattering results. 



Nuclear Correction Factors for F%(x,Q 2 ) and F^x^Q 2 ) A similar procedure is then used to 
extract the neutrino nuclear correction factor F 2 Fe /F 2 D . To apply this procedure now to u-A scattering, 
there are several data sets that can be considered. The earliest is the CDHSW u-Fe data, followed 
by the CCFR u-Fe data, the NuTeV u-Fe data and finally the CHORUS u-Pb data. However, if one 
demands a full covariant error treatment of the published data, yielding maximal discriminatory power 
of the data, there is essentially only the NuTeV u-Fe and CHORUS u-Pb scattering data for input. This 
would also be a natural selection by the fit itself. If the CDHSW and CCFR data, with their errors 
calculated via the sum of the squares of statistical and systematic errors, were added to the NuTeV and 
CHORUS data with their full covariant error matrix for the fit, the weight of the CDHSW and CCFR 
data points in the combined fit would be greatly reduced. This was shown in the earlier treatment of 
neutrino cross sections, when the NuTeV Fe and the CHORUS Pb data were combined. In this case, 
even though both data sets have full covariant error matrices, the relatively small NuTeV errors with 
respect to the CHORUS errors enable the NuTeV data points to dominate the combined fit. 

An additional input to the fit is the NuTeV/CCFR di-muon data (226J, which are sensitive to the 
strange quark content of the nucleon. However, no other data such as charged- lepton-nucleus {(^A) 
and DY data have been used. Because the neutrino data alone do not have the power to constrain all 
of the PDF components, a minimal set of external constraints [224J also has to be imposed and some 



of these external assumptions do indeed affect the behavior of the fit parton distributions at small x 
- the shadowing region. These include the Callan-Gross relation [F 2 A = 2xF" A ) as well as use of the 
assumption s = s and c = c. 

It is important to note that the nCTEQ fit is made directly to the double differential cross sections in 
order to extract the set of nPDF consistent with the NuTeV measured differential cross section points. 
The fit does not use the extracted NuTeV results of the average value of F 2 (x, Q 2 ), which contains all 
the assumptions made by NuTeV to extract them. The extracted nPDFs are then combined with the 
earlier-described free-nucleon PDFs to form the individual values of the nuclear correction factor R for 
a given x and Q 2 , separately for neutrino and anti-neutrino - not the average of both - as shown in 



Figure 47 for u— Fe and in Figure 48 for u-Fe. Since the difference between F 2 (uA) and F 2 (uA) is small, 



57 




(a) (b) 

Figure 47: Nuclear correction factor R for the structure function F 2 in charged current vFe scattering 
at a) Q 2 = 5 GeV 2 and b) Q 2 = 20 GeV 2 . The solid curve shows the result of the nCTEQ analysis of 
NuTeV differential cross sections (fit A2), divided by the results obtained with the reference fit (free- 
proton) PDFs; the uncertainty from the A2 fit is represented by the yellow band. Plotted also are NuTeV 
data points of the average F 2 (see text) to illustrate the consistency of the fit with the input points. 
For comparison the correction factor from the Kulagin-Petti model [75] (dashed-dot line), HKN07 [219J 
(dashed-dotted line), and the SLAC/NMC parametrization (dashed line) of the charged-lepton nuclear 
correction factor are also shown. 




the consistency of these extracted values of F 2 (x,Q 2 ) with the measured average values from NuTeV 



can be shown in Figure 47 



It is also of course possible to combine these fitted nPDFs to form the individual values of the 
average of F%(yA) and F 2 {i>A) for a given x and Q 2 , to compare directly with the NuTeV published 
values of this quantity. This was recently done and the nCTEQ preliminary results^] for low-Q 2 are 



shown in Figure 49 



The comparison between the nCTEQ fit that passes through the NuTeV measured points, and the 
charged-lepton fit is very different in the lowest-x, lowest-Q region and gradually approaches the charged- 
lepton fit with increasing Q. However, the slope of the fit approaching the shadowing region from higher 
x, where the NuTeV measured points and the nCTEQ fit are consistently below the charged-lepton A 
fit, make it difficult to reach the degree of shadowing evidenced in charged-lepton nucleus scattering at 
even higher Q 2 



6 We thank Karol Kovarik, (Karlsruhe, Inst. Technology) for allowing us to use these preliminary results presented at 
the 2012 DIS conference in Bonn 
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Figure 49: Nuclear correction factor R for the average F2 structure function in charged current vFe 
scattering at Q 2 =1.2, 2.0, 3.2 and 5.0 GeV 2 compared to the measured NuTeV points. The green 
dashed curve curve shows the result of the nCTEQ analysis of v A (CHORUS, CCFR and NuTeV) 
differential cross sections plotted in terms of the average F^ e divided by the results obtained with 
the reference fit (free-proton) PDFs. For comparison, the nCTEQ fit to the charged-lepton data from 
Figure 46 is shown by the solid blue curve. 



nCTEQ: Shadowing in i/-A Scattering As concluded above, the nCTEQ fit to t^A and DY data 
is generally consistent with the SLAC/NMC parameterization and HKN fits to calibrate and certify 
the nCTEQ method. However, the nCTEQ fits vFe scattering in Figures 47, 48 and 49 do not agree 
with the charged-lepton nucleus fits. Although there is a general tension over most of the x-range, the 
biggest difference is in the shadowing region, where the non-existent turn-over of the NuTeV vFe cross 
section data in the shadowing region at lower Q 2 is again evident. Since lower Q 2 data dominates the 
low-x region, this also explains the lack of any shadowing turnover at all in the ratio of cross-sections 
in Figure '43 



8.3.3 Global Fits: The Eskola-Paukkunen-Salgado (EPS) Nuclear Parton Distributions 
and Correction Factors 

Although the discussion of nuclear parton distributions began in the late 70's |228j . it was with the 
early work of K. J. Eskola and colleagues [229] that the more modern era of extracting nuclear parton 
distributions via global fits to a variety of lepton-nucleus data began. Their systematic consideration of 
uA scattering was treated in their work of 2009, published in 2010 [230], and then revisited to specifically 
address the consistency of vA DIS results and general nPDF results [23 lj. 

Their analysis [230J uses method 2 with the CTEQ6.6 free-nucleon PDFs as their basis set. It is 
important to note that instead of using the full covariant error matrix, they have chosen to add the 
statistical and systematic uncertainties in quadrature when computing \ 2 - This implies that the total 
errors they are using are larger than those used by nCTEQ and thus have less discriminatory power 
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when compared to other results. They form a ratio of data versus theory that is based on CTEQ6.6: 

pCTEQe.6 = (Experimental) 

K ~ a".* (CTEQ6.6) ' [i6b) 

When using vA as the experimental numerator and comparing it to a denominator calculated with 
free-nucleon PDFs (CTEQ6.6) a rather poor fit could be naturally expected. Then they form a second 
purely theoretical ratio: 

OCTEQ6.6XEPS09 = ^ (CTEQ6.6 X EPS09) 

^ (CTEQ6.6) ' { } 

using the nuclear correction factors from their EPS09 [232] analysis and introducing CTEQ6.6 PDF 
uncertainties based in turn on upper and lower cross section uncertainties based on the PDF uncertainty 
analysis. If the nuclear correction factors determined by EPS09, based on charged-lepton nucleus 
scattering, are correct also for vA scattering, then i^ CTE Q 6 - 6 should be consistent with ^ CTE Q 6 - 6xEPS09 . 

They report their conclusions as a function of additional necessary radiative and target mass cor- 
rections that are made either separately, together or not at all. An example, from Table 1 of their 
publication, is their findings for the \ 2 /iV- values when making neither and then both the radiative and 
target-mass corrections: 



Radiative and Target Mass corrections 


CTEQ6.6 


CTEQ6.6xEPS09 


NuTeV 


1.51 


1.05 


CHORUS 


1.15 


0.79 


CDHSW 


1.10 


0.71 


No Radiative or Target Mass corrections 


CTEQ6.6 


CTEQ6.6xEPS09 


NuTeV 


1.35 


1.08 


CHORUS 


1.23 


1.09 


CDHSW 


0.96 


0.86 



Table 3: The x 2 /iV-values computed using CTEQ6.6 with and without nuclear modification from EPS09. The 
numbers are given for calculations with and without the radiative and the target mass corrections. 



From their comparison of the change in \ 2 as they add these corrections, they conclude: 

• Whatever way they make the calculation, the one without nuclear corrections from EPS09 gives 
substantially larger y 2 . 

• For CHORUS and CDHSW data, the x 2 s get consistently smaller as the radiative and target mass 
are applied. 

• For NuTeV data, the y 2 remains practically unchanged whether radiative or target mass correc- 
tions are applied or not. 



They note that in their fit to the NuTeV data shown in Figure [50j their predictions are in perfect 
agreement for some NuTeV energy ranges, whereas there are significant and inconsistent deviations in 
other energy bins. 



Eskola-Paukkunen-Salgado: Shadowing in v-A Scattering Their overall conclusion of the anal- 
ysis of CDHSW, CHORUS and NuTeV data sets is that only the NuTeV neutrino data is inconsistent 
with nuclear PDFs determined from charged-lepton nucleus scattering. They attribute the NuTeV 
discrepancy to an unexplained, neutrino energy-dependent fluctuations in the data and assert that the 
NuTeV data should not be used in making conclusions about the universality of nuclear PDFs. Without 
the NuTeV data they conclude that there is no inconsistency between v-A and £-A scattering in any of 
the x-dependent nuclear effects regions - including the shadowing region. 
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Figure 50: The Q 2 -averaged NuTeV neutrino data compared to the EPS09*CTEQ6.6 fit. 

8.3.4 Global Fits: The Hirai-Kumano-Nagai Nuclear Parton Distributions and Correc- 
tion Factors 

This group extracts nPDFs using method 2, by fitting weight factors. They have a complete set of nPDFs 
for fits to data other than v-A scattering. They do comment on nuclear effects in v-A scattering [233J , 
but do not make explicit fits including v-A scattering results. 

8.3.5 Global Fits: The de Florian-Sassot (dFS) Nuclear Parton Distributions and Cor- 
rection Factors 

The de Florian-Sassot global analyses yielding nPDFs began in 2004 [220] . They were the first group 
to apply a NLO analysis to nuclear data. In their initial analyses they used a convolution relation to 
relate the "free-nucleon" PDFs to the bound nPDFs via a weight-function Wi(y, A, Z): 



tf(x N , Ql) = f ^Wi{y, A Z) f^, Ql) (138) 

Jx N y y 



that allowed them to perform evolutions in Mellin space for increased accuracy. They used the 
GRV98 [234] free- proton PDFs as their basis set and fit the weight functions to a variety of charged- 
lepton DIS and Drell-Yan large-A data. They performed no error analysis of the nPDF sets. Their 



extracted values of F^jF® = R^(fi) are shown in Figure 51. As can be seen they fit quite well the 
shadowing region for all values of A. 

The group latest global analysis |225] now also includes P. Zurita and M. Stratmann as co-authirs. 
There are considerable updates from the first version outlined above. They have replaced the proce- 
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Figure 51: F^/F® data. The lines interpolate the values obtained with the NLO nPDF set at the 
respective Q 2 , and extrapolate to low xjy 



dure based on the convolution integral Eq.(138) with the more conventional approach of a "nuclear 



modification factor" (nuclear correction factor) given by Eq.(135). They also have chosen to use the 



MSTW |235] NLO free-proton PDFs as their basis set. Their strategy is to parametrize the Rf(xN, Qc 



in Eq. (134) for individual quarks and gluons as, for example, the valence quarks, where both valence 
quark distributions are assigned the same nuclear modification factor R^), which they parametrize as: 



R${x,Ql) = e 1 x a "(l-x) /?1 x 

(l + e 2 (l-x) /J2 )(l + a„(l-x) 



(139) 



They have expanded the data set used in their fit to include not only DIS and Drell-Yan off larger-A 
targets, but now also u-A scattering and inclusive pion production in deuteron-gold collisions. 

Concentrating on their fit to u-A results, they include the data sets from NuTeV (Fe), Chorus (Pb) 
and CDHSW (Fe). They fit to the averaged structure functions i*2,3 = (F%f + -F 2 ^)/2 from appropriate 
linear combinations of neutrino and antineutrino CC DIS differential cross sections, as opposed to the 
differential cross section for v and v independently. Their treatment of data point uncertainties uses 
the statistical and systematic uncertainties of the data added in quadrature as compared to the nCTEQ 
analysis that uses the full correlated covariant error matrix provided by the NuTeV collaboration. Both 
of these choices make considerable difference when compared to the nCTEQ analysis. They conclude 
that the CC data for the averaged structure function F2 are quite well reproduced both in shape and 



in magnitude, within the experimental uncertainties, by their fit, as illustrated by Fig. 52 
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Figure 52: The de Florian et al.NLO fit compared to the averaged CC DIS structure function F 2 as 
a function of Q 2 in various bins of xn for NuTeV (Fe) and CHORUS (Pb) targets shown as solid and 
dashed lines, respectively. 

de Florian-Sassot: Shadowing in u-A Scattering On the basis of their analysis they conclude 
that there is no disagreement between neutrino and charged-lepton nuclear correction factors and that 
a common set of nPDFs fitting both neutrino and charged lepton data over the full x range, including 
the shadowing region, can be found. 

8.4 Comparison of the t^A and vA Nuclear Correction Factors 

There is an obvious difference in the conclusions drawn by the nCTEQ group compared to those drawn 
by the EPS and dFS groups with respect to the I A and vA nuclear correction factors, particularly for 
the shadowing region. 

For the nCTEQ analysis, the contrast between the charged-lepton {£ A) case and the neutrino 
(uA) case is striking. While their charged-lepton results that fit to charged-lepton and Drell-Yan 
data generally align with the SLAC/NMC and HKN determinations, the neutrino results clearly yield 
different behavior as a function of x, particularly in the shadowing/anti-shadowing region. In the V 
case, these differences are smaller but persist in the low-x shadowing region. They emphasize that both 
the charged-lepton and neutrino results come directly from global fits to the data, and there is no model 
involved. They further suggest that this difference between the results in charged-lepton and neutrino 
DIS is reflective of the long-standing "tension" between the light-target charged-lepton data and the 
heavy-target neutrino data in the historical global PDF fits [2361 1237] . Their latest results further 
suggest that the tension is not only between charged-lepton light-target data and neutrino heavy-target 
data, but also between neutrino and charged-lepton heavy-target data. In other words, a difference 
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between charged-lepton (^A) and the neutrino (uA) when comparing the same A. 

For both the EPS and dFS no significant difference is found between the charged-lepton (^A) and 
the neutrino {vA) cases, and particularly no difference in the shadowing region, where the nCTEQ 
group finds a significant difference and where differences are actually expected theoretically, although 
certainly not of the magnitude as found by the nCTEQ group. 

Concentrating on this interesting difference found by the nCTEQ group, if the nuclear corrections for 
the l^A and vA processes are indeed different there are several far-reaching consequences. Considering 
this, the nCTEQ group has performed a unified global analysis |221j of the £^A, DY, and vA data 
(accounting for appropriate systematic and statistical errors) to determine if it is possible to obtain a 
"compromise" solution including both UrA and vA data. Using a hypothesis-testing criterion based on 
the x 2 distribution that can be applied to both the total x 2 as well as to the x 2 °f individual data sets, 
they found that it was not possible to accommodate the data from vA and t A DIS by an acceptable 
combined fit. That is, when investigating the results in detail, the tension between the £ Fe and vFe 
data sets does not permit a compromise fit, which adequately describes the neutrino DIS data along with 
the charged-lepton data and, consequently, £^Fe and vFe based on the NuTeV results, have different 
nuclear correction factors. 

A compromise solution between vA and £ A data can be found only if the full correlated system- 
atic errors of the vA data are not used and the statistical and all systematic errors are combined in 
quadrature, thereby neglecting the information contained in the correlation matrix. In other words, 
the larger errors resulting from combining statistical and all systematic errors in quadrature reduces 
the discriminatory power of the fit, such that the difference between vA and l^A data are no longer 
evident. This conclusion underscores the fundamental difference |221] of the nCTEQ analysis with the 
analyses of EPS and dFS. 

On the other hand, a difference between vA and i^A is not completely unexpected, particularly 
in the shadowing region, and has previously been discussed in the literature |238[ I78| 12381 1239] . The 
charged-lepton processes occur (dominantly) via 7-exchange, while the neutrino-nucleon processes occur 
via H^-exchange. The different nuclear corrections could simply be a consequence of the differing 
propagation of the hadronic fluctuations of the intermediate bosons (photon, W) through dense nuclear 
matter. Furthermore, as stressed above, since the structure functions in neutrino DIS and charged 
lepton DIS are distinct observables with different parton model expressions, it is clear that the nuclear 
correction factors will not be exactly the same. What is, however, unexpected is the degree to which 
the R factors differ between the structure functions F% Fe and F| Fe . In particular the lack of evidence 
for shadowing in neutrino scattering down to x ~ 0.02 at low Q 2 is quite surprising. 



8.5 Evidence for Shadowing in v-A Scattering 

The ratio of F^(x, Q 2 ) with respect to the free-nucleon F 2 (x, Q 2 ) naturally reflects the general behavior 
of the cross section. However, instead of integrating over Q 2 to get a cross section in a given x interval 
as was done above, the x-behavior of the structure function within a specific Q 2 range can be examined. 



For the lower Q 2 regions of the average F 2 Fe (x,Q 2 ) shown already in Figure 49, the comparison 
between the nCTEQ u-A fit and the nCTEQ charged-lepton fit suggests very different behavior. There is 
no indication of shadowing in the lowest-x, lowest-Q regions, with the z/-A fit only minimally approaching 
the charged-lepton fit with increasing Q. However, even when examining the higher Q 2 range, as shown 



in Figure 53, where the low-x shadowing region becomes kinematically inaccessible as Q increases, the 
slope of the v-K fit approaching the shadowing region from higher x for 1.2 < Q 2 < 30 GeV 2 would 
make it difficult to mimic the shadowing evidenced in charged-lepton nucleus scattering. 

It is important to emphasize that this apparent disagreement between z/A and £ A shadowing 
behavior is based only on the NuTeV data and the use of their full covariant error matrix in the global 
analysis. This disagreement was not seen by the two other analyses that, however, did not use the full 
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Figure 53: Nuclear correction factor R for the average F2 structure function in charged current vFe 
scattering at Q 2 = 8.0, 12.5, 20.0 and 31.6 GeV 2 compared to the measured NuTeV points. The green 
dashed curve shows the result of the nCTEQ analysis of v A differential cross sections plotted in terms 
of the average _F 2 Fe divided by the results obtained with the reference fit (free-proton) PDFs. For 
comparison, the nCTEQ fit to the charged-lepton data from Figure 46 is shown by the solid blue curve. 



discriminatory power of the NuTeV data. Figure 54 summarizes the results of various global fits to both 
vA and l^A scattering. The curves marked nCTEQ £^A, EPS09 and HKN07 are fits to I A scattering, 
while the nCTEQ vA and DSSZ preliminary are fits to vA scattering. The DSSZ preliminary fit is 
made to the average F2 and does not use the full NuTeV covariant error matrix, while the nCTEQ 
fit is to the double-differential cross sections and does use the full covariant cross section. Before any 
far-reaching conclusions can be drawn, it is important to gather significantly more data from other 
experiments using other nuclei. 



9 Current Experimental Studies of u-A Nuclear Effects 

9.1 The MINERz/A Experiment at the FNAL Main Injector Neutrino 
Beam (NuMI) 

One such experiment that is currently taking data is the MINERVA (Main Injector ExpeRiment: v A) 
experiment [240] . a collaboration of elementary-particle and nuclear physicists. The experiment has 
installed a fully active fine-grained solid scintillator detector in the NeUtrinos from the Main Injector 
(NuMI) beam. The overall goals of the experiment are to measure absolute exclusive cross-sections, 
study nuclear effects in v - A interactions (with A varying from He to Pb), and perform a systematic 
study of the resonance-DIS transition and the full DIS region. A determination of these x-dependent 
nuclear effects with neutrino and antineutrino off a wide range of nuclear targets will start with the 
current exposure and will continue when the neutrino beam is tuned to higher energies starting in 2013. 



A schematic of the MINERz/A detector is shown in Fig. 55 The detector consists of five main 



regions: the fully active central detector, the upstream nuclear targets, a downstream electromagnetic 



65 



1.15 
1.10 
1.05 
1.00 
0.95 






(f- 


= 5GeV 






























































































































0.90 










































0.85 














n 


CTE 






















0.80 
0.75 














n 

b 


CTI 
PS0 c 


3 




I'A 


























r 

H 


SS/ 
KNO 


P'C 

7 


im 


















0. 


01 




0.1 1 




a) 

Figure 54: Nuclear correction factor R for the average F2 structure function in vA and £ A scattering 
for a) Q 2 = 5 and b) 20.0 GeV 2 compared to the measured NuTeV average F2 points shown with 
statistical errors. The blue curve shows the result of the nCTEQ analysis of v A (Fe and Pb) differential 
cross sections, plotted in terms of the average F 2 Fe divided by the results obtained with the reference 
fit (free-proton) PDFs. The green curve is the PRELIMINARY deFlorian et al fit to the NuTeV 
average structure functions. For comparison, the nCTEQ fit, the EPS09 fit and the HKN07 fit to the 
charged-lepton data are also shown, to emphasize the difference in the shadowing region 



and hadron calorimeter, and a surrounding electromagnetic and hadron calorimeter. 

The central fully-active detector serves as both the primary target and the tracking detector. The 
downstream electromagnetic calorimeter consists of alternating planes of Pb and scintillator planes, 
while the hadron calorimeter is similar, with planes of steel instead of Pb. The side electromagnetic 
calorimeter consists of Pb plates between tracking plane in the outer region of the central detector. The 
side hadron calorimeter consists of planes of steel with scintillator strips embedded. 

Upstream of the central detector are planes of passive targets, with two planes 2.5 cm thick of mixed 
Fe/Pb, one plane with 2.5 cm thick Fe/Pb and 7.5 cm C, a solid plane of Pb 0.80 cm thick, and a mixed 
plane of Fe/Pb 1.30 cm thick. The mixed Fe/Pb planes are split with part of the plane being iron 
and part of the plane being lead, such that the total mass is approximately equal. Tracking planes 
are placed between each plane of passive targets. Recently a 15 cm water target has been installed 
downstream of the Fe/Pb/C target and a tank of liquid 4He about 1 m in diameter has been installed 
upstream of the main detector. Figure 56 shows a schematic of the solid nuclear target region along 
with the placement of carbon, lead, and iron in the various targets. 

Charged current events originating in the central detector are fully contained, except for the muon, 
for neutrino energies of less than about 10 GeV. The MINOS near detector, directly downstream of 
MINERz/A, acts as a muon spectrometer and gives both muon energy and charge for forward going 
muons. For particles stopping in MINERz/A, particle identification can be determined from the dE/dX, 
but there is no charge determination. 

The NuMI beam is currently being converted to the Medium Energy (ME) configuration for the 
NOvA experiment. Figure 57 displays the expected produced (no corrections for acceptance or recon- 
struction) statistics in the fiducial volume of the fully-active scintillator tracker for 6 x 10 20 protons on 
target. The units of the z-axis are lk events and this represents the statistics we can expect per 1 year 
of running, when the Medium Energy beam configuration is running at expected power. The statistics 



in the Fe and Pb targets will be about 1/3 the statistics displayed in Figure 57 



Since MINERz/A plans to be running for many years in the ME configuration with NOvA, it will be 
able to determine the nuclear dependence of the inclusive cross section, and in particular the relative 
ratio of cross sections off lead, iron, water and carbon to a statistical precision of better than 1%. 
MINERvA has also proposed to fill the present He cryogenic target with deuterium providing a direct 
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Figure 55: Schematic side view of the MINERVA detector. 




ratio of A/D2. As can be seen, even with the ME beam, MINERVA will not be able to reach deep 
into the shadowing region. However, several data points integrated over Q 2 between x = 0.04 and 0.10 
should be obtainable with reasonable statistics. This should be sufficient to determine the direction of 
the ratio of nuclear to nucleon cross sections as it enters the shadowing region. 



10 Conclusions: Experimental Evidence for Shadowing in Elec- 
troweak Interactions 

Shadowing with t A scattering has been studied for decades and its experimental evidence is now 
overwhelming and conclusive. The A-dependence of the shadowing phenomena over a wide range of 
nuclei has been measured with both electron and muon beams, and has been fitted to various theoretical 
models, as summarized in the theoretical part of this review. 

What is new in lepton nucleus scattering is the introduction of statistically significant analyses of vA 
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Figure 57: Contour of the x vs Q 2 distribution of events produced (not corrected for acceptance or 
reconstruction efficiencies) in the ME beam configuration. The z-axis units are 1000 events and are for 
a 3-ton scintillator fiducial volume and 6 x 10 20 protons on target. 



scattering. The experimental evidence for shadowing in v A scattering is a comparatively new study 
and is therefore much less settled. In this review of experimental data the emphasis has thus been 
placed on these new vA scattering results and two very different and distinct conclusions have been 
reached among groups studying them. Whereas all three analyses - nCTEQ, EPS and dFS - agree with 
both the shape and the magnitude of the shadowing exhibited by l^A scattering, there is disagreement 
between the nCTEQ analysis and those of EPS and dFS when considering vA scattering. 

The nCTEQ analysis suggests an apparent difference in z/-A and £ A shadowing behavior, which is 
not seen by the analyses of EPS and dFS. However EPS and dFS did not use the full NuTeV covariant 
error matrix, yielding the full discriminatory power of the NuTeV data in their analyses, while nCTEQ 
did use the full covariant error matrix. 

It is important to again emphasize that any significant conclusions on shadowing with v and V 
nucleus scattering that we are able to draw here are based on ONE analysis of ONE experiment using 
ONE nucleus, the NuTeV v Fe scattering data. The other current data set used in these analyses, 
CHORUS v Pb scattering, has significantly larger systematic errors than NuTeV and thus plays a much 
reduced role in the analysis. Before any far-reaching conclusions can be drawn, it is important to gather 
significantly more data from other experiments, using other nuclei. It is expected that the MINERz/A 
experiment, soon to be taking data in the higher energy ME beam configuration at Fermilab, will 
provide data to compare with the NuTeV results. 

Should subsequent experimental results confirm the difference between charged-lepton and neutrino 
scattering in the shadowing region at low-Q 2 , it is interesting to speculate on the possible cause of 
the difference. A recent study of EMC, BCDMS and NMC data by a Hampton University - Jefferson 
Laboratory collaboration |241] suggests that anti-shadowing in charged-lepton nucleus scattering may 
be dominated by the longitudinal structure function F L . As a by-product of this study, their figures 
hint that shadowing in the data of EMC, BCDMS and NMC fi A scattering was being led by the 
transverse cross section with the longitudinal component crossing over into the shadowing region at 
lower x compared to the transverse. This is consistent with the expectations expressed in the equations 



and figures of 3.2.3 



In the low-Q 2 region, the neutrino cross section is dominated by the longitudinal structure function 
Fl via axial-current interactions, since Ft vanishes as Q 2 as Q 2 =^ 0, similar to the behavior of charged 
lepton scattering. If the results of the NuTeV analysis are verified, one contribution to the different 
behavior of shadowing at low-Q 2 demonstrated by v A and i A, in addition to the different hadronic 
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fluctuations in the two interactions, could be due to the different mix of longitudinal and transverse 
contributions to the cross section of the two processes in this kinematic region. 

Another interesting hypothesis to explain this possible difference was recently proposed by V. 
Guzey [242J. He suggests that at low-x and low-Q 2 the value of y is close to unity and the neutrino 
interactions primarily probe the down and strange quarks. This is very different from the situation with 
charged-lepton scattering, where the contribution from down and strange quarks are suppressed by a 
factor of 1/4 compared to the up and charm. If the strange quark is not subject to the same degree of 
shadowing as the valence quarks, then the neutrino cross section in this kinematic regime would also not 
exhibit strong shadowing. On the other hand, the charged-lepton cross section, with its comparatively 
small strange quark content in this same kinematic regime, would still exhibit strong shadowing. 

Acknowledgments 

This work was supported in part by Fondecyt (Chile) grants 1090291 and 1100287. The work of 
B.Z.K. was supported also by the Alliance Program of the Helmholtz Association, contract HA216/EMMI 
"Extremes of Density and Temperature: Cosmic Matter in the Laboratory". Fermilab is operated by 
Fermi Research Alliance, LLC under Contract No. De-AC02-07CH11359 with the United States De- 
partment of Energy. 



References 

R. J. Glauber, in Lectures in Theoretical Physics,Editois: W. E. Brittin et al. (New York, 1959). 

G. Antchev et. al. [TOTEM Collaboration ], CERN-PH-EP-2011-158 (Geneva,CERN,2011). 
P.V.R. Murthy et al., Nucl. Phys. B 92 (1975) 269 

V. N. Gribov, Sov. Phys. JETP 29 (1969) 483 [Zh. Eksp. Teor. Fiz. 56 (1969) 892]; Sov. Phys. 
JETP 30 (1970) 709 [Zh. Eksp. Teor. Fiz. 57 (1970) 1306]. 
V. Karmanov and L.A. Kondratyuk,Sov. Phys. JETP Lett. 18 (1973) 266 
B.Z. Kopeliovich and J. Nemchik, Phys. Lett. B 368 (1996) 187 

A. Gsponer et aX.,Phys. Rev. Lett. 42 (1979) 9 
J. Pumplin and M. Ross,Phys. Rev. Lett. 21 (1968) 1778 
E. Feinberg and I.Ya. Pomeranchuk,A^iiovo. Cimento. Suppl. 3 (1956) 652 
M.L. Good and W.D. Walker, Phys. Rev. 120 (1960) 1857 

B. Z. Kopeliovich and L.I. Lapidus,Sov. Phys. Sov. Phys. JETP Lett. 28 (1978) 664 

H. I. Miettinen and J. Pump\m,Phys. Rev. D 18 (1978) 1696 

B. Z. Kopeliovich, L. I. Lapidus and A. B. Zamolodchikov, JETP Lett. 33 (1981) 595 [Pisma Zh. 
Eksp. Teor. Fiz. 33 (1981) 612]. 

B. Z. Kopeliovich,P/i?/s. Rev. C 68 (2003) 044906. 

B.Z. Kopeliovich and B.G. Zakharov, Phys. Rev. D 44 (1991) 3466 

B.Z. Kopeliovich,J. Raufeisen and A.V. Tarasov, Phys. Lett. B 440 (1999) 151 

B.Z. Kopeliovich,J. Raufeisen and A.V. Tarasov, Phys. Rev. C 62 (2000) 035204 

B. Z. Kopeliovich,I. K. Potashnikova and I. Schmidt, Phys. Rev.C 73 (2006) 034901 

M. Alvioli,C. Ciofi degli Atti, B. Z. Kopeliovich, I. K. Potashnikova and I. Schmidt, Phys. Rev. C 
81 (2010) 025204. 

[20] C. Ciofi degli Atti,B. Z. Kopeliovich,C. B. Mezzetti,I. K. Potashnikova,I. Schmidt, Phys. Rev. C84 
(2011) 025205. 



69 



[21] N. Bianchi,E. De Sanctis,M. Mirazita and V. Muccifora, Phys. Rev. C 60 (1999) 064617 
[22] T. H. Bauer,F. Pipkin,R. Spital and D. R. Yennie,Rev. Mod. Phys. 50 (1978) 261 
[23] B. Z. Kopeliovich,I. Schmidt,M. Siddikov, Phys. Rev. D79 (2009) 034019. 

[24] B. Z. Kopeliovich,! Schmidt,M. Siddikov, Phys. Rev. D80 (2009) 054005. ibid Phys. Rev. D81 
(2010) 094013. 

[25] M. L. Goldberger and S. B. Treiman, Phys. Rev. 110,1178 (1958). 
[26] S. L. Adler, Phys. Rev. 135,B963 (1964). 

[27] M. Bell,K. Gottfried and M. J. G. Veltman, Singapore, Singapore: World Scientific (1995) 933 p. 
(World Scientific series in 20th century physics. 9) 

[28] C. A. Piketty and L. Stodolsky, Nucl. Phys. B 15,571 (1970). 

[29] B. Z. Kopeliovich and P. Marage, Int. J. Mod. Phys. A 8,1513 (1993). 

[30] J.S. Bell,Phys. Rev. Lett. 13 (1964) Phys. Rev. Lett. 13 (1964) 57 

[31] B. Z. Kopeliovich, Phys. Lett. B 227,461 (1989). 

[32] B. Z. Kopeliovich, Nucl. Phys. Proc. Suppl. 139,219 (2005). 

[33] P.P. Allport et al. [WA59 Collaboration], Phys. Lett. B 232 (1989) 417 

[34] A. A. Belkov and B. Z. Kopeliovich, Sov. J. Nucl. Phys. 46,499 (1987) [Yad. Fiz. 46,874 (1987)]. 
[35] B. Z. Kopeliovich,I. K. Potashnikova,I. Schmidt, J. Soffer, [arXiv: 1109.2500 [hep-ph]]. 
[36] R. T. Deck, Phys. Rev. Lett. 13,169 (1964). 

[37] K. Nakamura et al. (Particle Data Group), J. Phys. G 37,075021 (2010). 

[38] B. Z. Kopeliovich,I. K. Potashnikova,I. Schmidt,M. Siddikov, Phys. Rev. C84 (2011) 024608. 

[39] B. Z. Kopeliovich,I. Schmidt,M. Siddikov, Phys. Rev. D84 (2011) 033012. 

[40] B. Z. Kopeliovich,I. Schmidt,M. Siddikov, [arXiv: 1108.5654 [hep-ph]]. 

[41] A. Miicklich,Ph. D. thesis, University of Heidelberg (1995), in German. 

[42] M. Breidenbach et al.,Phys. Rev. Lett. 23 (1969) 935. 

[43] J. D. Bjorken,Phys. Rev. 179 (1969) 1547. 

[44] R. P. Feynman,Phys. Rev. Lett. 23 (1969) 1415. 

[45] B. Povh,K. Rith,Ch. Scholz,and F. Zetsche, Teilchen und 

iTerne,Springer,Heidelberg,Germany,1994. 

[46] C. G. Callan and D. J. Gross, Phys. Rev. Lett. 22 (1969) 156. 

[47] V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438. 

[48] L. N. Lipatov,Sov. J. Nucl. Phys. 20 (1975) 95. 

[49] G. Altarelli and G. Parisi, Nucl. Phys. B126 (1977) 298. 

[50] Yu. L. Dokshitzer,Sov. Phys. JETP 46 (1977) 641. 

[51] R. K. Ellis, W. J. Stirling,and B. R. Webber, QCD and Collider Physics, Cambridge University 
Press, Cambridge, UK, 1996. 

[52] M. Gliick,E. Reya,and A. Vogt,Z. Phys. C67 (1995) 433. 

[53] A. D. Martin,R. G. Roberts,W. J. Stirling,and R. S. Thorne,Eur. Phys. J. C4 (1998) 463. 

[54] J. Pumplin,D. R. Stump,J. Huston,H. L. Lai,P. M. Nadolsky and W. K. Tung, JHEP 0207 (2002) 
012 

[55] J. C. Collins,D. E. Soper,G. F. Sterman, Adv. Ser. Direct. High Energy Phys. 5,1-91 (1988). 
[56] O. V. Kancheli,Pisma Zh. Eksp. Teor. Fiz. 18 (1973) 465. 
[57] B. Kopeliovich,B. Povh, Phys. Lett. B367 (1996) 329 



70 



58] B. Kopeliovich,B. Povh, Z. Phys. A356 (1997) 467 

59] CERN NMC,P. Amaudruz et al.,Z. Phys. C 51 (1991) 387 

60] CERN NMC,M. Arneodo et al.,Nucl. Phys. B 441 (1995) 278 

61] K. Prytz, Phys. Lett. B 311,286 (1993); Phys. Lett. B 332,393 (1994). 

62] T. Gousset and H. J. Pirner, Phys. Lett. B 375,349 (1996) [hep-ph/9601242]. 

63] G. Ingelman,P. E. Schlein,Phys. Lett. B152 (1985) 256. 

64] Y. .M. Kazarinov,B. Z. Kopeliovich,L. I. Lapidus,I. K. Potashnikova, Sov. Phys. JETP 43 (1976) 
598. 

65] V. A. AbramovskyV. N. Gribov,0. V. Kancheli, Yad. Fiz. 18 (1973) 595 

66] Frankfurt L,Guzey V and Strikman M 2005 Phys. Rev. D 71 054001 

67] L. Frankfurt,V. GuzeyM. Strikman, [arXiv:1106.2091 [hep-ph]]. 

68] N. Armesto,A. B. Kaidalov,C. A. Salgado,K. Tywoniuk, Eur. Phys. J. C68 (2010) 447 

69] A. Schwimmer,Nucl. Phys. B 94 (1975) 445. 

70] Balitsky I 1996 Nucl. Phys. B 463 99 

71] Kovchegov Yu V, Phys. Rev. D 60 034008 

72] Melnitchouk W and Thomas A W Phys. Rev. D 47 (1993) 3783. 

73] Pillcr G and Weise W 2000 Phys. Rept. 330 1 

74] Kwiecinski J and Badelek B 1988 Phys. Lett. B 208 508 

75] Shaw G 1993 Phys. Rev. D 47 R3676 

76] Brodsky S J and Lu H j 1990 Phys. Rev. Lett. 64 1342 

77] Capella A,Kaidalov A, Merino C,Pertermann D and Tran Thanh Van J 1998 Eur. Phys. J. C 5 111 

78] Kulagin S A and Petti R 2006 Nucl. Phys. A 765 126 

79] Ioffe B.L. Khoze V.A. Lipatov L. N. Hard processes 1. Phenomenology Quark-Parton Model (El- 
sevier, 1984) (ISBN 0444862137) 

80] B. K. Jennings and B. Z. Kopeliovich, Phys. Rev. Lett. 70 (1993) 3384. 

81] J. D. Bjorken, In *Lake Louise 1995,Quarks and colliders* 61-90, SLAC-PUB-95-6949; 
http://www.slac.stanford.edu/pubs/slacpubs/6750/slac-pub-6949.pdf 

82] A. H. Mueller and J. -w. Qiu, Nucl. Phys. B 268 (1986) 427. 

83] J. -w. Qiu, Nucl. Phys. B 291 (1987) 746. 

84] N. N. Nikolaev and B. G. Zakharov,Z. Phys. C49 (1991) 607. 

85] J. B. Kogut and D. E. Soper, Phys. Rev. Dl (1970) 2901. 

86] J. M. Bjorken,J. B. Kogut,and D. E. Soper, Phys. Rev. D3 (1971) 1382. 

87] K. Golec-Biernat and M. Wusthoff,Phys. Rev. D59 (1999) 014017; Phys. Rev. D60 (1999) 114023. 

88] B. Z. Kopeliovich,A. Schafer and A. V. Tarasov, Phys. Rev. D 62 (2000) 054022. 

89] J. Bartels,K. J. Golec-Biernat and H. Kowalski, Phys. Rev. D 66 (2002) 014001 

90] S.J. Brodsky and A. Mueller, Phys. Lett. B206 (1988) 685. 

91] G. Bertsch,S.J. BrodskyA.S. Goldhaber and J.F. Gunion,Phys. Rev. Lett. 47 (1981) 297. 

92] J. M. Bjorken,J. B. Kogut,Phys. Rev. D8 (1973) 1341. 

93] L.L. Frankfurt and M.I. Strikman, Phys. Rept. 160 (1988) 235. 

94] R. Feynman and A.R. Hibbs, Quantum mechanics and path integrals (McGraw-Hill, New York,1965 
)• 

[95] NM Collab.,M. Arneodo et al, Nucl. Phys. B481 (1996) 23. 



71 



96] NM Collab.,P. Amaudruz et al. Nucl. Phys. B441 (1995) 3. 
97] E665 Collab.,M. R. Adams et al. Phys. Rev. Lett. 68 (1992) 3266. 
98] E665 Collab.,M. R. Adams et al. Z. Phys. C67 (1995) 403. 
99] J. Raufeisen, Ph. D. thesis, University of Heidelberg, (2000) 

100] A. Levy, Low-x Physics at HERA in Lectures on Q CD, Applications, F. Lenz, H. Griecrhammer and 
D. Stoll (eds.), Springer, Heidelberg, Germany,1997. 

101] NM collab.,P. Amaudruz et al. Z. Phys. C53 (1992) 73. 

102] L. W. Mo and Y. S. Tsai,Rev. Mod. Phys. 41 (1968) 205. 

103] A. A. Akhundov,D. Y. Bardin,and N. M. Shumeiko, Sov. J. Nucl. Phys. 26 (1977) 660. 
104] D. Y. Bardin and N. M. Shumeiko, Sov. J. Nucl. Phys. 29 (1979) 499. 
105] A. A. Akhundov,D. Y. Bardin,and N. M. Shumeiko, Sov. J. Nucl. Phys. 44 (1986) 988. 
106] J. J. van Hunen,Ph. D. thesis Utrecht 2000. 

107] The NM Coll.,M. Arneodo et al., Nucl. Phys. B 487 (1997) 3 [hep-ex/9611022]. 

108] K. Ackerstaff et al. [HERMES Collaboration], Phys. Lett. B 475 (2000) 386 

109] K. Ackerstaff et al. [HERMES Collaboration], Erratum: Phys. Lett. B 567 (2003) 339 

110] K. J. Eskola,V. J. Kolhinen and P. V. Ruuskanen, Nucl. Phys. B 535 (1998) 351; K. J. Es- 
kola,V. J. Kolhinen and C. A. Salgado, Eur. Phys. J. C 9 (1999) 61. 

Ill] B. Z. Kopeliovich,J. Raufeisen, A. V. Tarasov,M. B. Johnson, Phys. Rev. C67 (2003) 014903. 

112] S. J. BrodskyP. Hoyer,N. Marchal,S. Peigne and F. Sannino, Phys. Rev. D 65 (2002) 114025 

113] M. Arneodo,Phys. Rept. 240 (1994) 301. 

114] K. Saito,K. Tsushima,A. W. Thomas, Prog. Part. Nucl. Phys. 58 (2007) 1 
115] V. Barone et al., Z. Phys. C 58 (1993) 541. 

116] B. Z. Kopeliovich,I. K. Potashnikova,B. Povh,I. Schmidt, Phys. Rev. D76 (2007) 094020. 
117] T. Schafer and E.V. Shuryak,Rev. Mod. Phys. 70 (1998) 323. 
118] E. Shuryak and I. Zahed,Phys. Rev. D69 (2004) 014011. 

119] J. Hufner,Y. .P. Ivanov,B. Z. Kopeliovich,A. V. Tarasov, Phys. Rev. D62 (2000) 094022; ibid 
C66 (2002) 024903. 

120] B. Z. Kopeliovich,A. V. Tarasov, Nucl. Phys. A710 (2002) 180. 

121] V. Barone and M. Genovese,hep-ph/9610206. 

122] H. Fritzsch and P. Minkowski, Phys. Lett. 73B (1978) 80. 

123] G. Altarelli,G. Parisi,and R. Petronzio, Phys. Lett. 76B (1978) 351. 

124] G. Altarelli,G. Parisi,and R. Petronzio, Phys. Lett. 76B (1978) 356. 

125] P. L. McGaugheyJ. M. Moss,and J. C. Peng, Ann. Rev. Nucl. Part. Sci. 49 (1999) 217. 

126] G. Altarelli,R. K. Ellis,M. Greco,and G. Martinelli, Nucl. Phys. B246 (1984) 12. 

127] J. C. Collins,D. E. Soper,and G. Sterman, Nucl. Phys. B250 (1985) 199. 

128] P. Chiappetta and H. J. Pirner, Nucl. Phys. B291 (1987) 765. 

129] B. Z. Kopeliovich, proc. of the workshop Hirschegg '95: Dynamical Properties of Hadrons in 
Nuclear Matter, Hirschegg January 16-21, 1995, ed. by H. Feldmeyer and W. Norenberg, Darm- 
stadt, 1995,p. 102 (hep-ph/9609385). 

130] E772 collab.,P. L. McGaughey et al., Phys. Rev. D50(1994)3038; erratum ibid. D60(1999)119903. 

131] B. Z. Kopeliovich, A. Schafer, and A. V. Tarasov,Phys. Rev. C59 (1999) 1609, extended version 
in hep-ph/9808378. 



72 



[132 
[133 
[134 
[135 
[136 
[137 
[138 
[139 
[140 
[141 
[142 

[143 
[144 
[145 
[146 
[147 
[148 
[149 
[150 

[151 

[152 

[153 
[154 
[155 
[156 
[157 
[158 
[159 
[160 
[161 
[162 
[163 
[164 
[165 
[166 
[167 
[168 
[169 



S. J. BrodskyA. Hebecker,and E. Quack, Phys. Rev. D5519972584. 

B. Z. Kopeliovich,J. Raufeisen,A. V. Tarasov, Phys. Lett. B503 (2001) 91 

M. A. Betemps,M. B. G. Ducati,M. V. T. Machado,J. Raufeisen, Phys. Rev. D67, 114008 (2003). 

J. Raufeisen, J. -C. Peng and G. C. Nayak, Phys. Rev. D 66 (2002) 034024 [hep-ph/0204095]. 

M. B. Johnson,B. Z. Kopeliovich,A. V. Tarasov, Phys. Rev. C63 (2001) 035203. 

B. Z. Kopeliovich,J. Nemchik,A. Schafer,A. V. Tarasov, Phys. Rev. Lett. 88 (2002) 232303. 

B. Blattel,G. Baym,L. L. Frankfurt and M. Strikman, Phys. Rev. Lett. 70 (1993) 896. 

M. B. Johnson et al. [ FNAL E772 Collaboration ], Phys. Rev. Lett. 86 (2001) 4483. 

M. B. Johnson, et al, Phys. Rev. C65 (2002) 025203. 

D.M. Aide et al.,Phys. Rev. Lett. 64 (1990) 2479. 

B. Z. Kopeliovich,J. Nemchik,I. K. Potashnikova,M. B. Johnson, I. Schmidt, Phys. Rev. 
C72,054606 (2005). 
B. Z. Kopeliovich,J. Nemchik, J. Phys. G G38 (2011) 043101. 
M. A. Vasiliev,et al.,Phys. Rev. Lett. 83 (1999) 2304. 
M. Leitch,private communication. 

M. B. Johnson,B. Z. Kopeliovich,I. Schmidt, Phys. Rev. C75 (2007) 064905. 
B. Z. Kopeliovich,I. K. Potashnikova,I. Schmidt, Phys. Rev. C81 (2010) 035204. 
R. Brock et al. [CTEQ Collaboration], Rev. Mod. Phys. 67 (1995) 157. 
D. J. Gross and C. H. Llewellyn Smith, Nucl. Phys. B 14 (1969) 337. 

W. C. Leung,P. Z. Quintas,S. R. Mishra,F. J. Sciulli,C. Arroyo, K. T. Bachmann,R. E. Blair and 
C. Foudas et al, Phys. Lett. B 317 (1993) 655. 

J. H. Kim,D. A. Harris, C. G. Arroyo, L. de Barbara, P. de Barbara, A. O. Bazarko,R. H. Bernstein 
and A. Bodek et al, Phys. Rev. Lett. 81 (1998) 3595 

R. Fiore and V. R. Zoller, JETP Lett. 82 (2005) 385 [Pisma Zh. Eksp. Teor. Fiz. 82 (2005) 
440]; Phys. Lett. B 632 (2006) 87 

M. B. G. Ducati,M. M. Machado and M. V. T. Machado, Phys. Lett. B 644 (2007) 340 

J. -W. Qiu and I. Vitev, Phys. Lett. B 587 (2004) 52 

A. Airapetian et al. [HERMES Collaboration], hep-ex/0210068. 

J. J. Aubert et al. [European Muon Collaboration], Phys. Lett. B 123 (1983) 275. 

A. Bodek et al, Phys. Rev. Lett. 50 (1983) 1431. 

A. Bodek et al, Phys. Rev. Lett. 51 (1983) 534. 

R. G. Arnold et al, Phys. Rev. Lett. 52 (1984) 727. 

S. Dasu et al, Phys. Rev. Lett. 60 (1988) 2591. 

G. Bari et al. [BCDMS Collaboration], Phys. Lett. B 163 (1985) 282. 
A. C. Benvenuti et al. [BCDMS Collaboration], Phys. Lett. B 189 (1987) 483. 
J. Ashman et al. [European Muon Collaboration], Phys. Lett. B 202 (1988) 603. 
M. Arneodo et al. [European Muon Collaboration], Nucl. Phys. B 333 (1990) 1. 
P. Amaudruz et al. [New Muon Collaboration], Z. Phys. C 51 (1991) 387. 
P. Amaudruz et al. [New Muon Collaboration], Phys. Lett. B 294 (1992) 120. 
M. R. Adams et al. [Fermilab E665 Collaboration], Phys. Lett. B 287 (1992) 375. 
M. R. Adams et al. [E665 Collaboration], Phys. Rev. Lett. 68 (1992) 3266. 
M. R. Adams et al. [E665 Collaboration], Z. Phys. C 67 (1995) 403. 



73 



[170 
[171 
[172 

[173 

[174 

[175 
[176 
[177 

[178 
[179 
[180 

[181 
[182 
[183 



[184 
[185 
[186 
[187 
[188 

[189 

[190 

[191 



[192 
[193 

[194 
[195 
[196 
[197 
[198 
[199 



J. Gomez et al, Phys. Rev. D 49 (1994) 4348. 

D. F. Geesaman,K. Saito,A. W. Thomas, Ann. Rev. Nucl. Part. Sci. 45,337-390 (1995). 

J. Arrington,R. Ent,C. E. Keppel,J. Mammei and I. Niculescu, Phys. Rev. C 73,035205 (2006) 
[nucl-ex/0307012]. 

S. Dasu,P. deBarbaro,A. Bodek,H. Harada,M. W. KrasnyK. Lang,E. M. Riordan,L. Andivahis et 
al, Phys. Rev. D49,5641-5670 (1994). 

J. F. Owens, J. Huston, C. E. Keppel,S. Kuhlmann,J. G. Morfin,F. 01ness,J. Pumplin,D. Stump, 
Phys. Rev. D75,054030 (2007). [hep-ph/0702159 [HEP-PH]]. 

W. Melnitchouk,A. W. Thomas, Phys. Rev. C52,3373-3377 (1995). [hep-ph/9508311]. 

M. Arneodo et al. [New Muon Collaboration], Nucl. Phys. B 481,3 (1996). 

J. F. Owens, J. Huston, J. Pumplin,D. Stump, C. E. Keppel,S. Kuhlmann,J. G. Morfin,F. Olness, 
AIP Conf. Proc. 967,259-263 (2007). 

M. A. Parker et al. [ BEBC TST Neutrino Collaboration ], Nucl. Phys. B232,l (1984). 

A. M. Cooper-Sarkar et al. [ WA25 and WA59 Collaborations ], Phys. Lett. B141,133 (1984). 

J. Hanlon,R. A. Burnstein,H. A. Rubin, M. Kalelkar,E. B. Brucker,P. F. Jacques, E. L. Koller,R. J. Piano 
et al, Phys. Rev. D32,2441 (1985). 

T. Kitagaki et al. [ E-745 Collaboration ], Phys. Lett. B214,281 (1988). 
J. Guy et al. [ BEBC Collaboration ], Phys. Lett. B229,421 (1989). 

H. Abramowicz,T. Hansl-Kozanecka,J. May,P. Palazzi,A. Para,F. Ranjard,A. Savoy- 
Navarro,D. Schlatter et al, Z. Phys. C25,29 (1984). 



M. Tzanov et al.,PRD 74 (2006) 012008 

G. Onengut et &L,Phys. Lett. B 632 (2006) 65 

J. Altegoer et al. [NOMAD Collaboration], Nucl. Instrum. Meth. A 404,96 (1998). 
R. Petti [NOMAD Collaboration], hep-ex/0411032. 

G. T. Jones et al. [Birmingham-CERN-Imperial College-Munich(MPI)-Oxford- Collaboration], Z. 
Phys. C 62,601 (1994). 

BCDMS Collaboration: A. C. Benvenuti et al, Phys. Lett. B 223 (1989) 485; BCDMS Collabo- 
ration: A. C. Benvenuti et al, Phys. Lett. B 236 (1989) 592. 

HI Collaboration: C. Adloff et al, Eur. Phys. J. C 13 (2000) 609 [hep-ex/9908059]; Eur. Phys. 
J. C 19 (2001) 269 [hep-ex/0012052]; Eur. Phys. J. C 21 (2001) 33 [hep-ex/00 12053]. 

ZEUS Collaboration: S. Chekanov et al, Eur. Phys. J. C 21 (2001) 443 [hep-ex/0105090]; A.M. 
Cooper-Sarkar, Proceedings of International Europhysics Conference on HEP 2001, Budapest [hep- 
ph/0110386]. 

CDF Collaboration: T. Affolder et al, Phys. Rev. D 64 (2001) 032001 [hep-ph/0102074]. 

D0 Collaboration: B. Abbott et al, Phys. Rev. Lett. 86 (2001) 1707 [hep-ex/0011036]; and Phys. 
Rev. D 64 (2001) 032003 [hep-ex/0012046]. 

CDF Collaboration: F. Abe et al, Phys. Rev. Lett. 81 (1998) 5754 [hep-ex/9809001]. 

E866 Collaboration: R. S. Towell et al, Phys. Rev. D 64 (2001) 052002 [hep-ex/0103030]. 

E605 Collaboration: G. Moreno et al, Phys. Rev. D 43 (1991) 2815. 

A. Bruell, private communication. 

D.M. Aide et &l,Phys. Rev. Lett. 64 (1990) 2479 . 

S. A. Kulagin and R. Petti, Nucl. Phys. A 765,126 (2006) [hep-ph/0412425]. 



74 



[200 
[201 
[202 
[203 
[204 
[205 
[206 

[207; 
[208 

[209 
[210 
[211 

[212 
[213 

[214 

[215 
[216 

[217 

[218 

[219 

[220 

[221 

[222 

[223 

[224 

[225 
[226 



S. A. Kulagin and R. Petti, AIP Conf. Proc. 967,94 (2007) [arXiv:0711.3956 [nucl-th]]. 

S. A. Kulagin and R. Petti, Phys. Rev. D 76,094023 (2007) [hep-ph/0703033 [HEP-PH]]. 

M. Gluck and E. Reya, Nucl. Phys. B 130,76 (1977). 

D. W. Duke and J. F. Owens, Phys. Rev. D 30,49 (1984). 

A. D. Martin,R. G. Roberts and W. J. Stirling, DTP/87/34. 

J. G. Morfin,W. -K. Tung, Z. Phys. C52, 13-30 (1991). 

W. K. Tung,J. G. Morfin,H. Schellman,S. Kunori,A. Caldwell and F. I. Olness, Structure Functions 
and Parton Distributions," FERMILAB-CONF-89-026 

M. Guzzi,P. NadolskyE. Berger,H. L. Lai,F. Olness and C. P. Yuan, arXiv:1101.0561 [hep-ph]. 

R. S. Thorne,A. D. Martin,W. J. Stirling and G. Watt, PoSDIS 2010,052 (2010) [arXiv: 1006.2753 
[hep-ph]]. 

A. Guffanti [NNPDF Collaboration], AIP Conf. Proc. 1369,21 (2011). 
F. James and M. Roos, Comput. Phys. Commun. 10,343 (1975). 

R. Thorne and R. Roberts, Phys. Lett. B421 (1998) 303. A. D. Martin et. al. Eur. Phys. J. C18 
(2000) 117. 

L. W. Whitlow et. a/,Phys. Lett. B250 (1990) 193. 

K. J. Eskola,V. J. Kolhinen,and C. A. Salgado,T7ie scale dependent nuclear effects in parton 
distributions for practical applications, Eur. Phys. J. C9 (1999) 61-68, 

K. J. Eskola,V. J. Kolhinen,H. Paukkunen,and C. A. Salgado,j4 global reanalysis of nuclear parton 
distribution functions, JEEP 05 (2007) 002, 

H. Paukkunen, arXiv:0906.2529 [hep-ph]. 

K. J. Eskola,H. Paukkunen and C. A. Salgado, Nucl. Phys. A 855,150 (2011) [arXiv:1011.6534 
[hep-ph]]. 

M. Hirai,S. Kumano,and M. Miyama, Determination of nuclear parton distributions, Phys. Rev. 
D64 (2001) 034003, 

M. Hirai,S. Kumano,and T. H. Nagai, Nuclear parton distribution functions and their uncertainties, 
Phys. Rev. C70 (2004) 044905, 

M. Hirai,S. Kumano,and T. H. Nagai, Determination of nuclear parton distribution functions and 
their uncertainties in next-to-leading order, 

D. de Florian and R. Sassot, Nuclear parton distributions at next-to-leading order, Phys. Rev. D69 
(2004) 074028, 

K. Kovarik,I. Schienbein,F. I. Olness, J. Y. Yu,C. Keppel,J. G. Morfm,J. F. Owens and T. Stavreva, 
Phys. Rev. Lett. 106,122301 (2011) [arXiv: 1012.0286 [hep-ph]]. 

I. Schienbein,J. Y. Yu,K. Kovarik,C. Keppel,J. G. Morfm,F. Olness and J. F. Owens, Phys. Rev. 
D 80,094004 (2009) [arXiv:0907.2357 [hep-ph]]. 

I. Schienbein,J. Y. Yu,C. Keppel,J. G. Morfm,F. Olness and J. F. Owens, Nucl. Phys. Proc. Suppl. 
191,25 (2009) [arXiv:0812.3370 [hep-ph]]. 



I. Schienbein,J. Y. Yu,C. Keppel,J. G. Morfin,F. Olness and J. F. Owens, Phys. Rev. D 77,054013 
(2008) [arXiv:0710.4897 [hep-ph]]. 

D. de Florian, R. Sassot, P. Zurita and M. Stratmann, arXiv:1112.6324 [hep-ph]. 

NuTeV Collaboration, M. Goncharov et al., Precise Measurement of Dimuon Production Cross- 
Sections in v^Fe andv^Fe Deep Inelastic Scattering at the Tevatron, Phys. Rev. D64 (2001) 112006, 



75 



[227] http : // www . praktika . physik . uni-bonn . de/ conf erences/dis-2012?set_language=en] 

[228] G. Berlad,A. Dar and G. Eilam, Phys. Rev. D 22,1547 (1980). 

[229] K. J. Eskola,V. J. Kolhinen and P. V. Ruuskanen, Nucl. Phys. B 535,351 (1998) [hep-ph/9802350]. 

[230] H. Paukkunen and C. A. Salgado, JHEP 1007,032 (2010) [arXiv:1004.3140 [hep-ph]]. 

[231] H. Paukkunen and C. A. Salgado, PoS ICHEP 2010,166 (2010) [arXiv:1009.3143 [hep-ph]]. 

[232] K. J. Eskola,H. Paukkunen and C. A. Salgado, JHEP 0904,065 (2009) [arXiv:0902.4154 [hep-ph]]. 

[233] M. Hirai,S. Kumano and K. Saito, AIP Conf. Proc. 1189,269 (2009) [arXiv:0909.2329 [hep-ph]]. 

[234] M. Gliick,E. Reya,A. Vogt, Eur. Phys. J. C5,461 (1998). 

[235] A. D. Martin,W. J. Stirling,R. S. Thorne,and G. Watt, Eur. Phys. J. C63,189 (2009). 

[236] J. Botts et al. CTEQ parton distributions and flavor dependence of sea quarks. Phys. 
M.,B304:159-166,1993,hep-ph/9303255. 

[237] H. L. Lai et al. Global QCD analysis and the CTEQ parton distributions. Phys. i?ei>.,D51:4763- 
4782,1995,hep-ph/9410404. 

[238] S. J. Brodskyl. Schmidt, and J. J. Yang. Nuclear antishadowing in neutrino deep inelastic scat- 
tering. Phys. J Reu,D70:116003,2004,hep-ph/0409279. 

[239] J. -W. Qiu and I. Vitev, Phys. Lett. B 587,52 (2004) [hep-ph/0401062]. 

[240] D. Drakoulakos et al. Proposal to perform a high-statistics neutrino scattering experiment using 
a fine-grained detector in the NuMI beam, hep-ex/0405002. 

[241] V. GuzeyL. Zhu,C. Keppel,M. E. ChristyD. Gaskell,P. Solvignon and A. Accardi, "Impact of 
nuclear dependence of R=<tl/<jt on antishadowing in nuclear structure functions," arXiv:1207.0131 
[hep-ph] . 

[242] V. Guzey et al., "Nuclear shadowing in charged-lepton-nucleus and neutrino-nucleus scattering," 
work in progress. 



76 



